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Abstract 

Supernova remnants are beautiful astronomical objects that are also of high scientific interest, 
because they provide insights into supernova explosion mechanisms, and because they are the 
likely sources of Galactic cosmic rays. X-ray observations are an important means to study these 
objects. And in particular the advances made in X-ray imaging spectroscopy over the last two 
decades has greatly increased our knowledge about supernova remnants. It has made it possible 
to map the products of fresh nucleosynthesis, and resulted in the identification of regions near 
shock fronts that emit X-ray synchrotron radiation. Since X-ray synchrotron radiation requires 
10-100 TeV electrons, which lose their energies rapidly, the study of X-ray synchrotron radiation 
has revealed those regions where active and rapid particle acceleration is taking place. 

In this text all the relevant aspects of X-ray emission from supernova remnants are reviewed 
and put into the context of supernova explosion properties and the physics and evolution of su- 
pernova remnants. The first half of this review has a more tutorial style and discusses the basics 
of supernova remnant physics and X-ray spectroscopy of the hot plasmas they contain. This in- 
cludes hydrodynamics, shock heating, thermal conduction, radiation processes, non-equilibrium 
ionization. He-like ion triplet lines, and cosmic ray acceleration. The second half offers a review 
of the advances made in field of X-ray spectroscopy of supernova remnants during the last 15 year. 
This period coincides with the availability of X-ray imaging spectrometers. In addition, I discuss 
the results of high resolution X-ray spectroscopy with the Chandra and XMM-Newton gratings. 
Although these instruments are not ideal for studying extended sources, they nevertheless pro- 
vided interesting results for a limited number of remnants. These results provide a glimpse of 
what may be achieved with future microcalorimeters that will be available on board future X-ray 
observatories. 

In discussing the results of the last fifteen years I have chosen to discuss a few topics that 
are of particular interest. These include the properties of Type la supernova remnants, which 
appear to be regularly shaped and have stratified ejecta, in contrast to core collapse supernova 
remnants, which have patchy ejecta distributions. For core collapse supernova remnants I discuss 
the spatial distribution of fresh nucleosynthesis products, but also their properties in connection 
to the neutron stars they contain. 

For the mature supernova remnants I focus on the prototypal supernova remnants Vela and the 
Cygnus Loop. And I discuss the interesting class of mixed-morphology remnants. Many of these 
mature supernova remnants contain still plasma with enhanced ejecta abundances. Over the last 
five years it has also become clear that many mixed-morphology remnants contain plasma that 
is overionized. This is in contrast to most other supernova remnants, which contain underionized 
plasmas. 

This text ends with a review of X-ray synchrotron radiation from shock regions, which has 
made it clear that some form of magnetic-field amplification is operating near shocks, and is an 
indication of efficient cosmic-ray acceleration. 
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1 Introduction 

Supernovae play a central role in modern astrophysics. 
They are of prime importance for the chemical evolution 
of the Universe, and they are one of the most important 
sources of energy for the interstellar medium (ISM). Part 
of that energy is in the form of cosmic rays, which have 
an energy density in the Galaxy of 1-2 eV cm^^, thus con- 
stituting about one third of the total energy density of the 
ISM. Finally, supernovae, in particular Type la supernovae 
(see Sect. |2]i, play a central role in present day cosmology 
as their brightness al lows them to be d etected at high red- 
shifts (up to z ^ 1.7. lRiess et al. , 2007 ). Their use has led 
to the recognition that the expansi on of the Universe is ac- 



celerating instead of d ecelerating (IPerlmutter et al.L 11998 



Garnavichetal.. 1998). 



Supernova remnants (SNRs) are an important means 
to study supernovae. Since supernovae are relatively rare 
(2-3 per century in a typical spiral galaxy like our own), 
SNRs provide the best way to study the local population 
of supernovae. In addition, SNRs can reveal details about 
the explosions that are difficult to obtain from studying su- 
pernovae directly. For example, young SNRs can inform 
us about nucleosynthesis yields of individual SNRs and 
about the inherent asymmetries of the explosion itself, as 
revealed by the spatial and velocity distribution of heavy 
elements. Moreover, SNRs probe the immediate surround- 
ings of supernovae, which are shaped by their progenitors. 

Another aspect of SNRs concerns the interesting physi- 
cal process that shape their properties. SNR shocks pro- 
vide the best Galactic examples of high Mach number, 
coUisionless shocks. The physics of these shocks is not 
well understood, as heating of the atoms occurs collision- 
less, i.e. shock heating does not operate through particle- 
particle (Coulomb) interactions, but through electro- 
magnetic fluctuations and plasma waves. The SNR shocks 
are also thought to be the locations where part of the ex- 
plosion energy is converted to cosmic-ray energy. This is 
supported by observations that indicate that GeV and TeV 
particles are present in SNRs through their radiative sig- 
natures: radio to X-ray synchrotron emission from rela- 
tivistic electrons, and 7-ray emission caused by acceler- 
ated electrons and ions. 

X-ray observations are one of the most important means 
to study the many aspects of SNRs. X-ray spectroscopy 
is essential to obtain abundances of the prime nucleosyn- 
thesis products of supernovae, which are the so-called 
alpha-elements (O, Ne, Mg, Si, S, Ar, Ca) and iron-group 
elements (chiefly Fe, Ni, and some trace elements with 
20 < Z < 28). All these elements have prominent emis- 
sion lines in the 0.5-10 keV band for temperatures be- 
tween 0.2-5 keV, which happens to be the typical electron 
temperatures of plasma heated by SNR shocks. The hot 
plasmas of SNRs are also optically thin, so inferring abun- 
dances is relatively straightforward. X-ray spectroscopy of 
SNRs, therefore, provide us with a record of alpha- and 
iron-group element production by supernovae. 

X-ray spectroscopy can also be used to study several as- 



pects of SNR shock physics. First of all the line emission 
provides information about the temperature and ionization 
state of the plasma, but the absence of lines or weak line 
emission in young SNRs usually points toward the con- 
tribution of X-ray synchrotron radiation. Studying X-ray 
synchrotron radiation offers a powerful diagnostic tool to 
study electron cosmic -ray acceleration, as the presence 
of X-ray synchrotron radiation depends sensitively on the 
shock acceleration properties. Moreover, the size of the 
X-ray synchrotron emitting regions can be used to infer 
magnetic-field strengths. 

One of the major advances in X-ray spectroscopy over 
the last 15-20 years has been the emergence of X-ray imag- 
ing spectroscopy. Although presently the spectral resolu- 
tion offered is limited, it provides a direct way to map the 
spatial distribution of elements and temperatures in SNRs. 
For young SNRs it also helps to separate thermal from 
non-thermal (synchrotron) X-ray emission. 

The power of X-ray imaging spectroscopy is illustrated 
in Fig. [T] which shows Tycho's SNR (the remnant of the 
historical SN 1572) as observed by Chandra. The im- 
age shows beautifully the distribution of silicon and iron, 
but also reveals the synchrotron dominated emission from 
near the shock front. 

Over the last ten years high spectral resolution X-ray 
spectroscopy has gained in importance, due to the avail- 
ability of grating spectrometers on board Chandra and 
XMM-Newton. These spectrometers are also used to study 




Figure 1: Three-color Chandra image of Tycho's SNR 
(SN1572/G120.1H-1.4). The color red shows Fe L-shell 
emission, green Si XIII, and blue continuum (4-6 keV). 
Note the very narrow continuum rims near the shock 
front (blueish/purple in this image), caused by X-ray syn- 
chrotron radiation from electrons with energies up to 
100 TeV. (Image made by the author from Chandra data, 
see alsO iHwang et al.> .2002:, Warren et al.,.2005.) . 
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Figure 2: The classification of supernovae, based on opti- 
cal spectroscopy and light-curve shape. 

SNRs, but the contribution of grating spectrometers to our 
understanding of SNRs has been limited. The reason is 
that the slitless grating spectrometers are not particularly 
well suited for extended sources such as SNRs. Neverthe- 
less, some important studies of SNRs with properties that 
make them suitable for grating spectroscopy have been 
published. With the emergence of high spectral resolution 
imaging spectroscopy, as will be possibl e with the calori 
metric spectrometers on board Astro-H ("Takahas hi et al 



I2OO8) and IXO/ATHENA (Bookbinder, 2010), high reso- 
lution X-ray spectroscopy will become much more promi- 
nent. 

Here I review the current status of X-ray spectroscopy 
of SNRs. Since there is a lack of good text books or ex- 
tended reviews on SNRs in general, the first part of this 
review has a tutorial character, describing our current un- 
derstanding of the evolution and physics of SNRs. The 
second half of the review focusses on the X-ray observa- 
tions themselves, with an emphasis on the achievements 
obtained with data from the three most recent X-ray ob- 
servatories: Chandra, XMM-Newton, and Suzaku. 



2 Supernovae 

Supernovae are divided into two broad categories, reflect- 
ing our understanding of the explosion processes: core col- 
lapse supernovae and thermonuclear supernovae. In addi- 
tion an observational classification scheme is used (Fig.|2]l 
that goes back to Minkowski (1941), who observed that 
some supernovae do not show hydrogen absorption in 
their spectra (Type I) and others do (Type II). 

Type II supernovae are invariably core collapse super- 
novae, but Type I supernovae can be either core collapse or 
thermonuclear supernovae. The thermonuclear explosions 
are associated with spect roscopic class Type la (Elias 
et al. nT985l) . which have Si absorption lines in their spec- 
tra. Type lb and Ic supernovae are now understood to be 
explosions of stars that have lost their hydrog en-rich enve 



lope as a result of stellar wind mass loss (e.g. lHeger et al 



2003h. or through binary interaction (e.g. Podsiadlowski 
et al., 1992h . For Type Ic the mass loss seems to have re- 
moved even the helium-rich layers of the progenitor The 
fact that they are preferentially found in the most luminous 
regions of galaxies suggests t hat they are the e xplosions 
of the most massive stars (e.g. lKellv et al.L 120081) . Interest- 
ingly, also long duration gam ma-ray bursts are ass ociated 
with Type Ic supernovae (e.g. Galama et al.L 19981: Stanek 



et al., '20031). 

The sub-division of the Type II class in Type IIP 
(plateau). Type IIL (linear light curve) and Type lib is 
based on a combination of two observational criteria: op- 
tical spectroscopy and light-curve shape ( Fig. |2]i. Type 
IIP are the most common type of core collapse super- 
novae, and optical studies of potential progenitor stars 
confirm that their progenitors have initial masses in the 
^8 — 17 M0 range, and that they explode in the red su- 
pergiant phase, while still having a substantial hydrogen 
envel ope (see the discussions in Smartt . 2009; Chevalier , 
20051). Type IIL progenitors probably have a substantially 
less massive envelope, either due to stellar wind mass loss, 
or due to binary interaction. Type lib supernovae are a 
class intermediate between Type lb and Type II, in that 
their spectra would initially identify them as Type II ex- 
plosions, but at late times their spectra evolve into Type lb 
spectra. Also this can be understood as the result of sub- 
stantial, but not complete, removal of the hydrogen-rich 
envelope due to stellar wind mass loss, or binary inter- 
action. The prototypal Type lib supernova is S N 1993 J 



( Podsiadlowski etal., 19931: Woosley et al., 1994 ). Interest- 



ingly, recent identification and subsequent spectroscopy of 
the light echo of the supernova that caused the extensively 
studied SNR Cassiopeia A (Cas A) shows that it is the 
remnant of a Type lib supernova, as the spectrum shows 
both hydrogen and weak helium line absorption ( Krause 
et al.. l2008al) . 

Not listed in Fig. |2]is the Type Iln class. Type Iln su- 
pernovae are characterized by narrow hydrogen emission 
lines, which are thought to come from a dense circum- 
stellar environment, probably caused by substantial mass 
lost by the progenitor. Its place in the diagram is not quite 
clear, as at least one Type Iln supernova, SN 2001ic, was 
observed to be a Type la supernova whose s pectrum subse 



quent ly evolved into a Type Iln supernova (iHamuv et al 
2003h. 



2.1 Core collapse supernovae 

Core collapse supernovae mark the end of the lives of 
massive stars; that is, t hose stars with main seque nce 



masses M > 8 Mq (see IWooslev and Jankal I2005L for 



'But see lVinkI <2006l) ; lRevnoldsl )2008l) ; lBade"ne3 <2010l) for recent, 
more topical reviews. 



a review). Just prior to collapse the star consists of dif- 
ferent layers with the products of the different consecu- 
tive burning stages. From the core to the outside one ex- 
pects: iron-group elements in the core (silicon-burning 
products), then silicon-group elements (oxygen-burning 
products), oxygen (a neon burning product), neon and 
magnesium (carbon-burning products), carbon (a helium- 
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Figure 3: Left: Supernova yields for the most abundant X-ray emitting elements. The squares/black line indicates the 
mean yield for core collapse supernovae, whereas the circles indicate thermonuclear supernovae (the W7 deflagration 
model in red and the WDD2 delayed detonation model in magenta). The model yields were taken from llwamoto et al. 
(Il999h . Right: O xygen yield of core collapse s upernovae as a function of mai n sequence mass. The circles and squares are 
th e predictions of | Woosley and Weaved i the triangles are predictions o fI Chieffi and Limongil (l2004f) . and the crosses 

of Thielemann et al. I (ll996h . In general the oxygen yields obtained by the various groups are very similar, but above 30 M© 
one sees that certain models ( Wooslev and Weavei , 19951 using 10^' erg explosion energies) predict a diminishing oxygen 
yield. The reason is that above 30 M© stellar cores may collapse into black holes, and part of the oxygen falls onto the 
black hole. The amount of fall-back is governed by the explosion energy and the amount of pre-supernova mass loss, but 
it is also sensitive to the numerical treatment of the explosion. 



burning product), helium (a hydrogen-burning product), 
and, finally, unprocessed hydrogen-rich material. 

The creation of the iron-group core, which lasts about a 
day, is the beginning of the end of the star, as no energy 
can be gained from nuclear fusion of iron. The core col- 
lapses into a proto-neutron star, and for the most massive 
stars into a black hole. Most of the gravitational energy 
liberated {E ^ GM^/Ras ^ 10^^ erg, with R„, the neutron 
star radius) is in the form of neutrinos. This has been con- 
firmed with the d etection of n eutrinos from SN1987A by 
the Kamiokandej Hirata et aO , 1^87) and Irvine-Michigan- 
Brookhaven (iBionta et al r fl987h water Cherenkov neu- 
trino detectors. 

The supernova explosion mechanism itself, which re- 
quires that > 10^' erg of energy is deposited in the outer 
layers, is not well understood. The formation of a proto- 
neutron star suddenly terminates the collapse, and drives 
a shock wave through the infalling material. However, nu- 
merical simulations show that the shock wave stalls. It is 
thought that the shock wave may be energized by the ab- 
sorption of a fraction of the neutrinos escaping the proto- 
neutron star, but most numerical models involving neu- 
trino absorption are still unsuccessful in reproducing a su- 
pernova explosion (fjanka et al. , 2007 ). The most recent 
research, therefore, focusses on the role of accretion insta- 
bilities for the supernova explosion, as these instabilities 
help to enhance the neutrino absorption in certain regions 
just outside the proto-neutron star. A promising instabil- 
ity is the so called, non- spherical ly symmetric s t andins 



accretion shock instability (SASI, iBlondin et al 



2003) 



tions bv lBurrows et al.l (l2007h suggest that acoustic power, 
generated in the proto-neutron star due to g-mode oscilla- 
tions, ultimately leads to a suc cessful explosion (but see 
[Weinberg and Ouat aertL 2008 1). Finally, there have been 
suggestions that neutrino deposition is not the most impor- 
tant ingredient for a successful explosion, but that ampli- 
fication of the stellar magnetic field, due to differential ro- 
tation and compression, may lead to mag neto-centrifugal 



jet fo rmation, which drives the explosion dWheeler et al. 



2002) 



SASI, acoustic power, and magneto-centrifugal mod- 
els all predict deviations from spherical symmetry. In 
the magneto-rotational models one even expects a bipolar 
symmetry. Another reason why deviations from spherical 
symmetry has received more attention is that long duration 
gamma-ray bursts are associated with very energetic Type 
Ic supernovae (hypernovae). Given the nature of gamma- 
ray burst these explosions are likely jet driven. This raises 
the possibility that also more nor mal core c ollapse super- 
novae have jet components (Whee ler et al.l 2002). There 
is indeed evidence, based on optical polarimetry, that core 
colla pse supernovae, esp ecially Type Ib/c, are aspherical 



lap 



(e.g. IWang et all l2001h . In Sect. 19.11 1 will discuss evi 



which may also help to explain p ulsar kicks and rotation 
(IBlondin and MezzacappaL l2007h . Alternatively, simula- 



dence that also the remnants of core collapse supernovae 
show deviations from spherical symmetry. 

The ejecta of core collapse supernovae consist primarily 
of stellar material, except for the innermost ejecta, which 
consist of explosive nucleosynthesis products, mostly Fe 
and Si-group elements. These elements are synthesized 
from protons and alpha-particles, which are the remains 
of the heavy elements that have disintegrated in the in- 
tense heat in the the innermost regions surrounding to the 
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collapsing core ( Arnetl . 1996h . Some of the explosive nu- 



cleosynthesis products are radioactive, such as ^''Ni, and 
^'^Ti. In particular the energy generated by the decay of 
^^Ni into ^^Co, and finally ^^Fe, heats the ejecta, which 
leaves a major imprint on the evolution of the supernova 
light curve. The yields of these elements depend sensi- 
tively on the details of the explosion, such as the mass 
cut (the boundary between material that accretes on the 
neutron star and material that is ejected), explosion en- 
ergy, and explosion asymmetry. Since the mass of the neu- 
tron star/black hole, the location of energy deposition and 
the presence of asymmetries are not well constrained, the 
expected yields of these elements are uncertain, and vary 
substantially from one set of models to the other (F ig. [3] 
Wooslev an d Weaver, 1995^ Thielemann et al. ^ 19961: Chi- 
effi and Limongi, 2004). 

Overall the yields of core collapse supernovae are dom- 
inated by carbon, oxygen, neon and magnesium, which 
are products of the various stellar burning phases (e.g. 
Wooslev and Weave^ 19951: Thielemann et al.l 19961: Chi- 
effi and Limongi, '2004V These yields are a function of the 
initial mass of the progenitor (Fig.|3]l. It is for this reason 
that oxygen-rich SNRs (Sect. 19. il l are considered to be the 
remnants of the most massive stars. 



2.2 Thermonuclear supernovae 

Type la supernovae are generally thought to be thermonu- 
clear explosions of C/O white dwarfs, i.e. the explosion 
energy originates from explosive nuclear burning, rather 
than from gravitational energy liberated during the col- 
lapse of a stellar core. 

Although there is some variation in peak brightness of 
Type la supernovae, the variation is much less than that of 
Type II supernovae. This is in line with the idea that all 
Type la supernovae are explosions of similar objects: C/O 
white dwarfs with masses close to the Chandrasekhar limit 
(1 .38 M0). Moreover, an empirical relation exists between 
their peak b rightness and the post peak decline rate of the 
light curve dPhillips et al.L Il992h . which can be used to 
calibrate the absolute peak brightness of each event. This 
makes Type la supernovae excellent distance indicators, 
on which much of the evidence rests tha t the expansi on of 
the Univers e is acceleratin g (Perlmutte r et al.L 1998 : Gar- 
navich et al.. ll998l:lRie'sset al.L 120071) . 

There is no direct observational evidence that Type la 
progenitors are white dwarfs, but the fact that only Type la 
supernovae can occur among old stellar populations indi- 
cates that massive stars cannot be their progenitors. Their 
relative uniformity can be best explained b y assuming a 



proge nitor type with a narrow mass range dMazzah et al 



20071) . Moreover, C/O white dwarfs close to the Chan- 
drasekhar mass limit are very likely Type la progenitors, 
since their high den sity makes for an ideal "nuclear fu- 
sion bomb" (lArnettill996) : Once a nuclear reaction in the 
core is triggered, it will result in an explosion. The trigger 
mechanism itself is, however, not well understood. 
A more serious problem is that we do not know what 



the most likely progenitor systems are. It is clear that C/O 
white dwarfs have to accrete matter in order to reach the 
Chandrasekhar limit, so thermonuclear supernovae must 
occur in a binary system. However, this still leaves several 
possibiUties for the progenitor systems: double degener- 
ate systems (two white dwarfs), or white dwarfs with ei- 
ther a main sequence star, or an evolved companion (i.e. 
single degenerate). Only a limited range of mass trans- 
fer rates, 4 x 10"** - 7 x 10~^ Mpyr~ \ lead to stable 
gro wth of the white dwarf jNomotol 19821: Shen and Bild- 



sten, 



2007h . Novae, for instance, are relatively slow accre- 



tors, which are thought to blow off more mass than they 
accrete. The only stable white dwarf accretors seem to be 
supersoft sourcefl but their population seem s too small to 



account for the observed supernova rate (e.g. lRuiter et al 



12009) . Depending on the progenitor evolution and accre 
tion scenarios, the progenitor binary system may affect its 
circumstellar medium. For example, in the case of wind 
accretion, not all the mass lost from the donor w ill end up 
on the white dwarf. iHachisu et al.l (11996 , Il999h proposed 
a scenario in which Roche-lobe overflow from the donor 
star is stabilized through a fast, optically thick, wind from 
the white dwarf, which will also affect the immediate sur- 
roundings of the SNR. 

For some time it was thought that Type la supernovae 
were associated with population II stars, in which case 
they must evolve on time scales > 10^ yr, but it has re- 
cently been established that also a short "channel" exists, 
with an evolutionarv time scale of ^ 10^ vr (e.g. Mannucci 
et al.,|2006|;|AubourgetalJ,|2008|). 

Models for thermonuclear explosions come in three 
classes: detonation, deflagration, and delayed detonation 
models. In detonation models the explosive nucleosynthe- 
sis occurs due the compression and heating of the plasma 
by a shock wave moving through the star In deflagration 
models the burning front proceeds slower than the local 
sound speed. The nuclear fusion in the burning front is 
sustained by convective motions that mixes unburnt mate- 
rial into the hot burning zone. The classical deflagration 
model is the W7 of Nomoto et al. ( 1984). 

Pure detonation models predict that almost all of the 
white dwarf matter will be transformed into iron-group 
elements, whereas optical spectroscopy of Type la super- 
novae show that the ejecta c ontain significant am ount of 
intermediate mass elements (^ Branch et al. . 1982 ). On the 
other hand, the pure deflagration models overpredict the 
production of ^^Fe with respect to ^^Fe and predict too 
narrow a velocity range for the intermediate mass ele- 
ments compared to observations. For these reasons, the 
currently most popular model for Type la supernovae are 
the delayed detonation (DDT) models (lKhokhlovill99lh . 
in which the explosion starts as a deflagration, but changes 
to a detonation wave burning the remainder of the white 
dwarf into intermediate mass elements (IMEs), such as sil- 
icon. There is observational evidence that the fraction of 



"Soft X-ray point sources that are th ought to be white dwarfs that 
stably accrete material from a companion jKahabka and van den Heuvell 
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C/O that is burned is roughly constant M = 1.1 M0, but 
that the variation in peak brightness is either caused by 
the ratio of iron-group elements over IME products (Maz- 



zali et al.. ,2007,). or b v the ratio of stable iron over ^''Ni 



JWooslev et al.Ll2007h . ^""Ni decays into ^^Fe, and the heat 
that this generates determines the brightness of a Type la 
supernova (see Sect. l6.2l i. 

Figure [3] shows the predicted overall abundance pat- 
tern of deflagration and DDT models. Compared to core 
collapse supernovae, thermonuclear supernovae produce 
much more iron-group elements 0.6 Mp ). The explo- 
sion energy, E, of thermonuclear explosio ns is determined 



by th e mass of their burning products dWooslev et al 
2007h : 



£51 = 1 .56MNi + 1 .74MFe + 1 .24MiME - 0.46, ( 1 ) 

with the £51 the final kinetic energy in units of 10^' erg, 
and the masses of stable Fe, ^^Ni and IME in solar units. 



3 The classification of supernova 
remnants 

Given the fact that supernovae can be broadly classified in 
core collapse and thermonuclear/Type la supernovae, one 
would hope that SNRs would be classified as core collapse 
or Type la SNRs. Although it is possible now to determine 
the supernova origin of a given SNR using X-ray spec- 
troscopy (see Sect. [T] and |8]l, for many old SNRs, whose 
emission is mainly coming from swept-up gas, one has to 
rely on secondary indicators for determining their super- 
nova origin. The most reliable indicator is of course the 
presence of a neutron star inside the SNR, which makes 
it clear that the SNR must have a core collapse origin. 
But ev e n here one has to be aware of chance alignments 
(IKaspil [19981), in particular for SNRs with large angular 
diameters. A secondary indicator for a SNR to originate 
from a core collapse event is whether a SNR is located in a 
star forming reg ion or inside an OB association (e .g West- 
erlund, 1969). But this does not constitute proof for such 
an origin. In contrast, a position of a SNR high above the 
Galactic plane can be taken as supporting evidence for a 
Type la origin. Such is the cas e for, for example, SN 1006 
(IStephenson and GreenL 120021) . 

Because the supernova origin of SNRs is often diffi- 
cult to establish, SNRs have a classification of their own, 
which is mostly based on their morphology. Tradition- 
ally, this classification recognized three classes: shell type 
SNRs, plerions, and composite SNRs. As the blast wave 
sweeps through the interstellar medium (Sect. |4]) a shell 
of shock heated plasma is created. Therefore, in many 
cases the morphology of a SNR is characterized by a limb 
brightened shell, which classifies the SNR as a shell-type 
SNR. 

However, in case of a core collapse one may expect the 
presence of a rapidly rotating neutron star It loses energy 
with a rate of £ = IClCl = An^IP/P^, with £1, the angu- 
lar frequency, P the rotational period and / « 10^^ g cm^ 



the moment of inertia. This energy loss produces a wind 
of relativistic electrons and positrons, which terminates in 
a shock, where the electrons/positrons are accelerated to 
ultra-relativistic energies. These particles advect and dif- 
fuse away from the shock creating a nebula of relativistic 
electron/positrons which emit synchrotron radiation from 
the radio to the soft 7-ray bands, and inverse Compton 
scatt ering from soft y-rav to the TeV band ( Gaensler and 
Slane, 2006, for a review). Such a nebula is aptly named 
a pulsar wind nebula. The most famous pulsar wind neb- 
ula is the Crab Nebula (also known as Ml or Taurus A), 
which is powered by the pulsar B0531H-21. This nebula is 
associated wit h the historical supernova of 1054 ( Stephen- 
son and Green, l2002h . As a result the Crab Nebula, but 
also similar objects like 3C58, have the status of a SNR. 
Since the nebula has a morphology that is bright in the cen- 
ter, and does not show a shell, they are called ^Weii center 
SNRs or plerions, with the name plerion derived from th e 
Greek word for "full", pleres ( Weiler and Panagiai 1978 ). 

It is arbitrary whether one should call the Crab-nebula 
and related objects plerions or pulsar wind nebula. The 
radio and X-ray emission from plerions is powered by 
the pulsar wind, not by the supernova explosion, and as 
such the nebula should not really be referred to as a SNR. 
To complicate matters, the Crab nebula does show opti- 
cal line emission from supernova ejecta, which could be 
referred to as a SNR. Pulsar wind nebulae can be found 
around both young and old pulsars, even a few millisec- 
ond pulsar s have small nebulae around them (" Kargaltsev 
and Pavlov, 2008h . Since these older pulsars do not have 
a connection with recent supernova activity, one should 
not call them plerions. So the name pulsar wind nebula is 
more generic and informative than plerion, and is, there- 
fore, preferable. 

Energetic pulsars with ages less than 20,000 yr are 
expected to have blown a pulsar nebula while they are still 
surrounded by the SNR shell. One expects then a radio and 
X-ray morphology that consists of a pulsar wind nebula 
surrounded by a shell. Indeed a few SNRs have this char- 
acteristic (Fig. HI;) and are classified as composite SNRs. In 
fact, it is still puzzling why a young o bject like the C rab 
Nebula does not show a SNR shell (see iHesten, l2008l for 
a discussion). 

Since the 1980s the SNR classification has become 
broader Due to the imaging capabilities of X-ray satel- 
lites like Einstein and ROSATit became clear that there are 
many SNRs that display a shell-type morphology in radio 
emission, but whose X-ray emission mainly comes from 
the center of the SNR ( White and Long, 1 99 1 ; Rho and Pe- 
tre. ll998l) . as illustrated in Fig.gli. In general, these SNRs 
are relatively old, and are associated with dense interstel- 
lar medium. The X-ray emission from the center of these 
SNRs was not powered by a pulsar, but consisted of ther- 
mal emission from a hot plasma. These SNRs are referred 
to as mixed-morphology SN Rs (Rho and Petre, 1998) or 
as thermal-composite SNRs jshelton et all Il999h . I will 
discuss their properties in Sect. 110. 3] 

X-ray spectroscopy has greatly enhanced our ability 
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Figure 4: The SNR morphological classification illustrated with examples. From top left to bottom rig ht: a) T he Cygnus 
Loop, a shell-type SNR with a diameter of 3°, as observed by the ROSAT PSPC instrument (Levens on et al l[i998), red 
is very soft emission from ^0.1-0.4 keV, green ^ 0.5 — 1.2 keV, and blue ^ 1.2 — 2.2 keV. b) 3C58, a plerion/pulsar 
wind nebula, as observed by Chandra (Slane et al.|, 2004 ). The long axis of this object is ^7'. c) The composite SNR Kes 
75 as observed by Chandra (Helfand et al.. 2003i ) with the inner pulsar wind nebula, which has a hard X-ray spectrum, 
powered by the pulsar Jl 846-0258. The partial shell has a radius of ^1.4'. The colors indicate 1-1.7 keV (red, Ne and Mg 
emission), 1.7-2.5 keV (green, Si/S), and 2.5-5 keV (blue, mostly continuum emission from the pulsar wind nebula), d) 
The "thernial-composite" SNR W28 as observed in X-rays by the ROSAT PSPC (blue) and in radio by the VLA ( Dubner 
et al.. l2000l) . (Image credit: Chandra press office, http : //cheindra . harvard . edu/photo/2008/w28/more . html] ! 
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to derive the origin of SNRs (Sect. [T] and [8]l. So one 
now also encounters classifications such as Type la SNR. 
Other SNRs show optical and X-ray evidence for en- 
hanced oxygen abundances; these so-called oxygen-rich 
SNRs (Sect. 19. Il l are likely the remnants of the most mas- 
sive stars. Depending on the context these designations are 
often more helpful than the morphological classifications. 

Note that these new classifications are not mutually ex- 
clusive with the traditional morphological classifications. 
For example oxygen-rich SNRs such as G292.0H-1.8 and 
the Lai-ge Magellanic Cloud SNR B0540-69.3, also harbor 
a pulsar wind nebulae. They can therefore be classified 
both as composite SNRs and as oxygen-rich SNRs. 



4 The hydrodynamic structure and 
evolution of supernova remnants 

4.1 Evolutionary phases 

The evolution of SNRs is usually divided in four phases 
(IWoltien. Il972h : I) the ejecta dominated phase, in which 
the mass of the supernova ejecta, M^j, dominates over 
the swept-up mass, M^v,; II) the Sedov-Taylor phase, for 
which Msw > Mej, but for which radiative losses are not 



energetically important; III) the pressure-driven, or "snow- 
plough" phase, in which radiative cooling has become en- 
ergetically important; the evolution of the shock radius is 
now best described by using momentum conservation; IV) 
the merging phase, in which the shock velocity and tem- 
perature behind the shock become comparable to, respec- 
tively, the turbulent velocity and temperature of the inter- 
stellar medium. 

Although these discrete phases provide a useful frame- 
work to think of SNRs, it should be kept in mind that it is 
an oversimplification, and the phase of an individual SNR 
is not always that easily labeled. Moreover, different parts 
of a SNR may be in different phases. For example, a rem- 
nant like RCW 86 has radiative shocks in the southwest 
(phase III), whereas in the northeast it has very fast, non- 
radiative, shocks (phase I). This is the result of the com- 
plexity of the medium it is involving in, which is probably 
shaped by the stellar wind of the progenitor that created a 
cavity, surrounded by a dense s hell, with which part of the 
shoc k wave is interacting (e.g. Rosado et al.l 19941; Vink 



et al.,ll99j|). 

Phase I is sometimes called the free-expansion phase, 
and phase II the adiabatic phase. Both names are some- 
what misleading. Free expansion suggests that the shock 
velocity is described by Vs = Rs/t, with R^, the radius of 
the outer shock, and t the age of the SNR. However, as de- 
scribed below, even in phase I one has Vs < Rs/t- Since 
energy losses are according to standard models not im- 
portant in both phase I and II, it is also misleading to 
call exclusively phase II the adiabatic phase. Moreover, in 
the very early evolution of SNRs (almost the supernova 
phas e) SNR s may g o through a short radiative los s phase 
(e.g. iTruelove and McKeeL Il999l; ISorokina et all l2004 . 



Also escaping cosmic rays may add to the energy losses 
in phase I and II (see Sect. 15.1b . Despite these shortcom- 
ings, the labels have proven to be of some value, as they 
provide some framework to characterize the evolutionary 
phase of a given SNR. 

In the literature one also often finds design ations for 
SNRs like "young", "mature" and "old" Oones et al. 



19981) . These designations do not have a very precise 



meaning, but a general guideline is that young SNRs are 
less than ^ 1000 — 2000 yr old, and are in phase I or early 
in phase II, mature SNRs are in late phase II, or early 
phase III, whereas the label old SNRs is usually given 
to the very extended structures associated with SNRs in 
phase IV. These old SNRs hardly produce X-ray emis- 
sions, so in this review we will only encounter young and 
mature SNRs. 

4.2 Analytical models 

Several analytical models for the shock evolution of SNRs 
exist. The most widely used is the Sedov-Taylor self- 
similar solution (Sedov, 195Sj; Tavlorl Il950). It assumes 
that the explosion energy E is instantaneously injected 
into a uniform medium with uniform density po (i.e. a 
point explosion), and that no energy losses occur In that 
case the shock radius and velocity Vs will develop as 



dRs ^ 2 / £^ 
dt 5 \ Po 



1/5 _ 2^ 

5 f ■ 



(2) 



(3) 



The dimensionless constant ^ depends on the adiabatic 
index; ^ = 2.026 for a non-relativistic, monatomic gas 
(7= 5/3). An analytical solution exists for the density, 
pressure, and velocity profiles inside the shocked medium, 
which is shown in Fig. |5] The Sedov-Taylor solution 
can be generalized to a gas medium with a power-law 
density profile p(r) o<: r"*: R^ oc t^, Vs = pR^/t, with 
P = 2/(5 — s), the expansion parameter An astrophys- 
ically relevant case is 5 = 2, corresponding to a SNR 
shock moving through the progenitor's stellar wind (see 
below); a situation t hat likely applies to the y oung Galac- 
tic SNR Cas A (e.g. Ivan Veelen et all l2009h . This gives 



j3 =2/3, which should be compared to the experimental 
value found for Cas A b ased on the expansion found in 
X-rays: B = 0.63 ± 0.02 JVink et all Il9 98; Delaney and 
Rudnick ^003l;|Patnaude and FesenL 120091) . 

The Sedov-Taylor solutions do not take into the account 
structure of the supernova ejecta itself. This is a good ap- 
proximation once the swept-up mass exceeds the ejecta 
mass. However, in the early phase, only the outer layers 
of the supernova transfer their energy to the surrounding 
medium. As time progresses, more energy is transferred 
from the freely expanding ejecta to the SNR shell. This 
takes place at a shock separating hot ejecta from f reely ex- 
panding (cold) ejecta, the so-called reverse shock (IMcKee , 
Il974) . 
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Figure 5: The structure of a SNR as given by the self-similar models ofl Sedovl (Il959h (left) and IChevalien (1 1 982h (n = 
1 ,s = 2). The values for the parameters have been normalized to the values immediately behind the forward shock. For 
the Sedov model the radius is expressed in units of the shock radius, for the Chevalier model in units of the contact 
discontinuity Rc, the border between swept-up and ejecta material. The solid lines show the density, the dotted lines the 
velocity, and the short-dashed lines the pressure profiles. For the Sedov model also the temperature is indicated (long- 
dashed line). For this model the temperature goes to infinity toward the center. For the Chevalier model the density goes 
to infinity for r = Rc- Note that in the observer frame the velocity drops at the reverse shock. 



Two analytical models exist to describe the structure and 
evolution of SNRs taking into account t he initial velocity 
structure of the ejecta. The first one, bv lChevaheil(ll982h . 
describes the early evolution of SNRs, in which the freely 
expanding ejecta have approximately a power-law density 
distribution p, 



e] 



This is a reasonable approximation 
for the outer ejecta structure as found in numerical mod- 
els of supernova explosions, with n — 7 a situation that 
describes reasonably well the ejecta structure of Type la 
supernovae, whereas n = 9— 12isa valid approximation 
for the density str ucture of core c ollapse supernovae. 

As shown bv IChevaher ( 1982 ) the SNR evolution can 
be described by a self-similar solution of the form; 

Rs^t^, (4) 
with j3 the expansion parameter given b>0 



(5) 



For s ~0,n = 7 this gives j3 = 0.57, and for n — 9,s = 2 
this gives )3 — 0.86. The models describe a self-similar ve- 
locity, pressure, and density structure. 1 will not reproduce 
them here, but show as an example the density, pressure, 
and velocity structure ofa« = 7,5 = 2 model (Fig. |5]l. It 
important to realize that the Chevalier solutions describe 
the early evolution of SNRs, when the reverse shock has 
not yet reached the inner most ejecta. Once the inner ejecta 
are reached by the reverse shock the expansion param- 



The expansion parameter is often, but not uniformly, indicated with 
the symbol m, but since the characters i,...,n have an inte ger connotation, 
I opt here for j}. Note that lTruelove and McKed ) 19991) use n (but this 
conflic ts with its use in acceleration theory. Sect. 15. 5) . and IChevalieil 
dm uses 



eter will evolve toward the Sedov solution (thus the ex- 
pansion parameter evolves from (5 = (n — 3)/(n — s) to 
P 2/ (5 — . It is clear from the Chevaher solutions that 
from a very early stage on j3 < 1. There is also observa- 
tional evidence for this. For example, the initial expansion 
parameter of SN 1993 J in the g alaxy M81 has recently 
been determined to be j3 =0.85 dMarcaide et al.L l2009h : 
assuming an interaction of the shock with a circumstellar 
wind (s — 2) this implies « ss 8.5. 

An analytical model that takes into account the smooth 
transition from p hase I to phase II was obtained by Tru- 
elove and McKee ( 1999h . Their models employ the fol- 
lowing characteristic length, time and mass scales: 



tch 



(6) 



ej, 



with Mej the ejected mass, E the explosion energy, and 
Po the circumstellar medium density. These numbers can 
be used to construct a set of solutions, which now only 
depend on n and s, and the dimensionless variables R* = 
R/Rch and t* — t jt^-h. The models are continuous, but con- 
sist of two parts; one for the evolution in the ejecta domi- 
nated phase and one for the Sedov-Taylor phase, with t^j 
the dimensionless transition age. For example, the blast 
wave trajectory in the n^7,s — Q model is /?* 1 .06 = t*^^^ 
for t* < tsT and /?* = (1 .42f * - 0.312)2/5 f^j. ^* > ^^^^ ^-^^y^ 
tsT — 0.732. These solutions show that initially the expan- 
sion parameter is i dentical to the one derived bv Cheva- 
her ( 1982!) (Eq.|5]l, while it asymptotically approaches the 
Sedov-Taylor solution, j3 =2/5. 
Fig. |6]illustrates the evolution of SNR shocks as given 
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Figure 6; Shock radii and shock velocities as a function of 
time according to the Truel ove-McKee n = l,s = model 
(ITruelove and McKeeLll999h . The solid line represents the 
forward shock, whereas the dashed line represents the re- 
verse shock. In the bottom panel the reverse shock veloc- 
ity in the frame of the ejecta is shown as a dotted line. The 
dimensionless model was ad justed to fit the properties of 
Kepler's SNR (IVinkLl2008cl) . 



by one of the Truelove & McKee solutions. It shows that 
the reverse shock initially expands outward, despite its 
name, with /3 = (n ~ 3)/(« — s). But as soon, as the shock 
heated shell has more pressure than the ram pressure of the 
freely expanding ejecta, the reverse shock will move to- 
ward the center; i.e. its velocity in the observers frame will 
become negative. For the plasma temperature the shock 
velocity in the frame of the freely expanding ejecta is im- 
portant (Sect. lSTt . This is is given by 



Vr 



R„ 



rev,phys 



rev,obs5 



(7) 



with Rrev/t the velocity of the freely expanding ejecta just 
before entering the reverse shock, and Viev.obs — dRrev/dt, 
the reverse shock velocity in the frame of the observer 
Fig. |6] shows that initially IVrev.physl < V^, in which case 
one expects the shocked ambient medium to be hotter than 
the shocked ejecta. In later phases this situation reverses, 
as |Vrev,phys| > ^s, resulting in a hot core enclosed by a 
cooler shell. Once the reverse shock has reached the center 
of the SNR, it consists of a hot shell surrounding an even 
hotter, but very tenuous interior. 

It is worth mentioning that efficient cosmic -ray accelera- 
tion (Sect. 15.5b may alter the hydrodynamics of SNRs, be- 
cause it alters the equation of state, and because escaping 
cosmic rays may result in energy losses. This can result in 
a reverse shock radius at any given time that is closer to 
the outer shock t han indicated by models l ike that of Tru- 
elove & McKee dPecourchelle et al. [ |2000l) . I will discuss 
in Sect. 1 11. 3] the observational evidence for this. 

In phase 111, the momentum conservation or "snow- 
plough" phase, radiative energy losses have become dy- 



namically important. Instead of energy conservation, the 
radial expansion is governed by momentum conservation, 
MVs ~4-7i/3RlpQdRs/dt = constant. For the constant one 
can take the momentum at the time f^d, for which radia- 
tive losses have become imporant: frad = 4;r7J^^jj/3poVrad- 
Integration yields the following expres sion for the age, t, 
of the SNR as a function of radius (e.g. Toledo-Rov et al.L 
2009h : 



t = trad - 



R 



rad 



4yr, 



■ad 



R 



-1 



(8) 



As can be seen, the expansion parameter in this phase 
is j3 « 1/4. Generally speaking, radiative losses become 
important when the post-shock temperature falls below 
^ 5 X 10^ K, in which cas e oxygen line emis sion becomes 
an important coolant (e.g. lSchure et al.l l2009). This occurs 
w hen the shock velocitv is V« = V,aH = 200 km5~^ ( Wolt- 



ier. ll972l see also Sect. 15. 11 1. One can use the Sedov self- 
similar solution (Eq.|2]and[3]) to calculate when Vs = Vrad, 
thereby estimating Rad and frad- Eq. [3]leads to expression 

for frad in terms of Rmd' trad = ^rad/Vrad = 2.0 X 10"**/?rad S. 

Eq.|2] then provides the following estimate for the age and 
radius at which the SNR becomes radiative: 



frad =1.5 X 10 



(9) 



1.4xl0'2f^V/%«44,600(^) 



£51 \ 1/3 



yr 



Rrad =7.0 X lO'^f ^) cm « 23f ^) '^'pc, (10) 

with «H the pre-shock hydrogen density, £, ^ 2.026, and 
£51 the explosion energy E in units of 10^' erg. 

4.3 Supernova remnants evolving inside 
wind-blown bubbles 



The structure of the circumstellar medium is usually not 
as simple as used by the analytical solutions discussed 
above. The interstellar medium itself is not homogeneous, 
but also stellar winds from the supernova progenitors im- 
pose a specific structure on the ambient media. For ex- 
ample, a massive O-star will spend ^ 90% of its life on 
the main sequence, during which it blows a fast, tenu- 
ous wind with velocity v,^ « 1500 — 3000 kmi~\ with 
M,^ vr~' (e.g. lMokiem et al.Ll2007l) . 



M«5 X 10" 



10- 



This wind carves a low density region into the interstellar 
medium, surrounded by a shock heated shell containing 
the swept-up interstellar medium and wind material. 

In the subsequent red supergiant (RSG) phase the mass 
loss rate increases to as much as M « 10^^ M© yr^' (e.g. 



van Loon et al 



20051), which, together with a much slower 
10 kmi^i, results in a very dense 



wind velocity Vy 
wind, with a density at radius r given by 



M 



(11) 



a relation that follows from mass conservation. Stars 
more massive than 60 Mp, can even have extremely high 
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Figure 7: Left: The density (solid line) and pressure (dashed line) of the circumstellar medium around a red supergiant. 
Middle and right: The evolution of the s hock radiu s and v elocity for a SNR in a wind-blown cavity, surrounded by a dense, 
swept-up shell. All figures taken from iDwarkadasI (l2005h . 



mass loss rates, often accompanied by violent bursts 
of mass loss. These are so-called Luminous Blue Vari- 
ables (LBV s). An example is Eta Carinae ( Davidson and 
Humphrevs. ll997h . 

Most massive stars will explode in the RSG phase, but 
the highest mass stars may enter the Wolf-Rayet star 
phase. In this phase the stars have lost all of their hy- 
drogen envelopes. Their mass loss rates are high, M ss 
10"^ - 10""* Mb yr~', with a high wind velocity 1000 - 
2000 kms^' ( Nugis and LamersL I2OOO). This results in 
a tenuous wind, which will sweep up the RSG wind into 
a dense shell. It is not clear whether LBVs will become 
Wolf-Rayet stars, or whether they are the last stage of the 
stellar evolution. 

Progenitors that were part of a binary system may even 
have had more complicated mass loss histories. In some 
cases a common envelope phase may have existed, which 
could also have lead to a stellar merger In these case the 
circumstellar medium structure may be rather complex. 
See for example the binary evolution model that can ex- 
plain some of the features of the cir cumstellar medium of 
SN 1987A jMorris and Podsiadlow ski, 2007, Secl.W2\. 

So core collapse supernovae are, even in the simplest 
cases, surrounded by regions with different densities and 
temperatures, cre ated and shock heated during the vari- 



ous wind phases (iWeaver et all 119771 : iDwarkadasL 12005 



Fig.|2l). As long as the SNR shock moves through the un- 
shocked stellar wind of the progenitor, the analytical mod- 
els are valid, but if this is not the case, then the evolution- 
ary path may differ quite substantially from the standard 
SNR evolutionary models discussed above. 

These more complicated situations have been studied 
with nume rical simu lations. See for example Tenorio- 
Tagle et al. (Il99ll), |pwarkadas| (|2005b and Fig.|2] In the 
case that the shell swept up by the wind has a mass that is 
higher than the ejecta mass, the shock wave will first move 
through the low density region without much deceleration 
()3 = 0.87, for an n = 9 ejecta density profile, Eq.|5]l. In this 
phase the X-ray emission from the SNR is low, because of 
the low density inside the stellar wind. Once the shock 



reaches the dense shell, transmitted and reflected shocks 
develop. The transmitted shock heats the shell and pushes 
it outward, whereas the reflected shock further heats the 
interior of the SNR. In this phase, the X-ray luminosity in- 
creases rapidly by a factor ten to hundred, due to the high 
density in the shell dPwarkadasl l2005h . 



5 CoUisionless shock heating and 
particle acceleration 

5.1 Shock heating 

Shock waves transform part of bulk kinetic energy into 
thermal energy. Astrophysical shocks are characterized by 
low densities, which makes that the actual heating pro- 
cess cannot be established by particle-particle interaction 
(Sect. 15.2b . For SNRs part of the kinetic energy may not 
only be used for shock heating the plasma, but also for 
cosmic-ray acceleration. 

Irrespective of the microphysical details of particle heat- 
ing and acceleration, the coUisionless shocks still have to 
obey the conservation of mass, momentum, and energ 



flux across the shock (e.g. IZel'dovich and Raized 1196^ 



^ 



McKee and Hollenbachlll980l) . However, one should real 
ize that SNR shocks are not the idealized, infinitely narrow 
temperature transition zones treated in standard physics 
textbooks. Instead, coUisionless heating, and cosmic-ray 
acceleration imply that there is a broad transition zone 
from the unshocked medium (or upstream) to the shocked 
(downstream) region. In the case of efficient cosmic-ray 
acceleration (Sect. I5.5l l, this transition zone, the so-called 
cosmic-ray shock precursor, has a width of approximately 
the diffusion length scale of the cosmic rays with ener- 



gies near the peak in the energy spectrum (e.g. lBlasi et al 
2005l:IVladimirov et al.L l2008). 



Irrespective of whether cosmic -ray acceleration is im- 
portant, mass conservation across the shock remains valid. 
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and for a plane parallel case implies: 



POV'O = P2V'2 ^ P2 = — PO ; 

V'2 



ZPo, 



(12) 



with X the shock compression ratio, and vq and V2 taken 
to be in the frame of the shock. The subscript 2 refers to 
the downstream plasma, and the subscript the (far) up- 
stream medium, thereby symbolically ignoring "region 1", 
the transition zone containing the cosmic -ray shock pre- 
cursor. 

Cosmic-ray acceleration also hardly affects the conser- 
vation of momentum flux: 



Pi + Pivl =Po + Pov'o ^ ■P2 = - 



PoV/, (13) 



with the second part following from substituting Eq. [12] 
and relating vq to the shock velocity, vq = — Vs- For high 
Mach number shock waves the pressure far upstream can 
be neglected (Pq = 0). And if we make the assumption 
that the downstream plasma pressure consists of a thermal 
and a non-thermal (cosmic -ray) part, with the parameter 
w giving the non-thermal (NT) pressure contribution. 



Pi.. 



NT 



Pl.T+Pl. 



(14) 



NT 



we find from Eq. [13] and [14] that the post-shock plasma 
must have a temperature 



kT2^il~w)-{l--]^mpV^. 



(15) 



One of the consequences of efficient cosmic ray accel- 
eration is that the highest energy cosmic rays may escape 
far upstream (i.e. into the unshocked medium), whereas 
downstream (in the shock-heated plasma) the dominant 
flux of cosmic-ray energy is through plasma transport. 
This means that conservation of energy flux is not nec- 
essarily valid. A similar situation holds for so-called ra- 
diative shocks, for which radiative energy losses are im- 
portant, a situation that can occur for shocks with Vs < 
200 kms-'. 

We can take into account of the energ y flux losses. 



Pesc 

^esc 



by introducing flBerezhko and Elhson 



^ .With this equation the conservation of energy flux 
can be written as 

1 , 1 7 

(m2 +^2+ 2P2^2)^2 = (1 - £e.,c):^PQvtvQ, (16) 

with u — Y-^P^ the internal energy, and 7 the adiabatic in- 
dex. Taking for the non-thermal (cosmic ray) contribution 
7 = 4/3 and for the thermal plasma 7 = 5/3 we obtain 



"2+f2=(^W+^)P2 = rP2. 



Using Eq.[T3] we can rewrite Eq.ll6las 

-)+- 



2Y 



(17) 



(18) 




Figure 8: The impact of the partial pressure by non- 
thermal components (cosmic rays) and th e fractio n of en- 
ergy flux escaping as cosmic ravs (Vink, 2008a; Helder 
et al., 2009). The black lines indicate constant values of 
temperature (Eq. [TsT i, given as a correction factor with re- 
spect to post-shock temperature for the standard case (i.e. 
w = Q,£esc = 0, Eq.|20]): j3 = Yr/{^nmpV^). Also indi- 
cated are lines of constant post-shock compression ratio %■ 
The red line depicts the solution cosmic -ra y acceleration 



with the highest possible efficiency (Helder et al.L 12009 



VinketaLlEoioh : the combination of w and Ee: 
of this line are considered unphysical. 



to flie left 



which has the non- trivial root (e.g. iBerezhko and Ellison , 
1999l:IVinkLE008ab : 



^._ r+Vr2-(i-£,,,)(2r-i) ^ (i9) 

1 ^esc 

In most textbooks only the case of a single fluid gas 
without energy losses is considered. This corresponds to 
w = and ee.sc = 0. This gives a compression ratio of = 
4. Substituting this in Eq|15lgives 



kT = ^mpV,^ 



1.2 



Vs 



1000 kms 



-1 



keV, 



(20) 



with ji the average particle mass (jj. w 0.6 for a plasma 
with solar abundances). 

1 used here kT in order to indicate that this refers to 
the average temperature, as different particle species can 
have different temperatures. The reason is that for colli- 
sionless shocks it is not quite clear whether the micro- 
physics that governs the actual heating results in temper- 
ature equilibration among different, or whether different 
species will have different temperatures proportional to 
their mass: kTi °^ nii y2 (Sect. lis. The latter can be ex- 
pected if heating is predominantly governed by scattering 
isotropization of the incoming particles by plasma waves 
in the shock region. 

Finally, if cosmic-rays contribute substantially to the 
post-shock pressure (w ^ 0) and, if cosmic-ray escape 
plays an important role, then the compression ratios can 
be much higher than ;if = 4. In this so-called non-linear 
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cosmic-ray acceleration case, there is no longer a unique 
value for kT and but a range of values is possible. The 
possible values of the compression ratio and post-shock 
plasma temperature as a function of post-shock cosmic - 
ray pressure fraction, and escape flux are summarized in 
Fig. [8] In this figure the mean plasma temperature kT is ex- 
pressed in terms of the temperature for the w = QlEesc = 
case, using j3 = W /{'ipLm^V^ /IG). 

In all cases non-linear cosmic ray acceleration leads to 
lower plasma temperatures, in the most extre me cases per 



haps e ven quenching thermal X-ray emission (IDrurv et al 



12009'). On a macroscopic level this can be understood by 
considering that part of the incoming energy flux is used 
for particle acceleration. On a more detailed level it can 
be understood by considering that directly upstream of 
the shock the cosmic rays provide a non-negligible pres- 
sure, which pre-compresses and pre-heats the plasma. The 
presence of the cosmic -ray precursor results therefore in 
a lower Mach number gas-shock For a more thorough 

treatment of the shock equations, including a cosmic ray 

I — 1 I — - — I 

precursor, see Vink et al. (2010) for an analytical treat - 

ment, or yiadimii'ov et al. (.2008) ; ICaprioh et"aD (l2008h 

for a kinetic approach. In the latter two references also 

the effects of non-adiabatic precursor heating is taken into 

account. 

5.2 CoUisionless shocks 



Shock waves in the interstellar medium (ISM) are char- 
acterized by very low densities, typically « ^ 1 cm^^. In 
high density shocks, such as in the earth's atmosphere, the 
shock heating process is accomplished by particle-particle 
interactions. However, the mean free paths for particle- 
particle interactions in shocks in the ISM are long com- 
pared to the typical sizes of SNRs. This can be illustrated 
as follows. For charged particles the energy exchange is 
governed by Coulomb interactions. A particle is deflected 
by 90° in the center of mass frame, if the following condi- 
tion for the impact parameter b is met: 



1 



2 m\ + 1712 



m\m2 ^2 



(21) 



with mi, m2 the masses of the two particles, Zye, Z2e their 
charges, and v their relative velocity. The cross section 
for such a deflection by a single scattering is (Ocouiomb = 



■2-72„4 



O'couiomb ~ 4;r- 



ZfZ 



/mi +m2y 
V mim2 / 



(22) 



The corresponding collision time scale is Tcouiomb = 



l/{nav] 



,-3 



°^ E . Inserting for the mi and m2 the 



proton mass and a typical velocity of v = 1000 kms 



one finds for proton-proton collisions Tpp w 10 Wp s 
(about 32,000 Wp ' yr), and a mean free path Ap w 
lO^^Hp ' cm (32 pc). These are much larger than the ages 



^This is the actual shock that heats the plasma, but it is often called 
sub-shock in order to set it apart from the overall shock structure that 
includes both the precursor and the gas-/sub-shock. 



and radii of young SNRs. As these young SNRs clearly 
do have shocks and hot. X-ray emitting, plasma, as beauti- 
fully shown by Fig.lT] the formation of the shock and the 
plasma heating cannot be the result of Coulomb interac- 
tions. Instead "collective interactions", occurring through 
fluctuating electric and magnetic fields, must be respon- 
sible for t he plasma heating. Simulations (e.g. Bennett 
and EUison, [1995^ show that the heating in such shocks 
takes place over a distance of typically 10-100 times 
c/cOpe, with cOpe = (47re^ne/we)'''^ the plasma frequency. 
This corresponds to a shock thickness of roughly Ax = 

7 — 1/2 

10 He cm, thirteen orders of magnitude smaller than 
the range over which Coulomb collisions operate. For this 
reason SNR shocks and shocks in other low density media, 
such as the interplanetary medium, are called collisionless 
shocks. 

The heating mechanism is s omewhat analogous to the 
process of "violent relaxation" ( Lvnden-Belj 1967 ) in the 
formation of bound gravitational systems, such as galaxies 
and clusters of galaxies; i.e. large scale gravitational po- 
tential fluctuations rather than two-body interactions are 
important for the relaxation of the system. 

5.3 Temperature equilibration 

It is not a priori clear, whether collisionless shocks will 
result in temperature equilibration of all t ypes of parti- 



cles im mediately behind the shock front (see lBvkov et al 



2008al for a review). The two opposite cases are a) com- 
plete equilibration, in which all different particle popula- 
tions (electrons, protons, other ions) have the same tem- 
perature, which for a strong shock without cosmic-ray ac- 
celeration would mean T °^ jini^V^ (Eq. [TsT i, or b) full 
non-equilibration, in which case each particle species ther- 
malizes the bulk kinetic energy independently and one ex- 
pects 7] niiV^ (corresponding to replacing jim^ with m, 
in Eq. [T5]andl20]iPI There are many plasma instabilities 
that act to equilibrate the energies of different plasma con- 
stituents. Both hybrid simulations (Cargill and Papadopou- 
los, 1988) and part icle in cell (PIC) simulations (e.g. Shi- 
mada and Hoshino, 2000h show that high Mach number 
shocks may substantially heat the electrons, but not com- 
pletely up to the same temperature as the protons. So 
one may expect m ^ /nip <^ T^/Tp < 1. For example. Shi- 
mada and Hoshino (|2000|) report T^/Tp « 17%. Optical 
spectroscopy of SNR shocks show mixed results with 
indications that full equilibration is important fo Vs < 
50 km.;"', an d non-equilibration for higher velocities 



(iGhavamian et al ., 2007; Held ereta l., 201 1). In Sect. HTTl] 
I discuss the implications for X-ray spectroscopy. 

Note that these simulations have several limitations. For 
example, PIC simulations often have mp/mg ~ 0.01 — 1, 

' The latter case con'esponds to a velocity distribution that is indepen- 
dent of the particle mass. This is again similar to the outcome of violent 
relaxation for bound gravitational systems. For example, in galaxy clus- 
ters the velocity distribution is similar for massive and dwarf galaxies 
(with the exception of the giant ellipticals that reside in a large fraction 
of cluster centers). Also in this case the kinetic energy ("temperature") 
distribution is proportional to the mass of the galaxy. 
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for computational reasons. In contrast, hybrid simulations 
treat the electrons as a massless fluid. Moreover, the 
amount of equilibration of electron/proton temperatures 
may depend on the initial conditions, such as Mach num- 
ber, the ratio of magnetic over kinetic pressure and the 
magnetic-field orientation. This means that some variation 
in shock equilibration can be expected in SNRs. 

An important ingredient that affects electron heating, 
and one that is difficult to model, is cosmic -ray acceler- 
ation (see Sect. 15.5b . because the process of acceleration 
to high energies takes place on a longer time scale than 
the formation of strong shocks, and longer than the typi- 
cal time scales that can be modeled with PIC simulations. 

Neverthe less, hybrid and PIC simulations (e.g. Bennett 
and Ellison, 1995; Dieckmann and Breti 120091) do show 
evidence for the initial stages of particle acceleration: of- 
ten the particle distribution has non-thermal components, 
i.e. the particles have roughly a Maxwellian energy dis- 
tribution at low energies, but a non-thermal, quasi-power- 
law distribution at high energies. This non-thermal tail, 
may be the first stage of particle acceleration, but part of 
this non-thermal distribution, may disappear again further 
behind the shock, as the particles are slowed down through 
Coulomb interactions. This is especially relevant for the 
low energy part of the n on-thermal electron distribution 
(Sect. l6372llVinkil2008b[) . and amounts to a source of late 
time electron heating. 

At the shock front plasma wave coupling between elec- 
trons and accelerated ions may als o result in heating of the 
electrons (Rakowski et al. , Hogs'). Finally, very efficient 
particle acceleration by the Fermi process may alter the 
shock structure altogether, as it results in the formation of 
a cosmic -ray shock precursor that may pre-heat both elec- 
trons and protons before they enter the main shock (see 
Sect. |53] and Sect.lSTTT). In Sect. lU.ll l review the obser- 
vational measurements concerning electron/ion equilibra- 
tion. 

If the shock heating does not immediately lead to equal 
electron and ion temperatures. Coulomb interactions will 
slowly equilibrate them. It turns out that Eq. |22] is not a 
good approximation for the relevant cross sections, as it ig- 
nores the contributions for long range interactions. These 
result in less energy exchange per encounter, but they oc- 
cur much more frequently. The relevant expressions for en- 
ergy exchange between charged particles of various mass 
ratios can be found in Huba (2002), which gives for the 
tempe rature eauihbration rate (see also ZeFdovich and 
Raizer. ll966l:lltohLll984 : 



^=1.8xlO-"^x 
dt 



(23) 



Lj {kPi -kTi)eWs-\ 

[niikTi+mikTifl^ ^ ' 

with label 1 for the particle species of interest, ; for all the 
other species, and kT measured in units of eV. InA is the 
Coulomb logarithm, which is roughly similar for electron 



A 2 




Figure 9: An illustration of the effect of temperature non- 
equilibration at the shock front for a shock velocity of 
Vs — 3000 km ' . Shown is the temperature of electrons 
(dotted), protons (solid), helium (dashed-dotted) and oxy- 
gen ions (dashed) as a function of rigt, assuming that heat- 
ing at the shock front is proportional to the mass of the 
particle (Sect. 15.1b . The oxygen-proton equilibration is 
faster than the helium-proton equilibration, as the cross 
sections scale linearly with particle mass, but quadrati- 
ca lly with cha rge (Eq. |22] |. (Figure previously published 
in lVinkl(l2006l) .) 



proton interactions and for electron-electron interactions: 



InA = 30.9 -In 



VlkeV/ 



(24) 



For electron-proton equilibration Eq.|23]implies an equili- 
bration time scale: 



:3.1xlO'V-'f^)'/VM)"'s, (25) 
P VlkeV/ V30.9/ ' 



with T the mean temperature. Similarly, the electron self- 
equilibration time, which determines how fast the elec- 
trons can establish a Maxwellian distribution, is 

Tee«4.9xlOV'f^)'^V — r's. (26) 

" VlkeV/ V30.9/ 

Eq. I25] shows that it can easily take 10,000 years for 
electrons and protons to equilibrate, much longer than the 
ages of many of the young SNRs that will be discussed 
in the second part of this review. Note that the equilibra- 
tion time is roughly inversely proportional to np(?s We)- A 
relevant parameter is, therefore, not so much the age of 
the SNR, but the average «ef , a parameter that will resur- 
face in Sect. 16.1.21 where non-equilibrium ionization is 
discussed. In fact. X-ray spectroscopy provides a direct 
way of measuring «ef . thereby providing the means to esti- 
mate whether non-equilibration of temperatures could be 
important. 

The possible absence of full equilibration has impor- 
tant conseq uences for the interpretation of X-ray obser- 
vations (see I1978L for an early discussion), because 
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in general the X-ray spectra reveal the electron temper- 
ature (Sect. 16.11 ), and not the mean plasma temperature. 
The mean plasma temperature is energetically the relevant 
parameter. It is for example directly related to the shock 
velocity (Sect. 15. 11 1. Using the electron temperature, as 
an estimate for the mean temperature may lead to a seri- 
ous underestimation of the shock velocity. But, as shown 
in Sect. 15.11 efficient cosmic -ray acceleration will lead to 
lower mean temperatures. So even having a good estimate 
of both the electron and ion temperatures may not be suf- 
ficient to estimate the shock velocity. 

Fig.|9]illustrates the effects of non-equilibration of tem- 
peratures at the shock front. The equilibration time de- 
pends on the density, the mass ratio and the square of 
the charge of the particles (Eq. |23] ). The ion-ion equili- 
bration, therefore, proceeds faster than electron-ion equi- 
libration. Interestingly, iron-proton equilibration proceeds 
faster than proton-helium equilibration, due to the charge 
dependence. This becomes important for «ef > 10^ cm^^s, 
when iron has become stripped of all m-shell electrons. 
Broadly speaking one can recognize three phases: 

1 . Full non- equilibration; all species have different tem- 
peratures {net< 5 X 10^" cm^^s), 

2. Partial non- equilibration; all ions are equilibrated, 
but without electron-ion equilibration 

(5x 10'°<«ef< lO'^cm-^s), 

3. Full equilibration in^t> 10^^ cm^^s). 

It is worth noting that non-equilibration of ele ctron- 
ion te mperatures may also be caused by ionization fttohL 
19841) : Due to the fact that the cross sections for ionization 
peak near the ionization energy, the electrons originating 
from post-shock ionization will start relatively cool. The 
electron-electron equilibration is relatively fast, so this re- 
sults in a relatively cool electron population, even if the 
free electrons that entered the shock were equilibrated 
with the ions. This may be an important effect in metal- 
rich plasmas, which have been shock heated by the reverse 
shocks in young SNRs (Sect.|4|i, since in that case most of 
the electrons will originate from ionizations, rather than 
from entering the shock as free electrons. 



netic fields, or along magnetic fields, thermal conduction 
is usually mediated through electrons, because they have 
higher thermal speeds and are more strongly coupled to 
each other. Across magnetic-field lines the particles with 
the largest gyroradius, i.e. the ions, will dominate thermal 
conduction. 

The process of thermal conduction is described by 
Fourier's law: 

Fheat = -KVT. (27) 

In the absence of magnetic fields, or along magnetic field 
lines, the coeffcient K depends on the typical mean free 
path for energy exchange, X ~ TeeVth, with Vth — y/kTeJm^ 
and Tee the self-equilibration time (Eq |26] ). Numerically, 
X w Q.2{kT 1 1 keV)^ pc. The amount of heat being trans- 
ported relates to the energy density u = 3/2ne^7e and 
the transport velocity Vth, so that approximately K/k ^ 
2>/2nekTXvi^ = 2>/2nJ<:TXf:^v^^. This is correct up to a 
smal l factor that is needed to ensure zero electric current 
(e.g. Spitzer, '1965*). Taking this into account one obtains 
(.Huba.' 2002-1^ 



K«/k = 23 



npkTf. 



:2.0X 10 



,27 



(t 



kT, \5/2 



V 



keV/ 



(28) 

This shows that thermal conduction is strongly tempera- 
ture dependent. 

Eq. |27]is only valid when the temperature scale height 
is larger than the mean free path for energy exchange. If 
this is not the case then the heat flux is described by the 
so-called saturated heat flux (iCowie and McKeelll977h is 



/2;t7;\i/2 
f;at = 0.4 ^ n^kT,. 



(29) 



In this case heat is transported by electrons streaming 
away from hot regions, but heat exchange between elec- 
trons is not taken into account, i.e. the electron distribu- 
tion will not be Maxwellian. The factor 0.4 takes into ac- 
count that the heat flux should be electrically neutral, re- 
ducing the heat flow. Unlike the classical heat flux, the 
saturated heat flux is density dependent. The ratio of the 
classical over saturated heat flo w, involving a temperature 
scale height parameter 7?s is (c.f lCowie and McKeelll977D : 



5.4 Thermal conduction 

Thermal conduction is a topic that every now and then 
enters the discussion of SNRs, as significant thermal con- 
duction may alter the tempera ture and density struc tures 
of SNRs (see for example Cui and CoxL 1992 . and 
Sect. fTOJT l. 

The role of thermal conduction in SNR structure and 
evolution, like that of temperature non-equilibration, has 
never been satisfactory resolved. The reason is that ther- 
mal conduction in a plasma is likely to be anisotropic due 
to the inhibiting effects of magnetic fields. Note that the 
processes of thermal conduction and temperature equili- 
bration are related, because thermal conduction is a com- 
bination of transport and exchange of heat. Without mag- 



_ 7.2Tee y/kTe/r, 



= 1.5n 



-1 



1 pc 



-1/ kTe \2/lnA 

.IkevJ \3Q3 
(30) 

The saturated heat flux is, therefore, dominant for small 
scale sizes (/?s << 1 pc). 
The thermal conduction across field line s, which is pre- 



dominantly mediated by ions (protons), is (ISpitzen. 11965 
HubaLl2002h : 



2.8n„kTr, 7 
Kj_ k^ ^—^ = 1.4 X 10^ 



B 



lOliG 



cm-'s"', (31) 



with Ocp = eB/{m-pc) the proton gyrofrequency, and Tpp 
the proton self-equilibration time scale. This shows that 



Correcting for a slightly different definition of Tec • 
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there is a very large difference in conduction parallel and 
across field lines. Only for B ^ 10^^^ G does the perpen- 
dicular and parallel conduction become comparable (for 
n—\ cm^^). 

To get an order of magnitude estimate for the impor- 
tance of thermal conduction one should compare the ad- 
vective enthalpy flux F-^ = v^nkT (c.f. Eq [T6b with the 
thermal conduction flux, approximating the temperature 
gradient with kT /R^ and assuming a typical flow speed v: 



Fheat 2 k- V 1 pc 



1000 kms- 



IkeV 



-5/2 



(32) 

where for K the value along the magnetic field has been 
used. The fact that the numerical value is around one sug- 
gests that in SNRs the issue of thermal conduction is criti- 
cal. This in turn makes the regularity of the magnetic field 
a crucial parameter. Note that the advected over conduc- 
tive heat flow ratio (Eq[32]) is proportional to {kT)^^l^, if 
in addition one assumes v ^ VfcT, this ratio will be pro- 
portional to (kT)^. This means that thermal conduction is 
more important for hotter plasmas. However, there are in- 
dications that the magnetic fields in young SNRs, which 
have hotter plasmas, are more turbulent (Sect. 111. 3b . re- 
ducing the thermal conductivity. On the other hand little 
is known about the magnetic fields in the ejecta compo- 
nents of SNRs. This may, as we see be of interest for the 
interior components of some mature SNRs (Sect. [T03T l. 

5.5 Cosmic-ray acceleration by supernova 
remnant shocks 



SNRs are considered to be the prime candidate s ource s 
for cosmic rays (e.g. lOinzburg and SvrovatskiiL Il969h . 
at least for energies up to 3 x 10^^ eV, at which energy 
the cosmic-ray spectrum steepens from a power law with 
slope -2.7 to a slope of -3.1. This spectral feature is often 
referred to as the "knee". There is evidence t hat ar ound the 
"knee" the composition changes ( Horandel 2008 . for a re- 
view), suggesting a rigiditjQ dependent maximum cosmic- 
ray energy at the source, with more massive particles hav- 
ing their spectral breaks at higher energies. It is thought 
that only cosmic rays above ^ 10^^ eV may have an extra- 
galactic origin. 

The main reason to consider SNRs as the dominant 
sources of Galactic cosmic rays is that the inferr ed Galac- 

2 eV cm^^, Webber , 



tic co smic -ray energy density (^ 1 

19981) . combined with the inferred time an average cos- 
mic ray spends in the Galaxy, requkes a Galactic cosmic- 
ray production with a total power o f about lO'*' ergs^' 
(e.g. lGinzburp and Svrovatskiilll967[) . Supernovae are the 
only sources known to provide such a power: with their 
average explosion energy of ^ 10^' erg and an inferred 
Galactic supernova rate of 2-3 per century, they provide a 



'The rigidity of a particle is R = p/q, with p the momentum and 
q = Ze the charge. Massive particles, having a higher charge, will have a 
smaller gyroradius, and will therefore less easily diffuse away from the 
shock. 



total power of ~ lO^^ergs^'. This requires that on aver- 
age about 10% of the initial explosion energy is used for 
accelerating cosmic rays. Note, however, that the cosmic- 
ray energy density is very much dependent on the lowest 
energy part of the cosmic-ray spectrum, of which is little 
known, as low energy cosmic rays are shielded from the 
solar system by the solar wind. 

Although supernovae are the most likely energy source 
for cosmic-ray acceleration, there is a debate on how and 
when their kinetic energy is used to accelerate cosmic 
rays. The general view is that cosmic rays are mostly ac- 
celerated in the (early) SNR stage. An alternative view 
is that the collective effects of multiple SNR and stellar 
wind shocks inside star-forming regions are responsible 
for most of the Galactic cosmic-ray acceleration, and that 
this provides a mechanism to accelerate cosmic rays up 
to energies of lO'^ eV (e.g. Bvkov and FleishmanL 1992 : 



Parizot et all 12004 . 



A long standing argument to consider the SNRs as the 
main source of Galactic cosmic rays is that non-thermal 
radio emission from SNRs provides clear evidence for 
particle acceleration, although strictly speaking this only 
provides evidence for electron acceleration, whereas the 
cosmic-ray spectrum observed on Earth consists for 99% 
of protons and other ions. Over the last 10-15 years the 
evidence for cosmic -ray acceleration in SNRs has been 
reinforced by the detection of X-ray synchrotron radia- 
tion (discussed in Sect. II 1.2b . and TeV y-rays from sev- 
eral SNR s. See Sect. 1 1.4land the reviews bv Hinton and 
Hofmann ('2009') and|Reynolds| (|2008'). 

The acceleration of particles in SNRs is likely the re- 
sult of the so-called first order Fermi proc ess, also often 
referred to as diffusive shock acceleration ("Axford et all 
Il97 7; Bell 1978a: Blandford and Ostri ker. 1978: Krym- 
skii. ll977ir In this process some charged particles repeat- 
edly scatter back and forth across the shock front. On 
either side of the shock the particles scatter elastically, 
due to the presence of a turbulent magnetic field. After 
each crossing the particles will be pushed along with the 
plasma, resulting in a net energy gain given by 



4vo_ 
3 



-V2 



3V X 



(33) 



with vn, and V2 the plasma velocities upstream and down- 
strearrlj, seen in the frame of the shock, and vi —Vs — X'^i 
('Eq.ll2b. The particles downstream will be advected away 
from the shock by the plasma with a rate «crV2 ~ ncuVs/ %, 
but the rate at which particles scatter back into the up- 
stream medium is jncRC for an isotropic particle distri- 
bution]^ with c the velocity of the particle. The ratio of 
these two rates gives the chance that the particles escape 



'^To be more precise: if there is a cosmic-ray precursor one should 
write v{x) — V2 with v{x) the velocity at coordinate x within the precursor, 
and I'o < v{x) < V|, with subscripts as defined in Sect. 15. ll So pailicles 
of a given energy may not sample the total change in shock velocity over 
the shock structure, but only a fraction of it. 

' 1 /4c is the projected velocity to the shock normal for particles dif- 
fusing upstream. 
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downstream: 



«CRV2 

jncRC 



, ^2 



(34) 



The chance that a particle will be scattered up and down at 
least ^-times is P{n > A:) = (1 — Pesc)*^- A particle that has 
crossed the shock k time will on average have an energy 
E fa £0(1 +AE/E)''. Eliminating k gives 

The integrated probability function is therefore a power 
law with slope —^+1: 



dlnP In(l-fesc) 



dlnE ln(l+A£'/£') 

4ys/(cx) _ 3 



(36) 



4/3U-i)/zVs/c ;t-r 

The differential distribution function is then 

n(E)dE oc E-idE, q = ^-t^. 

Z-1 



(37) 



In order to give an idea about how many shock crossings 
are needed to reach 10^^ eV, assume that a proton is ac- 
celerated by a shock with Vs = 5000 kmi^\ starting with 
an initial energy of 100 keV, and assuming a shock com- 
pression ratio of X According to Eq.|33]the gain per 
shock crossing is 1.7%, therefore the number of scatter- 
ings is log(10'^ eV/lO^ eV)/log(1.017) w 1400. 

One of the outstanding uncertainties is how charged par- 
ticles are injected into the Fermi acceleration process, i.e. 
how do they reach sufficient velocity to make it back to the 
unshocked medium. It is often assumed that these are the 
particles in the high velocity tail of the Maxwellian distri- 
bution, but computer simulations show that immediately 
behind the shock a non-thermal distribution of particles is 
present that may act as seed pa rticles for further accelera - 
tion by the Fermi process (e.g. [ Bvkov and Uvarov . 1 19991 : 



Riquelme and Spitkovskv . 201 ll for the case of electrons'). 

The typical time accelerate a particle from an initial mo- 
mentu m p, to a final momentum Pf is (e.g. Malkov and 
Drurv. l200lh 



^acc — 



V2-V0 



J Pi ^ 



PffDo 



V2 / 



(38) 



with p the momentum (E ~ pc for ultra-relativistic parti- 
cles), D the diffusion coefficient, and v the plasma veloc- 
ity, with subscripts and 2 referring to the upstream and 
downstream regions. The diffusion coefficients are usually 
parameterized in terms of the so-called Bohm-diffusion 
coefficient, for which the particle mean free path corre- 
sponds to the gyroradius: 

Ec 

D^r,-, (39) 
with 77 expressing the deviation from Bohm diffusion. 



This suggest s a tvpical acceleration time scale of ( Pari- 
zot et al., 200^ 
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yr, 



with B_4 the downstream magnetic field in units of 
100 jUG and X4 the overall compression ratio in units of 
4. The factor 1.83 comes from taking into account the 
difference in diffusion coefficient between the upstream 
(ahead) and downstream (post-shock) region, under the as- 
sumption that upstream the magnetic field is highly turbu- 
lent and isotropic, whereas the downstream magnetic field 
is determined solely by the compression of the magnetic- 
field component perpendicular to the shock normal. Note 
that a higher compression ratio corresponds to a longer ac- 
celeration time, but that a high compression ratio itself is 
a result of efficient acceleration (Sect. lSTI ). 

The typical l eng th scale over which diffusion dominates 
over advectiorl^ is given by /^iff ~ 2D/v, which for the 
downstream (post-shock) region becomes 



/diff w 2 



DiX 



-SJxlQ'^BzlnxA 



5000 km s" 



-1 



(41) 
- I en 



lOOTeV/ 



The diffusion length scale gives approximately the region 
from which particles situated in the downstream (shock- 
heated) region are still able to cross the shock front. 

Eq|32] shows that a compression ratio of 4 corresponds 
to a particle index of ^ = 2. Note that this result is indepen- 
dent of the diffusion properties (slow or fast) as long as the 
diffusion is isotropic. This slope is close to what is needed 
to explain the cosmic-ray spectrum observed on Earth 
{q = 2.7), provided one takes into account that cosmic 
rays with higher energies escape faster out of the Galaxy 
than low energy cosmic rays. Note that for q = 2 each 
decade in particle energy contributes the same amount of 
energy to the overall cosmic -ray energy budget. For q <2 
the highest energy particles contain most energy, whereas 
for q > 2 the lowest energy particles contain most of the 
energy. For the lower energy limit one can take the parti- 
cle's rest mass energy, as around that energy a break in the 
cosmic -ray spectrum is expected (Bell, 1978b) . 

Diffusive shock acceleration becomes more compli- 
cated when particle acceleration becomes very efficient 
and the accelerated particles contribute a significant frac- 
tion to the overall pressure. This situation is usually re- 
ferred to as non-linear cosmic-ray acceleration. As ex- 
plained in Sect. 15.11 this results in a lower downstream 
plasma temperature and an overall compression ratio 
Xtot > 4. The main shock (subshock) will have a com- 
pression ratio Xgas < 4. The compression ratio sampled 



To see this consider tiiat the length scale associated with diffusion 
over a time t is /jiff = ^/2D^, whereas for advection /^jv = V2t. These 



two are equal for / = ID2IV2 = IDiXl'^s- 
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by a diffusing particle may therefore depend on its en- 
ergy/diffusion length scale. The highest energy particles 
experience a compression ratio of Xtoi and the low energy 
particles x$ii&- As a result the spectrum at low energies is 
steep, with q> 2 and flattens at higher energies. For very 
efficient acceleration Eq. [32] is not valid, but instead the 
spectrum approaches the limit ^=1.5 (MalkovLll997h . 

There is observational evidence that the particle spec- 
trum is indeed steeper at low energies than at high energies 
(Sect. III.2I 1. However, there is no proof yet that the spec- 
trum at high energies becomes as flat as ^ = 1.5. Some 
physical processes may prevent that, such as a possible 
tendency for Alfven waves to move away from the shock, 
reducing the velocity gradient tha t the particles experience 
( Zirakashvili and Ptuskin , 2008V Alternatively, an addi- 
tional pressure term in the precursor, for example caused 
by amplified magnetic fields and cosmic ray induced heat- 
ing, may prevent very high compression ratios in the pre- 
cursor, the reby limiting the ove rall compression ratio to 



7 4 - 7 (iMorlino et all 120071) . 



6 X-ray radiation from supernova 
remnants 

6.1 Thermal X-ray emission 

The hot. X-ray emitting, plasmas created by SNR shocks 
have two characteristic properties: they are to a very good 
approximation optically thin, and the ionization distribu- 
tion of atoms is often out of equilibrium. 

Optically thin. X-ra y emitting, plasmas are a l so fou nd in 
a clusters of galaxies (iBohringer and Werner , l2009h . col- 
liding stellar winds, and the coronae of cool stars like our 



sun (hence the name coronal plasmas, see iMewel 11999 



for details on emission processes). In some cases, when 
the lines are not too much broadened by thermal or tur- 
bulent motion, line emission may be affected by optical 
depth effects, caused by resonant line scattering ( Kaastra 
and Mewe, 1995). The overall spectrum will not be af- 
fected by this process, as resonant line scattering does not 
destroy photons, but merely changes the direction of pho- 
tons. In coronal plasma's excitations and ionizations are 
caused predominantly by the electrons colliding with the 
ions. Most of the spectral characteristics of the thermal 
emission (continuum shape, emission line ratios) are de- 
termined by the electron temperature, which, as explained 
in Sect. I5.2l is not necessarily the same as the proton/ion 
temperature. The most important spectroscopic effect of 
the ion temperature is thermal line broadening, which is 
difficult to measure with the current generation of X-ray 
spectroscopic instruments (but see Sect. 111. fl . 

Because SNR plasmas are optically thin for X-rays, X- 
ray spectroscopy is a powerful tool for measuring abun- 
dances in SNRs, as no detailed radiative transfer models 
are needed. For old SNRs this can be used to reliably 
measure the abund ances of the interstellar medium (e.g. 
Hughes et al.Lll998h . whereas for young SNRs X-ray spec- 



troscopy is an important tool for measuring abundances, 
and connect them to supernova ejecta yields for the vari- 
ous types of supernova. Here I discuss the basic radiation 
mechanisms. In depth r eview s of thermal X-ray em ission 
are provided bv lMewd (fl999h : lKaastra et all (l2008h . 



6.1.1 Thermal continuum emission 

Thermal X-ray spectra consist of continuum emission 
caused by bremsstrahlung (free-free emission), recombi- 
nation continuum (free-bound emission), and two-photon 
emission, the latter caused by the radiative electron transi- 
tion from a meta-stable quantum level. For a Maxwellian 
energy distribution of the electrons, the emissivity is given 
by 



2^Ke^ / 27r \ 1/2 , ^ _i/2 / hv \ 

X Me ^ njZf erg s^ 'cm^^Hz^ ' , 

with ^ff ~ 1, th e gaunt-factor, which has a fre quency de- 
pendence (e.g. iRvbicki and Lightma nL ll979h . The sub- 
script / denotes the various ion species, each with charge 
eZj. The emissivity at a given temperature is therefore de- 
termined by the factor n^Y^i^iZf. For solar or sub-solar 
abundances, bremsstrahlung is dominated by electrons 
colliding with protons and helium ions. In spectral codes 
one therefore usually parametrizes the normalization fac- 
tor in Eq. |43](neLi«i-Zf) with «e«H or n^. Since one ob- 
serves the emission from plasma of a given volume, spec- 
tral fitting codes usually use as a normalization factor for 
the plasma / n^niidV or / n^dV . This quantity is called the 
emission measure or EM. Of course one should also take 
into account the distance, in order to convert to flux units. 
A popular X-ray spectral fitting code like xspec uses the 
normalization factor / n^n^dV / {And^). 

Using this type of normalization is valid as long as the 
helium to proton ratio is close to the cosmic ratio, and the 
continuum contribution for elements beyond Z — 2 can 
be neglected. However, for shocked supernova ejecta elec- 
trons colliding with heavy ions can also be an important or 
even dominant source of continuum radiation. If one does 
not take this into account one can derive erroneous density 
and mass estirnates f rom the bremsstrahlung emissivities 
(e.g. lVinketal.lll996 ). Note that the higher emissivity of a 
metal-rich plasma is not only caused by the factor Z?, but 
also by the effect that ionized metals produce more free 
electrons. For example, each hydrogen atom contributes 
only 1 electron, whereas fully ionized oxygen produces 8 
electrons. 

Free-bound and two-photon emission are not always 
mentioned as sources of continuum radiation, but they 
can be a dominant source of continuum emission ( Kaas- 
tra et al., 2008h . in particular for the metal-rich plasmas in 
young SNRs. Free-bound emission arises as an electron 
is captured into one of the atomic shells. The energy of 
the emitted photon is then hv„ = Eg + %„, with E^. the en- 
ergy of the free electron, and %„ the ionization potential 
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Figure 10: The emissivity of a pure silicon plasma 
out of ionization equilibrium (kTs; — 1 keV, n^t— 5 x 
10'° cm- 
emission 

dashed line) and bremsstrahlung (free-free emission, red 
dotted line). The total emissivity is also shown, including 
Si-L and Si-K shell line emission (based on calculatio ns 
made with the spectral code ^PisX iKaastra et al.L 120031) . 



^s). Shown are the contributions of two-photon 
(red solid line), free-bound continuum (red 



for an electron in level n. An electron is more likely to be 
captured in a shell with high principal quantum number 
n, due to the high statistical weights of these shells. As 
the energy difference for levels with high n are small, free- 
bound emission shows a sharp edge near the series limits. 
For a given energy level the emissivity is given by 



efb 



/2V/2 g, 
gi+\ 

CCT(/.v)f^)f^) 

^ 'KrJWc^kTJ 



X exp 



■ Xn ■ 

[hV-Xn] 



3/2 



(43) 



ergs 'cm ^Hz ', 



with n,+i the density of a ion with charge z+ 1, gz.+iTgz 
the statistical weights of the ion before and after recombi- 
nation, and a{hv) the photo-ionization cross section of 
the ion in its final state. Since a{hv) °^ V^^, the spec- 
tral shape for hv ^ Xn resembles that of thermal brems- 
strahlung. However, if kT <C hv, a situation that will oc- 
cur in photo-ionized or overionized plasmas, free-bound 
emission results in narrow emission peaks near the series 
limits of lines. These line-resembling features are called 
radiative-recombination continua, which is usually short- 
ened to RRCs. RRCs have recently been identified in the 
X-ray spectra of a few mature SNRs, suggesting the pres- 
ence of overionized plasmas (Sect. [T03T l. 

Two-photon emission results from electrons in meta- 
stable states, such as the 2s state of a hydrogen-like atom. 
Since decay to the Is level is forbidden (because As = 0), 
it can either be coUisionally de-excited (unlikely in the 
rarified plasmas of SNRs), or it can de-excite by emitting 
two photons, with the associated energy distributed over 
two photons. 



The different types of thermal continuum processes are 
illustrated in Fig. [TO]for the case of a silicon-rich plasma 
out of equilibrium with kT^ = 1 keV. 

6.1.2 Non-equilibrium ionization 

The discussion of thermal emission so far pertains to all 
coronal plasmas, whether they are found in clusters of 
galaxies, SNRs, or cool stars, with the additional compli- 
cation that young SNRs can have very metal-rich plasmas. 
But there is another important difference between the opti- 
cal emission from SNRs and other hot astrophysical plas- 
mas: SNR plasmas are often out of ionization equilibrium. 
This is usually indicated with the term non-equilibrium 
ionization, or NEI. The plasmas of cool stars and clusters 
of galaxies are referred to as collisional ionization equilib- 
rium, or CIE. 

The reason that SNR plasmas are in NEI is simply that, 
for the low densities involved, not enough time has passed 
since the plasma was shocked, and per ion only a fe w ion- 
izing collisions have occurred for any given atom jltoh , 
1977). 

The number fraction of atoms in a given ionization state 
Fi is governed by the following differential equation: 

1 dF- 

--7^ = ai^i{T)Fi_i- [ai{T) +Ri-i{T)] Fi+Ri{T)F.+i , 

(44) 

with o:, (r) being the ionization rate for a given temper- 
ature, and ion /, and the recombination rateF*! 

For CIE the ion fraction of a given state remains con- 
stant, thus dFi/dt ~ 0. For NEI dFi/dt ^ 0, and the ion- 
ization fractions have to be solved using the coupled dif- 
ferential equations of all ion species i as a function of time, 
or as indicated by Eq. |44]as a function of n^. To compli- 
cate matters, in reality neither kT^, nor «e are expected to 
be constant in time. The parameter n^t is often referred to 
as the ionization age of the plasma. 

Equation|44]can be solved by direct integration, but this 
is CPU intensive and not very practical when it comes 
to fitting X-rav spectral data. A faster approach (Hughes 
and He lfand , 1985; Kaastra and Jansen, 1993; Smith and 
Hughes, l2010t) is to rewrite Eq.|44]in matrix notation. 



1 dY 

«e dt 



= A.F, 



(45) 



determine the eigenvalues and eigenvectors of A and then 
solve for the uncoupled equations 



1 dY' 

«e dt 



:A-F', 



(46) 



with X the diagonal matrix containing the eigenvalue, and 
F' = V^'F, with V containing the eigenvectors. 



"For this equation it is assumed that the fractions only increase by 
ionization from a lower ionization state or recombination from a higher 
ionization state. This may not be stiictly true, if one takes into account 
inner shell ionizations, which may result in multiple electrons being lib- 
erated. 
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Figure 1 1 : The effects of non-equilibration ionization (NEl) illustrated for oxygen. Both panels look very similar, but the 
top panel shows the oxygen ionization fraction as a function of electron temperature for coUisional ionization equilibrium 
(CIE), whereas the bottom panel shows the ioni zation fraction as function of rip t, and for a fixed temperature of kTs= 
1 .5 keV (ionization/recombination rates based on ShuU and van Steenberd, 1982 ). Approximate, mean nj values for the 
plasma in the young SNRs Cas A, Tycho and SN 1006 are indicated. (Figure earUer pubUshed in .Vink. 2006.) 
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The main effect of NEI in young SNRs is that the ion- 
ization states at a given temperature are lower than in the 
CIE situation. In Fig.[TT]the effect of NEI in young SNRs 
is illustrated. As an example consider the presence of O 
VIII in Fig. [nil, which peaks in CIE at kT^ = 0.2 keV. 
Seeing a plasma in a SNR with most of the oxygen in O 
VIII may lead to the conclusion that the plasma is rela- 
tively cool. However, as Fig.lllb shows, the temperature 
may very well be kT^ — 1 keV, but with a low ionization pa- 
rameter of n^t « 8 X 10^ cm^^s. For temperatures relevant 
for young SNRs {kT^ = 0.5 — 4 keV) the time needed to 
reach CIE is around n^t w lO'^ eV. Typical values for rigt 
in SNRs are between 10^ — 5 x lO'^ cm^^s, i.e. the plas- 
mas of most SNRs are not in CIE. Note that the relevant 
time scales for reaching CIE is similar to that for reaching 
electron-proton temperature equilibration (Sect. l53T l. 

Once a plasma has reached CIE one can have a reversal 
of the NEI situation; i.e. the plasma can become overion- 
ized instead of underionized, if the cooling rate is faster 
than the recombination rate. This will be discussed in con- 
nection to mixed-morphology SNRs in Sect. 110.31 

NEI models as found in X-ray spectral codes can have 
various levels of refinement, both concerni ng the atomic 
data (a strong point of the SPEX NEI code Kaastra et al 



6.8 



2003h . or in the additional physics taken into account, 
such as non-equilibration of electron and ion tempera- 
tures (Sect. 15.21 ). and the ensuing gradients in the elec- 
tron temperature b ehind the shock front (e.g . the vpshock 

model i n XSPEC, Borkowski et all 2001a). Some mod - 

I ' II II ^ ' ^ 1 1 

els (e.g . Hamilton et al..'1 983tlKaastra and Jansenlll993l : 

Borkow ski et al., 2001a) also take into account the temper- 



ature structure within a SNR in the Sedov-Taylor phase of 
its evolution (Sect.|4|l. 



6.1.3 Line emission diagnostics 

In SNRs line emission results from collisional excitation 
of ions, a process dominated by electron-ion collisions. 
This is either direct excitation or through recombination, 
which usually results in an ion in an excited energy level. 
Since the density is very low, most ions can be assumed 
to be in the ground state; and, once an excitation occurs, 
collisional de-excitation or further excitation or ionization 
can be neglected. This also means that the ionization bal- 
ance can be t reated indep endently of the line emission 
properties (seelMewe, 1999i, for a full treatment). 

An important aspect of line emission that is character- 
istic for NEI plasma (a nd unde r certa i n circumstances 
photo -ionized plasmas: iLie dahl 1 19991 : iKallman et al 



20041) is inner shell ionization. In this process an electron 
from an inner shell, for example the K-shell (n = 1) is re- 
moved, whereas higher level atomic shells (L- or M-shells, 
« = 2,n = 3) are still filled. This results in an ion with a 
hole in the K-shell, but plenty of electrons around from L, 
or M shells. The ion can then either de-excite without ra- 
diation, by filling the hole, and using the energy from that 
transition for further ionization (called Auger transitions), 
or it can adjust radiatively, which is called fluorescence. 
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Figure 12: Fe-K shell line emission ener gies, as deter- 



rgii 

mined theoretically (squares, Fe II-Fe XVII. IPalmeri et al 



2003b iMendoza et al.L 120041) and observa tionally (trian- 
gles, Fe XVin-Fe XXV. lBeiersdorfer et al.Lri993i1 . 



The likelihood for a radiative transition is called the fluo- 
rescence yield, and is higher for atoms with larger nuclear 
charge. For example the fluorescence yield for a K-shell 
ionization of neutral iron (Fe I) for a K-shell transition 
is 34%, for neutral silicon it is 5% and for neutral oxy - 
gen it is only 0.8% dKrausel [l979l : iKallman et alll2004i) . 
Because the radiative transition occurs in the presence of 
L-shell and M-shell electrons, the effective charge of the 
nucleus is reduced as compared to a K-shell transition in 
the helium-like state. As a result the Ka line energy is 
smaller, and rises slowly as a function of ionization state 

(Fig.inj. 

Iron line emission is an important diagnostic tool for 
the state of a SNR plasma, i.e. the average electron tem- 
perature and ionization age, n^t. This is even true for 
medium energy resolution spectroscopy as provided with 
the CCD instruments on board Chandra, XMM-Newton, 
and Suzaku. Because of its high fluorescence yield and 
high abundance, Fe K-shell emission can be observed for 
all ionization states of iron, provided that the electron tem- 
perature is high enough (^T^e > 2 keV). The average line 
energy of the Fe-K shell emission provides information 
about the dominant ionization state (Fig. fT2l) . For ioniza- 
tion states from Fe I to Fe XVII the average Fe-K shell 
line is close to 6.4 keV. 

Fe has also prominent Fe-L-shell transitions in the 0.7- 
1.12 keV range. These transitions occur for ionization 
states with electrons present in the L-shell {n = 2), i.e. Fe 
XVII to Fe XXIV, with each ionization state having its 
own specific line transitions, which increase on average 
in line energy for higher ionization states. Fe-L-line emis- 
sion occurs for lower temperatures/ionization ages than 
Fe-K, i.e. kT^ > 0.15 keV. Taken together, Fe-L- and Fe- 
K-shell emission can be used to accurately determine the 
ionization state of the plasma. For example, Fe-K emis- 
sion around 6.4 keV could be caused by Fe XVII-XIX, 
or by lower ionization states, but in the latter case no 
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Figure 13: Simplified Grotrian diagram of He-like ions, 
with the transitions fro m an excited n = 2 state to the 
ground state (taken from Porquet et al. [ l200lh . The triplet 
consists of the resonant line r or w, the intercombination 
line, i — x + y, and the forbidden line / or z- 



Fe-L emission should be present. The presence of Fe-K- 
line emission around 6.7 keV indicates the presence of 
Fe XXV (He-like Fe), and around 6.96 keV Fe XXVI 
(H-like). High resolution X-ray spectroscopy greatly im- 
proves the value of Fe-line diagnostics, as it allows one 
to resolve the Fe-L shell emission in individual lines. And 
for Fe-K shell emission it is possible to detect different 
ionization stage individually instead of relying line emis- 
sion centroids. 

It is worth mentioning that Ka line emission from low 
ionization states of iron can be the result of dust grains em- 
bedde d in hot SNR plasmas (iBorkowski and SzvmkowiakL 
1997h . Hot electrons can penetrate dust grains, giving rise 
to inner shell ionizations inside the grains, whereas the 
emitted photon can escape from small grains. In addition, 
dust grains are slowly destroyed in hot plasmas due to dust 
sputtering, on a time scale of ~ lO'^/we s. This results in a 
slow release of near neutral iron into the hot plasma. This 
should give rise to the presence of a broad range of ioniza- 
tion stages of Fe inside the plasma. 

Another line diagnostic, available for high resolution 
spectroscopy {E /AE > 200), is provided by the "triplet" 
line emission from He-like ions, consisting of three promi- 
nent (blended) lines designated w (resonance), x + y (in- 
tercombination), and z (forbidden) (Fig. fT3l see Gabriel 
and Joi-dan, 1969; Mewe, 1999; Liedahl, 1999; Porquet 

(/ 



et al., 201 



|The ratio G= {z + x + y)/w = {f + i)/r 
is of interest for SNRs, as it is sensitive to the tempera- 
ture and ionization state of Jhejplasma (Mewe and Schri- 
jver. 



19781 ; lLiedahill999t) . For example, inner shell ion- 



ization of the Li-like state, which has ground state config- 
uration 15^25, will result in an excited He-like ion in the 
ls2s^Si state. Since de-excitation to the ground state is 
forbidden for electric dipole radiation (A/ = 0), this gives 
rise to forbidden, magnetic dipole, radiation. Thus, inner 
shell excitation enhances the forbidden line transition, and 



acts as a measure for the fraction of Li-like ions. This can 
be seen in Fig. [T3] showing that the G-ratio is high for 
He? < 3 X 10^ cm^^s for oxygen and ngf < 2 x lO"' cm^^s 
for silicon. This figure also shows that for higher «ef val- 
ues G-ratio may fall below the CIE values. The reason is 
that Li-like ions are no longer present, but at the same time 
the recombination rate to the He-like state is low, because 
the atoms are still underionized, and collisions are more 
likely to result in ionizations than recombinations. Only 
when n^t> lO'^cm^^s are the G-ratios similar to those for 
CIE at the same temperature. 

Another potential influence on the G-ratio is resonant 
line scattering. Resonant line scattering only influences 
the intensity of the resonant line, as it scatters the reso- 
nant photon out of the line of sight. A very high G-ratio 
(i.e. relatively strong forbidden line emission) may there- 
fore help to diagnoze resonant line scattering, and its im- 
portance can be estimated, if rtet and kT^ can be estimated 
independently from the G-ratio. Note that, since resonant 
line scattering does not destroy the photon, from the other 
parts of the SNR, in particular from the edges, there may 
be enhanced resonant line emission. 

The ratio R = z/{x + y) = f/i is for CIE plasmas a sen- 
sitive diagnostic tool to measure electron densities, if the 
densities are in the range «e = 10^ — lO'^ cm^-', which is 
much higher than encountered in SNRs. For SNRs it is of 
more interest that the R-ratio is also sensitive to the ion- 
ization age, and to a lesser extent on temperature (Fig. 1141 
right). The reason is that the R-ratio in SNRs is mainly 
determined by inner shell ionization, which enhances the 
forbidden line emission for low «ef ■ Once the fraction of 
Li-like ions is negligible, and inner shell ionizations are 
therefore no longer important, the R-ratio becomes rela- 
tively flat. Note that the R-ratio may help to disentangle 
resonant line scattering effects from pure NEI effects; res- 
onant line scattering does not influence the R-ratio. 

6.2 Line emission associated with radioac- 
tivity 

Apart from collisional processes one may also have line 
emission caused by radioactivity. Radioactive elements 
are produced during the lifetime of stars, and, of partic- 
ular interest here, during the supernova explosion. The 
most important radioactive element is ^^Ni, which decays 
m 8.8 day£3 into ^^Co, which subsequently decays into 
f^Fe (see Table 16. 2] for the details, and Diehl and Timme4 
1998L for a review). This decay chain forms the dominant 
source of energy for the expanding supernova ejecta dur- 
ing the first year after the explosion. It has a significant 
imprint on the light curves of supernovae. Moreover, its 
final product, ^^Fe, is the most abundant Fe isotope in the 
Universe. A large part of the production comes from Type 
la supernovae (Sect. l2.2l l. which produce typically 0.6 M© 
per explosion. 



The word "triplet" denotes here the presence of three closely spaced 
lines, and is not used in the usual sense indicating a mixed quantum state. 



'^^The decay time T is the e-folding time, as opposed to the half-life 



'1/2 - 



■ Tin 2. 
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Figure 14: Left: The G-ratios G = (/ + /) / r for He-like oxygen and silicon. Three different temperatures are show 0.5 keV 
(solid line), 1 keV (dotted) and 2 keV (dashed). Right: Idem, but now for the R-ratios R = f /i. Note that these ratios are 
based on pure He-like ion lines only, in reality some blends may occur with L-like satellite lines (e.g. Vink et al., 2()03l). 
which tend to increase the G-ratio at low «ef values. In these figures the effects of blends with line emission from Li-like 
ions are not taken into account. These enhance the G -ratio and R-ratio at low n^t. (This figure is based on calculations by 
the author using the SPEX code. lKaastra et al.M2003h . 



An element that is less abundantly produced is ^^Ti, 
with typical expected yields in core co llapse super- 
novae of 10"^ - 10""* Mq (e.g. IPrantzosl [2011). Type 
la supernovae probably produce less, unless the explo- 
sion is triggered by an explosion at the surface of 
the white dwarf (the so-c alled doub l e deto nation sub- 
Chandrasekhar model, e.g. Fink et al. . 2010h . '*'*Ti has a 
mucMonger decav time t han ^^Ni. namely 85 vr ( Ahmad 
et al., 2006). This makes it of interest also for the SNR 
phase. ^^Ti is an alpha-rich freeze-out product, which 
means that its yield is determined by the amount of a- 
particles left o ver after the i nitial stages of explosive nu- 
clear burning ( ArnetlL 1996h . It is therefore sensitive to 
the speed of expansion of the inner layers of the ejecta, as 
this determines how rapidly the density drops, which pre- 
vents further build up to more massive elements like ^^Ni. 
In addition, the '*'*Ti yield is sensitive to the mass cut (the 
boundary between what is accreted onto the proto-neutron 
star and what w i ll be e jected) and explosion asymmetries 
(INagataki et al. .11998'). This sensitivity to explosion con- 
ditions explains the astrophysical interest in '*'*Ti. 

The most unambiguous, direct signature of radioactiv- 
ity are the y-ray lines associated with their decay (as pre- 
dicted by Clavton and Silk, 1969). The daughter prod- 
ucts of radioactive elements are usually not in their nu- 
clear ground states, but rapid de-excitation occurs through 
the emission of 7-ray lines (Table |6^ . The nuclear decay 
lines of ^^Co at 847 keV and 1238 keV have been detected 
from SN 1987 A by several balloon and satellite experi- 
ments (Cook et al., 1988; Mahoney et al., 1988; Sandie 
et al.. 1 19881; i Teegarden et all il989i see also Sect. |9^ . 
Kurfess et al.l(ll992h later report the detection of 122 keV 
line emission from ^^Co by the Oriented Scintillation 
Spectrometer Experiment (CGRO-OSSE) on board the 
Compton Gamma Ray Observatory (CGRO). No super- 
nova Type la has exploded near enough in modern times to 
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Figure 15: The decay chain of ^""^Ti with the indicated chan- 
nels, electron capture (EC) and beta-decay. 



firmly detect 7-ray lines, but there was a hint of ^^Co emis- 
sion in CGRO-Comp tel observations of t he bright Type la 



supernova SN1991T (IMorris et al.Lll997h . which had a dis 



tance of about 13.5 Mpc. 

Apart from 7-ray line emission, the decay chain of ^'^Ti 
also results in hard X-ray line emission at 67.9 keV and 
78.4 keV, which are caused by the nuclear de-excitation 
of ^^Sc (Fig. [TSl) . This line emission, and the the 1157 
keV 7-ray line of ^Ca, have been detected for the 
young SNR Cas A (Sect. lU) with the CGRO-Comptel, 
Beppo-SA X-¥BS, and INTEGRAL-imS i nstruments (re- 
spect ivelv. [ivudin et allll994l;l'Vink et allE oOl; Renaud 



et al., '2006b'). The flux per line of (2.5 ± 0.3) 10"= cm 
s^'cm"^ impl ies a rather large ^^ Ti yield of (1.6 ± 
0.3) 10""^ Mr:, (lRenaud et all l2006bl) . This yield is com- 
parable to what has been inferred for SN 1987A, based on 
the late time light curve of SN 1987 A (iJerkstrand et al 
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Table 1 : Decay chains of the most important radioactive explosive nucleosynthesis products and some of their hard X-ray 
and 7-ray signatures. 









Decay time 


Process 


Lines (keV) 








8.8 d 


EC 


158,812 




56Co^ 


sope 


111.3d 


EC, e+ (19%) 


847, 1238 


5^Ni^ 






52 hr 


EC 


1370 




"Co^ 


"Fe 


390 d 


EC 


122 


44Ti^ 






86.0 yr 


EC 


67.9, 78.4 




'^^Sc^ 


44Ca 


5.7 hr 


e+,EC(l%) 


1157 



201 ih . A ^^Ti detection was also reported for the SNR 
RX J0852-4622 ("Vela jr"), based on measurements of the 



1 157 keV line with CGRO-Comptel dlvudin et al.L Il998h 



But INTEGRAL-imS measurements of the 68 keV and 
78 k eV lines are inconsistent with this detection ( Renaud 



2006ah . 



et al. 

The radioactive processes in the decay chains of ^^Ni, 
"Ni and '^^Ti are electron capture (EC) or beta-decay 
(e+ in Table 16.21 ). The beta-decay process is a source 
of positrons and may be partially responsible for the 
511 keV ele ctron-positron annihila tion hne emission from 
the Galaxy dPrantzos et al.L 1201 ih . Most of the positrons 
released by the decay of ^^Co probably do not survive the 
supernova phase, but in the subsequent SNR phase not 
many positrons may annihilate due to the low densities 
and high plasma temperature. No 5 1 1 keV line ernissiori 
has yet been detected from a SNR dKalemci et al. ■ l2006l : 



Martin etalll2010h . 



Apart from nuclear decay line emission, radioactivity 
also gives rise to X-ray line emission. The reason is that 
electron-captures results in a daughter product with a va- 
cancy in the K-shell. This vacancy may result in a K-shell 
transition that could be detected using more sensitive fu- 
ture X -r ay telescopes su ch as IXOIATHENA dBookbinder 



2010h 



LeisingI (2001) lists in total 17 radioactive ele- 
ments, but the most promising of them are ^^Fe (— >^^Mn, 
T = 3.9 yr, Ka energy 5.888 & 5.899 keV), "^^Ti (->44sc, 
T = 85 yr, 4.086 & 4.091 keV), ^^Ni (^^'^Co, t 108 kyr. 



6.915 & 6.930 keV) , and =^Mn (->^^Cr, T = 5.4 Myr, 
5.405 & 5.415 keV). 

For SNRs, Ka line emission of ^^Ti and ^^Ni is of par- 
ticular interest, since the decay times are compatible with 
the lifetimes of SNRs. For '*'*Ti the K-shell transitions 
form a complement to the nuclear decay line emission. 
The advantage of the K-shell transition measurements is 
that the spatial resolution of soft X-ray telescopes is in 
general superior to hard X-ray or 7-ray telescopes. But it 
should be noted that the line transitions are not unique. In 
principle, other Sc isotopes may also produce the K-shell 
lines. But stable Sc is not a very abundant element. A more 
severe problem is that the K-shell transitions are also a 
function of the ionization state of the atom (Sect. I6.1.3I I, 
which may not be a priori known, complicating the search 
for line emission. Nevertheless, future more sensitive tele- 



scopes are exp ected to detect these lines ( Leisinj, 2001 



Hughes et al.L 12009.) . And one may then be able to de- 



termine the ionization state of the radioactive elements, 
which allows one to position the "^^Ti in the shocked or 
unshocked parts of the eje cta. 

Borkowski et al. d2010h reported the detection of veloc- 
ity broadened Ka line emission from ''^^Ti (or technically 
speaking from its daughter product Sc) fo r the youngest 
known Galactic SNR G1.9H-0.3 (~ 100 vr lCarlton et al . 
201 ih . The measured flux of (0.35 - 2.4) x 10"^ ph 



cm s implies a Ti production of (1—7) x 10 M©. 
G1.9H-0.3 is possibly a Type la SNR (Sect.[8]l, but the ^^Ti 
yield is consistent with either a core collapse or Type la 
origin. A search for '* '*Ti associated K a line emission from 
Ca s A (iT^e iling an d LeisingLl2006h and SN 1987A ( Leis- 



ing, 



2006h did 



not yet result in detections. 



6.3 Non-thermal emission 

6.3.1 X-ray synchrotron radiation 

Relativistic electron s gyrating in a magnetic field em it syn- 
chrotron radiation (Ginzburg and Svrovatskiilll967l)ln In 



most cases encountered in astrophysics, synchrotron radi- 
ation is caused by a non-thermal population of electrons, 
which often well approximated over a large energy range 
by a power-law distribution, certainly if first order Fermi 
acceleration is responsible for the production of the rela- 
tivistic electron populations. 

Sources of synchrotron radiation usually produce this 
radiation over a very wide range in frequency, from the 
low frequency radio bands (^10 MHz) up to X-rays 
(^ 10'** Hz). Most SNRs are radio synchrotron sources, 
with synchrotron radiation coming from the SNR shell, 
but for composite SNRs also from the embedded pulsar 
wind nebulae X-ray synchrotron emission from SNRs 
was traditionally associated with the presence of a pul- 
sar wind nebula, but only relatively recently has X-ray 



'^For more massive charged particles synchrotron radiation is still 
emitted, but at a much lower efficiency. 

'^See for example the catalogue by Green D. A., 2009, 'A 
Catalogue of Galactic Supernova Remnants', Astrophysics Group, 
Cavendish Laboratory, Cambridge, United Kingdom (available at 
"hrefhttp://www.mrao. cam. ac.uk/surveys/snrs/" 
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synchr otron emission from young SNR shells been estab- 
lished (iKovama etal.[ll995L see Sect. lll.2b . 

A relativistic electron with energy E moving in a mag- 
netic field with strength B will emit synchrotron radiation 
with a maxi mum emission at the characteristi c frequency 
Vch given by ( Ginzburg and SvrovatskiiL 1965h 



Vch = 1.8 X 10'**Bx(-r^y Hz 



/zVch = 13.9 



lergy 

—]( ) keV, (47) 

lOO^G/ VlOOTeV/ 



with Bx ~ Y^2/3B the magnetic-field component perpen- 
dicular to the motion of the electron. For the typical mag- 
netic fields inside SNRs, B = 10 - 500 /xG, the emitting 
electrons have energies of 10-100 TeV. These electrons 
may therefore also be responsible for many of the de 



tected TeV gamma-ray emission (e.g. lAharonian et al 



200 IL |2004|). However, pion production by high energy 
ions also gives rise to TeV gamma-ray emission, and there 
is an on going debate on which radiation mechanism dom- 
inates the TeV emission fro m SNRs (see Sect. II 1.41 and 
Hinton and Hofmannll2009l for a review). 

An electron emitting synchrotron radiation suffers radia- 
tive losses with a rate: 



dE 
dt 



( E x2Bl 



(48) 



-4.05 X 10 



-7 



100 TeV/ VlOO^uG 



Bi 



ergs 



-1 



with Or the Thomson cross section. Usually one includes 
into B± a term to include the effect of losses due to inverse 
Compton scattering, replacing Bx by B^ff — B^^^ + B^, 
with Biad = 3 .3 jiG, the value for which the magnetic-field 
energy density equals the radiation energy density of the 
cosmic microwave background^ This corresponds to a 
synchrotron loss time scale of 



'^syn — 



dE /dt 



n 



= ^2.5 (,^) (^) yr. 

(49) 

For a population of relativistic electrons with a power- 
law energy distribution, the spectral index, a, of the syn- 
chrotron flux density spectrum is related to the spectral in- 
dex of the electron energy distribution, <7, as a = (<7 — 1 )/2. 
The typical spectral radio index of young SNRs is a « 0.6, 
indicating that electrons have an energy distribution with 
<7 w 2.2, close to the predictions from first order Fermi ac- 
celeration (Sect. l53T l. For X-ray emission it is more com- 
mon to use the photon index, F = a + 1 . 

Typically, X-ray synchrotron spectra of young SNRs 
have rather steep indices, F = 2 — 3.5, indicating a rather 
steep underlying electron energy distribution. This steep- 
ness indicates that the synchrotron X-ray emission is 
caused by electrons close to the maximum energy of the 
relativistic electron distribution. This maximum is usually 
for energies at which acceleration gains are comparable 



'« BlJ(&7t) = HCMB = 4.2 ; 



10 " erg cm ^. 



to the radiative losses, the so-called loss-limited case. Al- 
ternatively, the electron energy spectrum breaks off, be- 
cause the shock acceleration process has not acted for long 
enough time for the electrons to reac h higher energies . 
This is the so-called age-limited case (IRevnoldsL Il998h . 
In the age-limited case the underlying electron spectrum 
has an exponential cut-off in energy o= exp(— Zi/finiax), 
whereas in the loss-limited case the cut-off is expected to 
be steeper ex p(— i£'/£'iTiaTl^) ( Zirakashvili and Aharo- 
nian, .2007.) . Because the synchrotron emissivity function 
for a given electron energy has a rath er broad spectral pro- 
file (Ginzburg and Svrovatskii, 1965), the rather steep cut- 
off in electron energies results in a much more gradual roU 
over of the synchrotron spectrum: 



n(hv)<^(hv) "^expl-f , — \ ^ \ 



(50) 



The cut-off photon energy for the loss-limited case can 
be calculated by combining Eq. [41] and Eq. |49l which 



gives 



100 TeV 



:^0.32t7 



-1/2 



VIOO^G/ 



(51) 



5000 km s 



-1 



X4 



Inserting this in Eq. |47] shows that the cut-off photon en- 
ergy, /zVcut-off, is indepe ndent of the magnetic field ( Aha- 
ronian and Atovan. .l999f) : 



/iVcut-off = 1-47? 



5000 km s- 



keV. 



(52) 

See lZirakashvili and Aharonianl ( l2007h for a more detailed 
derivation. This equation indicates that X-ray synchrotron 
emission requires the relatively high shock velocities (> 
2000 km s^ ' ) only encountered in young SNRs. Moreover, 
the detection of X-ray synchrotron emission from SNR 
shells suggests that the diffusion coefficients cannot be too 
large compared to the Bohm-diffusion case (Sect. l53T l. i.e. 
J? ^ 10. 

6.3.2 Non-thermal bremsstrahlung 

Apart from synchrotron radiation, non-thermal X-ray 
emission can also result from lower energy electrons 
through bremsstrahlung and inverse Compton scattering. 
Inverse Compto n scattering is for SNRs important in GeV- 
TeV 7-ray band ( Hinton and HofmannL 12009) . but for the 
magnetic fields inside SNRs, B ft! 5 — 500 jxG, it is gener- 
ally not expected to be important in the soft X-ray band. 

Bremsstrahlung caused by the presence of a non- 
thermal electron distribution has been considered as a pos- 
sible sourc£_of2£:rax£ontinuun^ for some time 



(e.g.lAsvarov et al.Lll990l:IVink et allll997l:lBleeker et al 



20011 iLaminj, l2001b ). Since the electrons that produce 
X-ray bremsstrahlung have non-relativistic, to mildly rela- 
tivistic energies, identifying non-thermal bremsstrahlung 
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Figure 16: The fate of a suprathermal electron distribu- 
tion exchanging energy with thermal e lectrons thro ugh 
Coulomb collisions, as a function of n^t (IVinkL l2008bh . 



would help to obtain information about the low energy 
end of the electron cosmic-ray distribution. This relates 
to the poorly known electron cosmic-ray injection pro- 
cess (Sect. 15.5b . However, it is unlikely that non-thermal 
bremsstrahlung contributes enough radiation to identify 
it with the current generation of hard X-ray telescopes. 
The reason is that the Coulomb interactions between ther- 
mal and non-thermal electrons thermalizes the low en- 
ergy tail of the electrons on a relatively short timescale, 
net ~ lO'" cm^^s, so that non-thermal bremsstrahlung 
can only be e xpected in a n arrow region close the shock 
front (Fig.[T6] vin3.l2008bh . or for those SNRs that have 
an overall low ionization age. Such a low ionization age 
is for example encountered in SN 1006, but here the X- 
ray continuum is overwhelmed by synchrotron radiation 
(Sect. 111. 21 ). The low ionization age is caused by a low 
density of the plasma, which also makes that the brems- 
strahlung emissivity is not very large (Eg. 143b. 



7 X-ray spectroscopy with Chandra, 
XMM-Newton, and Suzaku 

Before reviewing the most recent results obtained by the 
current generation of X-ray missions capable of either 
high resolution spectroscopy and/or imaging spectroscopy, 
Chandra, XMM-Newton, and Suzaku, 1 describe some of 
the characteristics of these telescopes and their instru- 
ments in the context of studying SNRs. 

All the aforementioned X-ray telescopes have multi- 
ple instruments on board, with the CCD-detectors be- 
ing the work horses for all three missions. X-ray de- 
tecting CCDs were first used by t he Br oad Band X-ray 
Telecsope (BBXRT, 'Petre et all Il993h . flown for ten 
days on the space shuttle Col umbia, but on l y with the 



launch of the ASCA satellite jTanaka et all 1 19941) did 



the results of CCDs imaging spectroscopy become avail- 
able to the broad community of X-ray astronomers. X- 



ray CCD-detectors are suitable for imaging spectroscopy 
with medium spectral resolution (R = E/AE w 10 — 43 
FWHM for 0.5-6 keV) and imaging resolution close to 
the telescope resolution^ The Chandra telescope pro- 
vides the highest spatial resolution, which is sUghtly un- 
dersampled by the CCD pixel size of 0.49". The CCDs 
themselves come in various types, even within one instru- 
ment. The CCD-detectors used in current observatories 
go by the names ACIS (AXAF CCD Imaging Spectrome- 
ter) on board Chandra, EPIC (European Photon Imaging 
Camer a) behind the XMM- Newton telescopes, with EPIC- 



MOS (ITurner et all 1200 Ih beh ind two of the telesc opes. 



and EPIC-pn behind the third dStruder et all 1200 Ih . and 



XIS (X-ray Imaging Spectrometer) for Suzaku. XIS con- 
sists also of four separate detectors behind four indepen- 
dent telescopes. 

CCD-detectors are excellent for studying the spatial dis- 
tribution of elemental abundances inside SNRs, or identi- 
fying regions with likely synchrotron radiation contribu- 
tions. They cannot be used to resolve Fe-L complexes, or 
the He-like triplets (Sect. 16.1.3b . which offer more precise 
diagnostic power to measure temperatures and ionization 
parameters. They are also not good enough to measure 
Doppler shifts for velocities below a few thousand km s^ ' . 

In addition to CCD camera's, Chandra and XMM- 
Newton carry grating spectrometer j*^ which give high 
spectral resolution for point sources, but, compared to the 
CCD instruments, have a smaller effective area. 
XM M-Newton has a Reflectio n Grating Spectrometer 



(RGS, den Herder et al 



2001'), consisting of two grat- 
ings behind two of the X-ray three telescopes. It covers 
the soft part of the X-ray spectrum from 6-38 A (0.32-2 
keV). About 50% of the X-rays go through the gratings, 
and the remaining fall on the EPIC-MOS instruments. As 
a consequence, the RGS gathers data simultaneous with 
EPIC. The spectral reflection angle for the RGS is large 
(^ 4.8°) compared to those of the transmission spectrom- 
eters of Chandra. As a result an extended source in a 
monochromatic line gets distorted (squeezed) in one di- 
rection. The design of the RGS is such that the spectrum 
is stretched over a large detecting area, with the spectrum 
being detected by nine CCDs for both RGS 1 and RGS20 
The CCDs' intrinsic energy resolution are used to separate 
spectra of different spectral orders. The large physical size 
of the spectral extraction region has as a disadvantage that 
the background induced by energetic particles is relatively 
large. But it has the advantage that the RGS is less sensi- 
tive to the source extent, with the additional line width 



A predecessor of the CCD-detectors was the Solid State Spectrome- 
ter (SSS) on board the Einstein satellite, which did not provide spatial in- 
formation, but an energy resolution similar to present day CCD-detectors 
I Joyce et al.. 1978). 

Suzaku has another high resolution X-ray spectrometer on board, 
XRS, which is based on the caloiimetric detection method. Unfortunately 
this instrument failed shortly after the launch oi Suzaku. 

" Both for RGSl and RGS2 there is one CCD not working, luckily 
intended to cover different wavelengths. So the current number of CCDs 
is 8 per grating. 
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Energy (keV) Energy (keV) 

Figure 17: CCD {XMM-Newton-MOS) and Reflective Grating Spectrometer {XMM-Newton-RGS) spectra of an oxygen- 
rich (1E0102.2-7219) and a Type la SNR (0519-69.0). The figure iUustrates the gain in spectral resolution of the RGS 
instrument over CCD-detectors, even for mildly extended objects; both SNRs have an extent of about 30". The figure also 
illustrates the differences between Type la spectra and core collapse SNR spectra, with the former Type la SNR being 
dominated by Fe-L emission (see Sect.[8]and |9.1t . 
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G292.0+1.8 (East) 




Supernova remnants observed with HETGS 



Figure 18: The O VII triplet (Sect. |6T3] ) as observed 
with the XMM-Newton RGS in the oxygen-rich SNR 
G292.0+1.8 (Sect. |93J- The RGS still offers a good res- 
olution despite the fact that G292.0H-1.8 has an angular 
size of nearly 8'. In this case the presence of a bright, nar- 
row, bar-like feature helped to obtain a good spectral res- 
olution. The red line indicates the emission from the reso- 
nance line, magenta the intercombination line, and green 
the forbidden line. (Figu re first produced in Vink. Bleeker, 
Kaastra, and Rasmussen ( 20041) .) 



caused by source extent, A0, given by 
AA =0.124(A^/1') A. 



(53) 



So RGS spectra from typical young SNRs in the Magel- 
lanic Clouds, which have A0 sa 0.5', are only mildly dete- 
riorated compared to those of point sources with an addi- 
tional broadening of AA = 0.06 A. This is comparable to 
the intrinsic resolution of the RGS, AA = 0.06 A. Adding 
these numbers quadratically shows that the spectral resolu- 
tion for those SNRs is A /AA = E /AE w 265 at A = 22.5 A 
(i.e. around the O VII triplet, Sect. l6.1.3l l. Even for a source 
extent of 4', one still obtains a spectral resolution of 40, 
outperforming CCDs at long wavelengths and compara- 
ble to the spectral resolution of CCD-detectors at 6 keV. 
If the extended object is dominated by some particularly 
bright regions, the spectral resolution may be very good, 
as illustrated with the O VII triplet line emission of SNR 
G292.0H-1.8 in Fig. [18] The gain in spectral resolution ob- 
tained by using the RGS compared to the CCD data is 
illustrated in Fig.fTT] 

The Chandra transmission gratings come in two vari- 
eties, the Hi gh Energy Trans r nissiq n Grating Spectrome- 



ter (HETGS. ICanizares et al.Ll2005h and the Low Energy 



Transmission Grating Spectrometer (LETGS), which pro- 
vide spectral coverage in the range 1 .5-30 A (HETGS) and 
1.2-175 A (LETGS). The gratings for both instruments 
can be put into the optical path of the mirror assembly. The 
HETGS consists of two grating arrays, a Medium Energy 
Grating (MEG) and High Energy Grating (HEG),placed at 
a small angle with respect to another This results in an X- 
shaped dispersion pattern on the detectors, with the zero 




Figure 19: An overvie w of Chandra hi^h energy gratings 
observations of SNRs (IDewevl 120021) . (Figure kindly pro- 
vided by D. Dewey.) 



order image in the center of the X. Since they are trans- 
mission gratings, the HETGS and LETGS have plus and 
minus order spectra. Only one order is shown in Fig. [19] 
The detectors onto which the spectra are dispersed are the 
normal imaging detectors of Chandra, i.e. either the ACIS 
(CCD) instrument or the High Resolution Camera (HRC), 
a microchannel plate detector. For the HETGS the ACIS-S 
is the instrument of choice as it allows for the separation of 
different spectral orders. The LETGS is either used with 
the ACIS-S array, if the long wavelength range is not of in- 
terest, or with the HRC, which offers a longer wavelength 
range and smaller pixel size, but does not offer an intrinsic 
energy resolution and has lower quantum efficiency. 

The dispersion angles of the HETGS and LETGS are 
much smaller than those of the XMM-Newton-RGS. As a 
result lines are much more narrowly spaced on the detec- 
tor, with smaller detector pixels to match. However, for 
extended sources this means that the monochromatic im- 
age of a source in a given line will easily start overlap- 
ping with the images from nearby lines. This is even true 
for moderately extended sources, such as the already men- 
tioned Magellanic Cloud SNRs. In numbers: the degra- 
dation in spectral resolution is (Eg. l53Tl 0.67 A/', 1.33 A/', 
and 3.33A/ fo r the HEG, MEG and LETGS, respectively 
(lDewevLl2002h . The effect of source extent on the spectra 
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is illustrated for the HETGS in Fig. [19] which gives an im- 
pression of the challenges faced when analyzing HETGS 
SNR data. Note that the complications are more severe 
for Type la SNRs (Sect. [8]i, because they have bright Fe- 
L hnes, which, due to the richness of the hne emission, 
easily blend together for extended s ources. This is illus- 
trated for N103B in Fig.[T9](see also lDewevll2002h . The 
situation is much better for the spectrum of the core col- 
lapse SNR 1E0102. 2-7219, whose intrinsic spectrum is 
less complex (Fig. 117b. 

Given that small SNRs are already difficult to analyze, 
one would assume that the situation is almost hopeless 
for large SNRs. However, HETGS observations of a large 
SNR like Cas A (5') still provide important results. But in- 
stead of obtaining high resolution spectra of the SNR as a 
whole, one can obtain spectral properties of small. X-ray 
bright spots, embedded in a backgr ound of diffuse emis - 
sion from other parts of the SNR (iDewev et "ZL l2007h . 
In general, for extended sources, the RGS performs bet- 
ter for obtaining high resolution spectra of the SNR as a 
whole, whereas the HETGS will provide spectral informa- 
tion that relates to spectral properties in localized regions 
of the SNR. 

The problems of obtaining high resolution X-ray spectra 
of extended sources will be over once the next generation 
of spectrometers will be put into space. These will in all 
likelihood be based on the calorimetric principle; an X- 
ray photon heats a small metal element, and the increase 
and subsequent decrease in temperature is measured with 
high sensitivity, using, for example. Transition Edge Sen- 
sors (TESs). TESs are superconducting devices operated 
near the critical temperature, where a small change in tem- 
perature leads to a strong response in conductivity. The 
spectral resolution can be as good as A/i = 2 eV at 6 keV. 
The small elements can be built into arrays, thereby allow- 
ing for imaging, high resolution spectroscopy. Currently, 
microcalor imeters are planned for the Japanese/US Astro- 
H mission (iTakahashi et al. and the proposed ESA 

mission Athenal^*! 



8 Type la supernova remnants 



X-ray spectroscopy of Type la super- 
nova remnants 



8.1 



As discussed in Sect. 12.21 thermonuclear or Type la SNRs 
produce about a factor ten more mass in Fe-group ele- 
ments than core collapse supernovae. This is the reason 
that Type la SNRs are often characterized by strong Fe-L- 
shell emission, which cannot be resolved with CCD instru- 
ments, but it can be partially resolve d by the current grat - 
ing instruments (Sect.|7] Fig. [17] and iHughes et al.Lll995h . 
Many SNRs are now identified as Type la remnants based 
on the prominence of the Fe-L emission (see Table [2] for 
a list of likely Type la SNRs and relevant references). In 



the case of two SNRs, Tycho/SN1572 (Fig.[T]l and SNR 
B0509-67.5 the Type la identificatio n has been confirmed 
by optical spectra of the light echoes dKrause et al l l2008bl: 



Rest et al.L 12008. truly spectacular results). 



Not all Type la SNRs display prominent Fe emission. 
For example the X-ray spectra of SN 1006 show ha rdly 



any evidence for Fe-L emission (IVink et all l2003r) . al 



though recently Fe-K e mission has been r e ported based 



on Suzaku observations (lYamaguchi et al.L 120071 120081) . 



For SN 1006 the lack of Fe-L emission is partially due 
to th e low ionization age (w^ffti 2 x 10^ cm~^s. e.g. Vink 
et al., '2003*), which makes that the dominant Fe ioniza- 
tion stage is below that of Fe XVII. In that case hardly 
any Fe-L line emission is produced, but some Fe-K fluo- 
resence line emission is still to be expected (Sect. I6.1.3I I. 
However, optical/UV spectroscopy of a bright UV star be- 
hind SN 1006 also shows a lack of absorption features 
from cool, unshocked, Fe (Hamilt on et al.Ll2007h . This is 
somewhat discomforting, but it may indicate that SN 1006 
produced a relatively large amount of intermediate mass 
elements and relatively small amounts of iron-group ele- 
ments. This would sti ll fit in with the opt ical results for 
Type la supernovae of iMazzali et al.l (l2007h . which allows 
for the occurrence of fainter, Fe/Ni-poor supernovae. It is, 
however, more difficult to expla in the lack of Fe with the 
models of IWooslev et al. I (l2007h that suggest that the total 
mass of all Fe-group elements does not vary much among 
Type la SNe, but that there is only a variation in the the 
relative amounts of radioactive ^^Ni and stable Fe. This is 
clearly of interest, and both models and the observational 
data should be scrutinized further 

Another reason why a Type la SNR may not be very 
bright in Fe-L or Fe-K line emission is that the Fe-rich 
layers has not yet been completely shock heated by the 
reverse shock. For example, Tycho/SN 1572, although far 
from devoid of Fe-L and Fe-K emission, shows less Fe-L 
emission with r espect to Si XIII emission th an most other 
Tvp e la SNRs (iHwang and Gottheld.ll997l: Decourchelle 
et al.. l200ll:lHughes et al.lll99 iK 

The effects of the progression of the reverse shock into 
the ejecta is nicely illustrated using the Large Magellanic 
Clould (LMC) Type la SNRs in Fig. [20] For SNR B0509- 
67.5 ( Warren and Hug hes, 2004; Kosenko et al. 



Badenes et al. ■ l2008bh the Fe-L emission is confined to 



a small, partial shell that is most pronounced in the west. 
This is the youngest LMC remnant (Rest et al., 2005|), 
and the one with the largest X-ray hne Doppler broaden- 
ing, Ov ~ 5000 km ' , as measured using XMM-Newton- 
RGS high resolution spectra ( Kosenko et al., 2008). These 
spectra also indicate that the Fe-L emissio n is mainly com' 



ing from Fe XVII. In SNR B05 19-69.0 dKosenko et al 



Athen a is a simplified ve rsion of tfie International X-ray Observa- 
tory (IXO. lBookbindeAl2010l) 



[2OIO) and N1 03B tLewis et all 120031; van der Heyden 
et al., l2002bh the Fe-L line emitting shell is more pro- 
nounced. The XMM-Newton-RGS specti-a of B05 19-69.0 
show less line broadening, Oy ~ 1900 kmi^', and Fe-L 
emission is now also coming from the higher Fe ioniza- 
tion stages Fe XVII-Fe XXI. The Chandra data show that 
the lower Fe ionization stage, Fe XVII, is located closer to 
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Table 2: Likely Type la supernova remnants 



Name 


Location 


Distance 
(kpc) 


Radius 
(pc) 


High Res. Spec. 


Other studies 


Gl.9+1.3^ 


MW 


8(?) 


2 


(1,2,3) 




Kepler (SN 1604/G4.5+6.8) 


MW 


6 (4,5,6) 


3 




(7,8,9,10,5) 


Tycho (SN 1572/G 120. 1+1.4) 


MW 


3 


3.8 


(11,12) 


(13,14,15,16,12) 


SN 1006 (G327.6+14.6) 


MW 


2.0(17) 


9.3 


(18,19)) 


(20,21,22,23,24,25) 


B0509-67.5 


LMC 


50 (25) 


3.6 


(27) 


(28,29) 


N103B 


LMC 


50 


3.6 


(30) 


(31) 


BOS 19-69.0 


LMC 


50 


3.9 


(32) 




DEM L 71 


LMC 


50 


8.6 


(33) 


(34) 


B0548-70.4 


LMC 


50 


12.5 




(35) 


DEM L3 16 A" 


LMC 


50 


15 




(36,37) 


B0534-69.9 


LMC 


50 


16 




(35) 


DEM L238 


LMC 


50 


21 




(38) 


DEM L249 


LMC 


50 


23 




(38) 


B0454-67.2 


LMC 


50 


27 




(39) 


IKT 5 (B0047-73.5) 


SMC 


60 (26) 


15" 




(40) 


IKT 25 (BO 104-72.3) 


SMC 


60 


18 




(40,41) 


DEM S 128 (B0103-72.4) 


SMC 


60 


26 




(38) 



MW=Milky Way, LMC=Large Magellanic Cloud. SMC=Small Small Magellanic Cloud. 

References: 1) Reynolds et al. ( 2008); 2) Borkowski et al. ( 2010) 3) Carlton et al. ( 201 1 );4) Sankrit 'etaiH200l: 5)rwiigi2008cl): 6^IChiotellis et aUlMl: 7) Kinugasa and 

Tsunemi : 1999); 8 Cassam-Chenai et al. :2004); 9) Badenes et al. (2007); 10) Reynolds et al. (2007); ll)IGhavamian et al.lf2001l); 12)'Havato et al.' (2010^; 13)Hwang and 

Gotthelf ( 1997); 14) Decourchelle et al. : 2001 ); 15) Hwang et al. (2002); 16) Tamagawa et al. ( 2009); 17) Ghavamian et al. ( 2002); 18) Vink et al. ( 2003); 19) Vink ( 2004a); 

20).Long et al. (2003); 21) Yamaguchi et al. (2007); 22) Acero et al. (2007); 23) Yamaguchi et al. (2008); 24) Katsuda et al. (2009); 25) Miceli et al. (2009); 26) Sch aefeti 
20081); 27) Kosenko et al. (2008); 28) Warren and Hughes : 2004); 29) Badenes et al. (2008b);30) van der Heyden et al. (2002a); 31) Lewis et al. (2003); 32) Kosenk o eta!] 
2010) ; 33) van der HeydenetaL (_2003); 34) Hughes et al. ( 2003a); 35) Hendrick et al. ( 2003^; 36) .Nishiuchi et aLaOOl,) ; 37) , Williams and Chaa005.) ; 38) ,BorkowskiTtal] 
2003) ; 39) ISeward et'an ' iBoog) r40) lvan der Heyden et alj 120041) ; 41) lLee etalllMl li) 

Notes; 

^ Type la origin debatable, and largely based on morphology 'Borkowski et al.V2010|). 
Only the Fe-L emitting interior is visible in X-rays, so the radius is not necessarily the shock radius for IKT 5. 

^ D^M^^^^^a^^^^^ssociate^S^^^She^^^^iveiyt^^bundaM and 

Chu T200^r The^wTS?rarniMvno^^ 
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Figure 20: Chandra X-ray images of LMC Type la SNRs. These are three-color images, and in all cases red represents 
the O VII/O VIII ba nd (0.5-0.7 keV), the color green represen ts Fe-L emission line 1 ke V). From left to righ t, top to 
bot tom: B0509-67.5 dWarre n and Hughes, '20 04|), B0519-69.0 jKosenko et al.l l2010b. N103B jLewis et al.Ll2003h . Dem L 
71 (iHughes et al 1 l2003ah . and B0534-699 ( Hendricket al.L l2003h . The order of the figures is approximately indicative of 
the relative dynamical age. (Images generated by the author using the Chandra archive.) 
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the center. This is expected, becaus e this is the loca t ion of 
the most recently shock-heated Fe ( Kosenko et al. . 2010L 
see Fig.lSTb. 

Co ncerning the more mature SNRs Dem L7 1 ( Hughes 



et al. I2003al: Ivan der Hevden et all l2003h and B0534- 
69.9 (* Hendrick etalll2003h . one observes that Fe-L emis- 
sion is coming from the whole center of the SNR, in- 
dicating that (almost) all ejecta have been shock-heated. 
There is a clear indication from recent X-ray data from 
both Galactic and LMC SNRs that the shocked ejecta 
of Type la remnants are well stratified; from the outside 
in one first encounters oxygen, then a silicon-rich layer, 
and finally the iron-group elements (Fig. 1211 1. testifying 
of a more or less spherically symmetric explosion. This 
is cle arly the case in SNR BOS 19-69.0 ("Kosenko et all 
2010h and in the southweste rn region of Kepler's SNR 
dCassam-Chenai et al. . l2004 . whereas the more mature 



SNR Dem L7 1 has an Fe-rich core (Ivan der Hevden et al., 
2002bl: iHughes et all l2003ah. Onl v Tvcho's SNR seems 



less stratified dDecourchelle et all 1200 Ih . Note that the 



radial X-ray line emission is also influenced by ioniza- 
tion age differences in different layers, as indicated by 
ASCA and X MM-Newton observations of Tvcho ( Hwang 
and Gotflielf. 1 19971: 'Decourchelle et all I2OOII). Kepler 
dCassam-Chenai' et al., .2004) , and the already discussed 
SNR BOS 19-69.0. 

The more or less layered structure is unlike that of core 
collapse SNRs (Sect. 19. Il l, where in some cases plasma 
with more massive elements seems to have overtaken 
plasma with less massive elements. In addition, the mor- 
phologies of Type la SNRs are more symmetric than those 
of core collapse remnants (more on this in Sect.sec:orich). 
The reason for this may be a more spherically symmet- 
ric explosion and that the cuxumstellar medium (CSM) 
has been less disturbed, possibly because Type la SNe oc- 
cur more frequently outside turbulent starforming regions. 
Nevertheless some imprints of CSM interaction can be 
found in Type la SNR (see next section). 

Both for core collapse and for Type la SNRs the varia- 
tion of elemental abundances inside the remnant poses a 
problem in modeling the overall X-ray spectra of SNRs. 
The reason is that one cannot realistically expect to fit the 
spectra with one or two NEI (Sect. 16.1.21 ) plasma models. 
On the other hand, adding many more NEI components 
leads to unconstrained results. Imaging spectroscopy of- 
fers some help, as one can extract spectra from distinct re- 
gions and apply single or double NEI components to these 
individual spectra. But projections effects play a role as 
well. Recently Kosenko et al. (2010) stretched the capa- 
bihties of multi-component fitting of SNR BOS 19-69.0 to 
the maximum, by using seven NEI components, each cor- 
responding to a layer that could also be identified in the 
stratification inferred from narrow band Chandra imag- 
ing (Fig. I2TI) : an O-rich layer, a Si/S-rich layer, a Ar/Ca 
layer, a low ionized Fe layer, and a highly ionized Fe layer 
The continuum component was constrained by setting the 
abundance pattern to that of the LMC. Such a procedure 
works better for the radially stratified Type la SNRs than 




4 5 6 

Energy (keV) 

Figure 23: Detail of the Suzaku spectrum of Ty- 
cho/SNlS72, showing Ka emission from Cr, Mn and Fe 
(reprinted from .Tamagawa et al...2009.) . 



for the more chaotic spatial abundance distribution of core 
coflapse SNRs. 

An alternative method that has already given interest- 
ing results is to model the whole hydrodynamical evolu- 
tion of the SNR, from the explosion to the SNR phase, 
including the temperature/ionization evolution and X-ray 
emission, using the theoretically predicted initial strati- 
fication and velocitv distribu tion of the ejecta (Badenes 
et al., 2003, 200S, 2 008h:,Sorokina et al.i..2004:iKosenko . 



I2OO6: Kosenko et al.l l2008l ] 201 ih . This procedure works 



relatively well for Type la SNRs, because the explosion 
properties of thermonuclear supernovae are much better 
constrained and allow for less variation than for core col- 
lapse supernovae. The grid of explosion models for Type 
la SNe is, therefore, much more limited. The method has 



been applied to Tycho /SN 1S72 (Fig. [22] B adenes et al 



2006t Kosenko, 2006) and SNR BOS 19-6 9.0 (Badenes 



et al..|2006 : .Kosenkai2006) . For the latter. .Badenes et al 
(l2008bl) found that the most energetic Type la explo- 
sion models, which synthesize more Fe, fit the data best. 
Kosenko et al ] (12008 ) reported for this SNR a very high 



plasma velocity based on the Doppler broadening mea- 
sured with the XMM-Newton-RGS instrument. In addition, 
the RGS detected N VII line emission, which helps to con- 
strain the amount of shocked CSM, because nitrogen is 
not a Type la nucleosynthesis product. In this case the in- 
ferred circumstellar density is «h = 0.4 — 0.8 cm^-'. 

An uncertainty in modeling the combined 
hydrodynamics/X-ray emission from Type la SNRs 
is the influence of efficient cosmic-ray accelerati on, as 
discussed in Sect.l4land lS.SllKosenko et al.l (l2008l) incor- 
porated the cosmic -ray acceleration effects by altering 
the equation of state of the shocked plasma, which gave 
a good fit to the XMM-Newton (EPIC-MOS) spectra for 
a cosmic-ray en ergy density that is about 40% of the gas 



densi ty (see also lDecourchelle et al. [ |2000l :l iKosenko et al 
201 ll) . 



Abundance studies of SNRs in X-ray are usually con- 
fined to the alpha-elements (C,0, Ne, Mg, Si, S, Ar, Ca) 
and Fe-group elements (Fe/Ni). These elements are the 
most abundant elements in the Galaxy, and are the dom- 
inant nucleosynthesis products of supernovae (Sect. |2]i. 
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Figure 21: The stratification of O VII- VIII, Si XII, and Fe XVII-XVIII (Fe low), Fe XIX-XXI (Fe high) in SNR BOS 19- 
69.0, based on Chandra data. The left panel shows the radial surface emi ssion profile, whereas the right panel shows the 
deprojected (emissivity) profile. (Figure taken from Kosenko et al. . 2010h . 
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Figure 22: XMM-NewtonMOS>\ spectrum of west ern part of Tycho/SN 1572. The two panels show two models (solid 
lines) from an extensive grid of Type la models (iBadenes et al. . l2006h . The model characteristics are indicated: the 
explosion model (here DDTc and DDTc), the circumstellar density, Pam in units of 10^^^ g cm^^ and absorption column 
density, A^h in units of 10^' cm^^. (Courtesy C. Badenes) 
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and all have prominent line transitions in the X-ray band. 
Other elements, for example those with uneven atom num- 
ber, ar e usually too rare to detect in X-rays. However, re- 
cently iTamagawa et al. I (I2OO9,) reported the detection of 
chromium (Cr, Z = 24) and manganese (Mn, Z = 25) K- 
a line emission from Tycho's SN R using Suzaku data 



(Fig. |23] see also lYang et al.L 12009.) . The nucleosynthe 



sis of Mn is influenced by the presence of the neutron 
rich element ^^Ne, whose abundance is determined by the 
metallicity of the white dwarf progenitor An explanation 
for the relatively high flux ratio of Mn/Cr is, therefore, 
that SN1572 had a white dwarf progenitor that started as 
a main sequence star w ith solar or supersolar metallicity 
(IBadenes et al.ll2008ah . In other words, SN 1572 was not 
the explosion of a Population II star. 



Evidence for Type la progenitor im- 
prints on the circumstellar medium 



8.2 



As explained in Sect. 12.21 there is little doubt that Type la 
supernovae are thermonuclear explosions of white dwarfs, 
but there is much uncertainty about the binary system that 
leads to the growth of a C/O white dwarf close to the Chan- 
drasekhar limit. The mass transfer from the companion of 
the white dwarf can occur during various phases in the evo- 
lution of the companion; during the main sequence, red 
giant (RG), or asymptotic giant branch (AGB) phase, or 
when the companion itself is a white dwarf. The mode by 
which the mass transfer occurs may shape the CSM of the 
system. For example, Roche-lobe overflow is considered 
to lead to unstable mass transfer, but it may be stabilized 
by a fa st, ^ 1000 krn s^^, wind em anating from the white 
dwarf oHachisu et al. [ Il996lll999h . This fast wind will cre- 
ated a large (> 10 pc), low density bubble surrounded by 
a dense shell. Evolved stars, like RG and AGB stars, have 
slow winds 10 km 5^') which result in a rather dense 
CSM with a density profile scaling as r^^ (Sect.|4]i. 

The specific characteristics of the binary are, therefore, 
expected to have an important imprint on the CSM, and, as 
a result, on the evolution, morphology and X-ray spectra 
of the SNR. Models for Type la usually ignore the possi- 
ble effects of pre-supernova evolution, and instead assume 
a uniforni CSM (i.e. an 5 = model, in the notation of 
Chevalieil(ll982h . see Sect.©. 

One important constraint on the progenitor system pro- 
vided by SNRs comes from optical observations. These 
show that many of the SNRs mentione d in TablelUhave so- 



called B aimer-domi nated shocks (e.g. lSmith et al.L 11991 



Sankrit et all l2005h . Optical spectra from their shock re- 
gions show Ha line emission, without [Nil] line. This is 
caused by collisional excitation of neutral hydrogen en- 
tering the shock region. In addition, broad Ha emission 
arises from c harge transfer of neutra l hydrogen with shock 
heated ions (IChevalier et al. . 119801: Ivan Adelsberg et al 



Balmer-dominated shocks are Type la SNRs. The rea- 
son for this is that progenitors of core co llapse supernovae 
create large Stromgenspheres (> 10 pc IChevaheii[l990l) . 
And even if they explode as red supergiants the bright 
flash of ionizing photons when the s hock reaches the 
surface, creates a large ionized region (iChevahen. l2005h . 
Type la supernovae, on the other, do not produce a UV 
flash at shock break-out bright enough to ionize a large 
region. 

The presence of neutral hydrogen within a couple of 
parsecs from the explosion, therefore, is easy to explain 
for Type la SNR, provided that the progenitor itself was 
not a powerful source of ionizing photons. This excludes, 
therefore, an important class of Type la candidate progen- 
itors, namely the so-called supersoft sources. These are 
soft X-ray sources that probably consist of a white dwarfs 
accreting material from a companion at a rate close to the 
Eddington limit, and bu rn the accreted material in stable 
way at the surface (Li and van den Heuvelll 19971) . But be- 
cause supersoft sources are luminous UV/X-ray sources, 
10^^ — 10^** ergs^^ they will easily ionize a region of 
30 pc in radius. Many of the Balmer-dominated Type la 
SNRs, including Tycho, Kepler, B0509-67.5, B05 19-69, 
DEML 71, can , therefore, not have had supersoft sources 
as progenitors (Ghavamian et al.. 2003i) . 

There are several other recent studies that cast a light 
on the structure of Tvpe la SNRs. In one of them Badenes 
et al. (2007) investigates whether Galactic and LMC Type 
la SNRs show any evidence for having evolved inside 
large, tenuous cavities, as implied by the fast-wind model 
of Hachisu et al, ( 1996). This study shows that this is not 
the case, as the radii of the SNRs are too small and their 
ionization age, n^t, too large compared to model predic- 
tions that incorporate the effects of fast progenitor winds. 
Possible exceptions are SN 1006 and RCW 86. However, 
SN 1006 is located high above the Galactic plane, where 
the density is expected to be low irrespective of the pro- 
genitor model, and RCW 86 is not a bona fide Type la 
SNR, as it is located in an OB association (IWesterlund . 
il969) . in a region with a high star forming rate, as indi- 
cated by outbreaks of g as from the Galactic plane region 
dMatsunaga et al. [ |200lh FI 

Another study that may be revealing something about 
Type la progenitor systems concerns an analysis of 
ChandralXMM-Newton of two o ld, ~ lO'^ yr, LMC SNRs , 
DEM L238 and DEM L249 by iBorkowski et all (l2006h . 
These two SNRs are rather faint in X-rays, but are both 
characterized by a relatively bright interior dominated by 
Fe-L emission, somewhat similar to DEM L7 1 and SNR 
B0534-699 displayed in Fig. |20] The most surprising re- 
sult of the analysis is that of the spectra from the inte- 
rior indicate high values for the ionisation parameter, nj. 



20081) . The presence of Balmer-dominated shocks is only 
possible when the neutral fraction in the CSM is relatively 
high. It is therefore of interest that most of the SNRs with 



^'The only exceptions are the large SNRs the "Cygnus Loop" and 
RCW 86, which are both likely core collapse SNRs I Ghavam ian et all 
I2OOII) . but see WilUams et al. 1 20 11 ) for an alternative view on RCW 86. 

^^However, recently W illiams et al,. 1,201 1.) argued, nevertheless, for a 
Type la origin, because there is no evidence for a neutron star inside the 
SNR, and there is evidence for large amounts of Fe. 
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Figure 25: Results of the hydro dynamical study of Ke- 
pler's SNR lChiotellis et aL ( 201 1 ). showing the expansion 
parameter in the SNR for a model in which the progeni- 
tor systems moved with respect to the interstellar medium, 
and created a one side shell due to interaction of the pro- 
genitor wind and the ISM. The expansion parameter in the 
densest regions of the s hell are as low as i3 = 0.35, com- 



parable to obs ervations ( Dickel et ami988HVinkll2008ct 
Katsuda et al.ll2008c) . (Figure provided by A. Chiotellis.) 



- lO'^cm-^s for DEM L238 and - 4 x lO^'cm^^s for 
DEM L249. Such high values require relatively high den- 
sities in the interior. This is in contrast to the low densities 
implied by the X-ray spectra of the shell, and, indeed, the 
overall faintness of these SNRs in X-ravs. lBorkowski et al. 
§006) argue that in order to account for the high central 
density the blast wave must have encountered a high den- 
sity CSM early during the SNR evolution. Such a dense 
medium is usually a signature of a slow wind outflow, in- 
stea d of the fast winds im plied by the models of Hachisu 
et al. I.I999I), or a simple, uniform density medium. 

Finally, evidence that winds from Type la progenitors 
affect the CSM is provided by the historical SNR of SN 
1604 (also known as Kepler's SNR, or Kepler for short). 
Kepler is located high above the Galactic plane, with lat- 
itude b = 6.8°, which transla tes to 446 pc for a distance 
of 4 kpc d Sankrit et all l2005h FI This high latitude would 
make Kepler an obvious Type la SNR, if it were not for 
the presence of dense, nitrogen-rich material, mostly in 
the northern part of the remn ant, which indicates the pres- 
ence of a circumstellar shell (Ivan den Bergh and Kampei , 



Or 670 pc for a distance of 6 kpc, the distance preferred based on 
mod elin g the expIosionjtf^he^inetic_energ^k^^O^^_e^_^A 



1977 : Bandieral 1987 ). During the 1990s a core collapse 
origin for SN 1604, was t herefore preferred. X-ray spec - 
troscopy, first with ASCA ("Kinugasa and Tsunemi','l999') 
and later with XMM-Newton (Cassam-Chenai et al.. .2004) 



and Chandra (iRevnolds et al.L 120071) . revealed that the 



spectrum is dominated by Fe-L emission around 1 keV, 
and bright Fe-K shell emission around 6.7 keV, character- 
istic for Type la SNRs. Additional evidence for a Type la 
origin is the low oxygen content, and the lack of an ob- 
vious X-ray point source that may be the neutron star that 
should in most cases re sult from a core collapse supernova 
dRevnolds et al.Ll2007 '). 



The idea that the SNR encounters a dense shell in the 



north is reinforced by two recent X-ray studies dVink , 
2008d :' Katsuda et al., 20 08cr). w hich confirm an earlier ra- 
dio measurement (.Dickel et all [l988) that the expansion 
of this SNR in the north is relatively slow with a small 
value for the expansion parameter (Eq. |5]l of j3 sa 0.35, 
whereas the average value is « 0.5. The value for the north 
is unusual in that it falls below the value 0.4 for a SNR in 
the Sedov phas e of is evolution. On the basis of expansion 
measurements, IVinkI (l2008c,) estimated that the mass in 
the swept-up, nitrogen-rich shell must have been ^ 1 M©. 

This result, and the morphology of Kepler, resembling 
a bow shoc k in the north (Fig.l24la nd Fig.|25]for a model), 
is used by Ichiotellis et al. ( 201 lb to argue that the pro- 
genitor system was a binary consisting of a C/O white 
dwarf with an asymptotic giant branch (AGB) star as a 
companion. Such a system would be labeled as a symbi- 
otic binary. Part of the slow wind (^15 km 5^') of the 
AGB star was accreted onto the white dwarf, until it ex- 
ploded. In order to account for the nitrogen enrichment of 
the shell one has to assume a main sequence mass for the 
secondary star of 4-5 Mq . The bow-shock shape of Kepler, 
and its height above the Galactic plane (^ 400 — 700 pc), 
can be explained by a proper motion of the system with 
about 250 km^~^(c.f. Bandiera and van den B ergh. 199ll : 
iBorkowski et al. . 19941) . Nevertheless, within a Type la 
scenario, it is not quite clear, why a binary system was 
ejected with such a high velocity, as usually interaction 
of two binaries leads to the ejection of fast single stars 
([Leonard and Duncanl I1990I) . ■ Jut it does occur, see for ex- 
ample the relatively high velocity ( ^ 130 km^~') of th e 



AGB star Mira A and its companion (iMartin et al. L l2007h 



mod elin g the explosionj_^^he^inetic_enere^^s; 
et al.' l2008a : ,Vinl^i2008d : IChiotellis et alifeoilh . 



9 Core collapse supernova remnants 
9.1 Oxygen-rich supernova remnants 

As explained in Sect. |2] the nucleosynthesis yields of core 
collapse supernovae peaks for the element oxygen. The 
oxygen yield itself is strongly correlated with the main 
sequence mass of the progenitor, although for masses 
> 30 Mq the poorly known fraction of material that falls 
back onto the central object (neutron star or black hole) 
may reverse this trend. An important class of SNRs are 
the so-called oxygen-rich SNRs (Table [3]). As the name 
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100 




Energy (keV) 



Figure 24: Left: Chandra image of Kepler's SNR, with red indicating Si-K a emission (1.75-1.95 keV), green Fe- 
L emission (0.8-1.6 keV), a nd blue continuum emission (4-6 keV). The image is based on a deep, 750 ks, Chandra 



observation (IRevnolds et al.L 12007) . Right: raw count rate spectrum from Kepler's SNR, based on the same observation. 



Note the dominant Fe-L emission between 0.7-1.5 keV (c.f. Fig.fTTb). 



Table 3: Oxygen-rich supernova remnants 
Name Location Distance Radius High Res. Spec. Other studies 







(kpc) 


(pc) 






Cas A (Gil 1.7-2.1) 


MW 


3.4(1) 


2.6 (2) 


(3,4) 


(5,6,7,8,9) 


MSH 11-54 (G292.0H-1.8) 


MW 


>6(10) 


>7.5 


(11) 


(12,13) 


Puppis A (G260.43.4) 


MW 


2.2(14) 


17.6 


(15) 


(16,17) 


B0540-69.3 


LMC 


50(18) 


8.4 


(19) 


(20,21) 


N132D 


LMC 


50 


11 


(22,23) 


(24) 


IE 0102.2-7219 (IKT 22) 


SMC 


60(18) 


5.9 


(25,26) 


(27) 


B0049-73.6(IKT6) 


SMC 


60 


21.2 




(28,29) 


B0103-72.6(IKT23) 


SMC 


60 


26.1 


(28) 


(30) 



MW=Milky Way, LMC=Large Magellanic Cloud, SMC=Small Small Magellanic Cloud. 

Selected references: 

l)jReed et al. : 1995); 2)'Gotthelf et al.' JlOOit); 3'Bleeker et alJ V2m%; 4)|Lazendic et alJ (Jopl; 5)IWillingale et alj l2003); 6)IHwang and Lamins) f2O03l); 7)IHwang et alj 
(|004); 8) Maeda et al. ( 2009); 9;iDeLaney et al. 1 2010); 10) Gaensler and Wallace : 2003); 11) Vink et al. ( 2004); 12) Gonzalez and Safi-Harb ( 2003); 13) Park et al. : 2007); 
14)Reyno,soetal. (1995); 15) Winkler et al. 1 1981); 16) Hwang et al. ( 2008); 17) Katsuda et al. ( 2008a); 1 8) Schaefet i2008); 19). v an de r Heyden et al.. i,200U; 20) Hwang 
et al. (2001); 21) Park et al. (2010); 22) Behar et al. (2001); 23) Canizares et al. (2001); 24) Borkowski et all boOTj) ; 25) lRasmusse'n"et alj 120011) ; 26) IFlanagan et aljfSOOi) ; 
27) ,Gaetz et al, 1200(1) ; 28) . van der Heyden et al.. i2004.) ; 29) .Hendrick et al.. a005-) ; 30) .Park et al. ( .2003i) 
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suggests, these SNRs show large overabundances of oxy- 
gen, indicating that these are remnants of the most mas- 
sive stars {Mms ?J 18 M©). Initially these sources were 
mainly identified based on t heir optical oxygen- line emis- 



sion pro perties ("rOIin. e.g. Goss et all 19791: Ma thew- 
son et al.. Il980l: Laskeil ll979l: lDopita et al.1 Il98ll) . but 



in many cases the X-ray spectral properties confirm the 
elevated oxygen abundances, in particular for young rem - 
nants like Cas A (Vinket al.', 1996; Willingale et al., 2002*), 
1E01 02.2-7219 (.Hughes and Singh, 1994; Gaetz et al 



2OOOI: iRasmussen et all l20oT[ Flanagan et al.L l2004 . and 
G292.0H-1.8 (MSH 11-54, e.g.l Parket al. '.'2007S. Even for 
some of the more mature oxygen-rich SNRs, like N132D 
and I'uppis A, it is still possible to identify regions with 



oxygen-ric h ejecta ( Behar et al 
E007 : Katsuda et alll2008al) ^ 



2001: Borkowski et al 



The neon and magnesium yields of core collapse su- 
pernovae trace the oxygen yield (Fig. O, one therefore 
expects also that the Ne and Mg abundances of oxygen- 
rich SNRs are high. This appears to be indeed the case 
for young oxv gen-rich SNRs like IE 0102.2-72 19 (Ras- 
mussen et al., .2001; Gaetz et al., 2000. Fig. [TTl l26l l28Tl 
and G292.0+1.8"(lHug , hes et al.[l200ll: Gonza lez and Safi- 
Harb.l 2003l:|Park et all 120071: 1 Vink et al.ll2004 ). But an im- 
portant exception to this rule is Cas A, which has Ne/Mg 
abun dances a factor 5-10 below the expected values (" Vink 
et al.7T99^ In general, the most striking feature of the 
X-ray spectrum of Cas A is the very bright silicon and sul- 
phur emission. It is usually assumed that the Cas A ejecta 
are dominated b y oxygen, which is almost fully ionized 
(IVink et al.l 1996 : Hwang and Laming, 2 003 )a.but Dewey 
et al. (l2007h has argued that the most abundant element in 
the ejecta is silicon. 

If the ejecta indeed are almost devoid of hydrogen, this 
greatly affects the ejecta mass estimates ("Sect. l6.ll Eq.l43Ti. 
For example, for Cas A it implies a rather low eje cta mass 
of 2-4 Mp„ with about 1-2 of oxygen (.Vink et al.l 



19961: IWillingale et al.Ll2003h . This oxygen mass is consis- 
tent with a progenitor main sequence mass of 18 ± 2M0. 
This led Chevalier and Oishi (2003) to the conclusion that 
Cas A must have been a Type lib SNR, given the simi- 
larity between the apparent mass loss and main sequence 
mass of the Type lib supernova SN1993J and Cas A. This 
conclusion has been spectacularly confirmed by the opti- 
cal light echo s pectrum of the Cas A supernova (" Krause 
et al.^008a). 

Another property of Cas A is the strong bipolarity of the 
silicon/sulphur-rich eiecta, as indicated bv optical (Fesen 
et al.l2006). X-ray (V ink, 2004b; H wang et al., 2004), and 
infrared observation (iHines et al. [ l2004h . These observa- 
tions show the presence of Si/S-rich material outside the 
radius of the main shock in the northeast and southwest. 
The presence of ejecta outside the shock in the northeast 
was already known for a long time, and its optical feature 
has long been referred to as "the jet". In X-rays the pres- 
ence of two opposite jets is best shown by mapping the 
ratio of Si XI II over Mg XI Ka emission ("Vinkl \ 



of the Si XIII emission /Hwang et all I2OO0I : lYang et al, 



2008) 



The hydrodynamical simulations by Schure et al.l ( 120081) 
show that these jets cannot have survived the interaction 
with a circumstellar shell. Such a shell would be present 
if there has been a fast Wolf-Rayet star wind ploughing 
through a red super giant wind (Sect. I4.3l l. This strongly 
suggests that the Cas A progenitor exploded while in the 
red supergiant phase. This is consistent with a main se- 
quence mass below ~ 25 M©, for which stars do not go 
through a Wolf-Rayet star phase. Since there was likely no 
high mass loss Wolf-Rayet star phase, the low ejecta mass 
is best explained with a binary star scenario, in which a 
high mass loss is caused by a common enve lope phase 



(lYoung et all 120061 : Ivan Veelen et all l2009h . A similar 



scenario was propose d for the Type lib supernova 1993 J 
(IWooslev et all 1 19941) . whose optical spectrum closely re- 
sembles that of the l ight echo spectrum of Cas A (Krause 
et al., 200 8a). There is no candidate secondary star in 
Cas A, so the common envelope phase may have led to 
a merger of the cores of the two stars. 

The apparent bipolarity of the Cas A explosion may be 
related to the mass loss history of the progenitor, which 
may have included a common envelope phase. It is not 
clear whether all core collapse supernovae develop strong 
bipolarities either due to the accretion process onto the 
collapsing core, or perhaps due to rapid rotation (Sect.|2]i. 
In more massive progenitors, these bipolarities may be 
smothered while the explosion ploughs through the ex- 
tended envelope of the star. Indeed, for supernovae the 
amount of asphericity indicated by optical polarimetry 
seems to depend on how much of the hy drogen-rich en- 
velope was removed (I Wheeler et alll2002h . It may also be 



related with the common envelope merger, as it may have 
increased the rotation of the stellar core. 

A puzzling property of the jets in Cas A is their sili- 
con richness. If the jets really come from the core one 
expects the jets to contain more iron-rich plasma. Cas A 
does contain many iron-rich knots, but not so much in the 
jets, as well as a t the outer edge of the southeastern part of 
the main shell (iHughes et all I2OO0I: iHwang and Lamingl 



20031) : see the blueish colored regions in Fig.l26fa. The pro 



2004b 



jected radius of these h' on knots, together with th e inferred 
age of Cas A 330 yr (IThorstensen et al. [ l200lh . indicate 
that the mean velocity of these iron knots are higher than 
7000 km ' . This is very high inde ed, as core collapse 
simulations bv lKifonidis et al. [ (l2006h show iron-group el- 
ements with velocities only up to 3300 km 5^'. 

The morphology of Cas A as seen in Fe-line emission 
does not have an obvious axis of symmetry. As Fig. l26h . 
shows iron-rich regions can be found in the southeast, the 
north, inside the silicon-rich shell, and in the west. Al- 
though in the northern part the Fe-rich ejecta seem to be 
on the inside of the Si-rich shell, Doppler measurements 
using XMM-Newton EPIC-MOS data show that in this re- 
gion the velocity of iron is high er, ^ 2500 km5~\ than 
that of silicon, - 2000 kmi"' dWillingale et all l2002h . 



Hwang et al. , I2OO4 Fig. |26h), or the equivalent width It is therefore likely that also in the northern part the Fe 
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Figure 26: Multi-color Chandra X-ray images of four oxygen-rich SNRs (see Table[3]for the sizes). Top left: Cas A, with 
red i ndicating O VIII Lv-i3 emission, green Si XIII Hea. and blue Fe XXVI Hea (based on the 1 Ms observation Hwang 



et al.. 120041) . Top right: A different view of Cas A, based on a combination of radio (blue), X-ray Si XIII emission (green), 
and the ratio of Si XIII over Mg XI emission (red). The la tter brings out the jet - counter- jet system, which protrudes the 
radio/X-ray shell in the northeast-southwest region dVinld i2004b: Hwang et all l2004l : ISchure et al.L 12008.1 . Bottom left: 
G292.0H-1.8, with red O VIII Ly-g emission, green is Ne XLya, and blue Si XIII HEa emission. Image based on a 516 ks 
long observation dPark et alll2O07h . Bottom right: The Small Magellanic Cloud SNR IE 0102.2-7219, with red O VII and 
O VIII emission, green Ne IX and Ne X emission, and blue emission above 1 .27 keV, which includes Mg XI, Mg XII, and 
Si XVIII emission (see also Fig.[T7]andl28Tl. 
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Figure 27: Left: Doppler maps of Cas A in the bands of He-like silicon, sulphur and iron, based on XMM-Newton-MOS 
data. Right: Reconstructed "side view" of Cas A, based on the Doppler maps, using a combination of north-south position 
and Doppler velocity to reconstruct the third dimension. Red represent silicon, green si Ucon, and blue iron. No te that in 
the North iron is in front of the silicon-rich layer. Both figures are based on figures from lWillingale et al.l (l2003h . 



rich ejecta are outsi de the Si-rich shell. The study b v Will- 
ingale et al. (i2002h also showed, by combining Doppler 
information with the spatial structure, that the overall mor- 
phology of the Si-rich resembles more closely a donut-like 
structure than a spherical shell (Fig. |27] ). as already sug- 
gested by Einstein Fo cal Plane Crystal Spectrometer da ta 



(IMarkertet all 1 19831) and ASCA data (IHolt et al.Lll994l) . 



Finally, it is worth pointing out that Cas A is special in 
the sense that the detection of ^'^Ti 7-ray and hard X-ray 
lines from this SNR suggest a rather large production of 
this element (1.6 x lO^'^Mrri. see also Sect. 16.21 Renaud 
et al.. l2006bh This is much higher than model predictions, 
but consistent with what has been inferred for SN 1987 A. 
This higher than expected yield may be re lated to the fact 



that t he Cas A explosion was asymmetric dNagataki et al 
1998h . Somethi ng that ni a y also apply to SN 1987A. On 



the other hand, iPrantzosI (1201 L) shows that the solar ra- 
tio of '^^Ca/^^Fe, the stable products of, respectively, ^^Ti 
and ^^Ni, is more consistent with the inferred yields for 
Cas A and SN 1987A than with model predictions. This 
suggests that asymmetric explosions are more common 
than assumed, or that the production processes of these 
elements is not yet well understood. 

High resolution X-ray spectroscopy with the Chandra 
HETGS suggest that the SMC SN R IE 0102.2-7219, ha s 
a similar, "donut-like" structure dFlanagan et al.l l2004t) . 



This has been established using a special property of trans- 
mission gratings, namely that the plus and minus orders 
are each others mirrors in terms of wavelengths, whereas 
the projection of the object in a single wavelength is not 
mirrored between the plus and minus orders. If for exam- 
ple one side of the object has on average a redshift and 
the other side a blueshift, this results in the blueshifted 



side to be shifted toward the center (zeroth order) and the 
redshifted side to be shifted away from the center The 
monochromatic images in the plus or minus orders are, 
therefore, squeezed in one order and a stretched in the op- 
posite order 1E0102. 2-7219 shows exactly this behavior 
as shown in Fig.|28]for the Ne Ly-a line. 

If we now turn the attention to another well studied 
oxygen-rich SNR G292.0H-1.8 (MSH 11- 54) one also sees 



a non -spherical distribution of elements (IPark et al.L 12002 



20071 Fig. 126b). There is the enigmatic bar, running from 



east to west, first seen in X-ray im ages taken with the Ein- 
stein satellite ( Tuohv et al. . 1982h . The origin of this struc- 
ture is not so clear, but may be a pre-existing feature in the 
CSM that was overtaken by the forward shock. The idea 
that the bar is of CSM origin is supported by the lack of ev- 
idence for enhanced metal abundances , as determined by 



Chandra X-ray imaging spectroscopy ( iPark et al.Ll2002r) . 



The overall metal abundance distribution shows a strik- 
ing asymmetry between the northwest, which contains 
more Si-group elements, and the southwest, which con- 
tains more oxygen, neon, magnesium. There is no evi- 
dence that this asymmetry is caused by an asymmetry in 
the CSM. T hus an asymmetric explosion seems to have 



taken place (IPark et al.L l2007h . iGonzalez and Safi-Harb 



(I2OO31) analyzed Chandra X-ray spectra and came to the 
conclusion that the abundance pattern matches that of a 
star with a main sequence mass of 30-40M,3. 

If there is one thing that one that can be concluded from 
X-ray observations of core-collapse SNRs, it is that their 
explosions are very irregular, with their SNRs showing 
large deviations from spherical symmetry and a lack of 
elemental stratification. Material from the core (Fe-group 
elements) can be ejected with speeds in excess of those 
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Figure 28: C/ianfira-HETGS spectrum of the oxygen-rich 
SNR 1E0102.2-7219. The top panel shows part of the 
spectrum in the negative order (m = — 1) as seen with the 
Medium Energy Grating of the HETGS. The SNR in in- 
dividual emission lines shows up as images ordered by 
wavelength. Bottom panel: a comparison between Ne X 
Lya images in the negative, zeroth order, and positive or- 
der Differences in radial velocities over the face of the 
SNR result in distortions of the image with respect to the 
zeroth order images. Deviations from spherical symmetry 
show up as distortions of the images that are different in 
the negative and positiv e orders, as can be c learly seen 
here. (Figure taken from lFlanaean et aL ( 2004 ).') 
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Figure 29: Results from a statistical analysis of the Si 
XllI emission images (around 1.85 keV) as observed with 
the Chandra-AClS instrument. The numbers Pq,P2,P3 are 
linked to different symmetry properties of the SNR emis- 
sivity distribution, with Pj indicating an overall deviation 
from circular symmetry, whereas P2 is sensitive to mirror- 
symmetry. Type la SNRs (red) clearly occupy a differ- 
ent region of the d iagram than core co llapse SNRs (blue). 
(Figure taken from lLopez et al.l (l2009h .) 



of layers higher up in the star. This behavior is markedly 
different fro m that of the well stratified Type la SNRs. 



Recently, iLopez et al.l (120091 1201 Ih made a quantita 



tive analysis of the morphological differences between 
core collapse and Type la SNRs, using a a multipole sta- 
tistical analysis of the X-ray surface brightness distribu- 
tion of SNRs. Deviations from circular symmetry show 
up as large values for P3 /Pq (normalized third moment), 
whereas a mirror symmetry gives large values for P2/Po- 
As Fig. |29] shows. Type la SNRs are on average rather 
more spherical than core-collapse SNRs; their mean val- 
ues for both f 2 / Pq and f 3 /Pq are smaller. In addition, there 
is the trend that if an individual Type la SNR shows devi- 
ations from spherical symmetry (large P3/P0) it tends to 
have less mirror symmetry. These difference makes that 
Type la and core collapse SNRs occupy different areas in 
the P3 / Pq versus P2 /Pq diagram. 

That there is a structural difference between core col- 
lapse and Type la SNRs is not too surprising, since the 
energy of Type la SNRs comes from explosive nucleosyn- 
thesis, and requires a regular burning front to propagate 
through the white dwarf. Core collapse supernovae are 
driven by a the gravitational energy from the collapsing 
core, and much of that energy is deposited in the inner re- 
gions. Irregularities will not impair the explosion of the 
rest of the star. 

However, there are some signs of aspherical structures 
that may reveal us something about the explosion pro- 
cess: there are the donut-Uke morphologies of Cas A and 
1E0102. 2-7219, and, in particular, the bipolar Si-rich jets 
in Cas A. Currently there is no theory to explain these 
features, but it is tempting to relate the jets of Cas A to 
those of long duration gamma-ray bursts. In fact, the less 
energetic gamma-ray bursts, the X-ray flashes, have prop- 
erties rather similar to what has been infer red for Cas A . 
For examp le, for the X-ray flash SN 2006aj ("Maz zali et al. , 
[2007; L aming et al.l 2006) the inferred explosion energy 
is 2 X 10^' erg, and the inferred main sequence mass is 
~ 20 M0 with an ejecta mass of ~ 2 Mq. 

9.2 The X-ray emission from SN 1987A 

The Type 11 supernova SN 1987 A was discovered on 
February 23, 1987 in the LMC. Due its proximity, it re- 
mains, to date, the best studied supernova. The detection 
of ~ 20 neutri nos during a short interval of ^ 13 s by 



Kamiokande 11 dHirata etal.[ll987h and 1MB dBionta et al. 



1987h confirmed beautifully the theory that Type 11 super- 
novae are the result of the collapse of a stellar core into 
a neutron star, and that most of the energy is released 
through a cooling population of neutrinos. The detection 
of neutrinos also meant that the collapse of the core was 
very accurately timed, and that the core collapse preceded 



1989 



the optical detections by about 3 hr (see lArnettet al. 
for a review). 

In many respects SN 1987 A was a peculiar supernova. 
Its maximum luminosity was unusually low, correspond- 
ing to Mg = —15.5, compared to ^ — 18 for typical Type 
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II supernovae. This low luminosity is the result of the com- 
pactness of the progenitor star at the moment of explosion: 
the progenitor was not a red supergiant, but a B3 I su- 
pergiant, identified as the star Sanduleak -69°202. Such 
an evolved, blue giant is expected to have a fast, tenu- 
ous wind, creating a low density cavity within a red super- 
giant wind of a previous mass loss phase. Instead, optical 
spectroscopy dFransson et al. ,ll989|) provided evidence for 
dense (we ^ 3 x 10^ cm^^) circumstellar material enriched 
in nitrogen. L ater optical imaging revealed the presence 



ring itself has considerable small-scale structure, with re- 
gions of differe nt densities, and fingers sticking inward 
(IMcCravi 120071 Fig. [3TJ. As a result the ring is not light- 
ing up all at onc e, but several "hot spo ts" turned on; first 
in the northeast (iBurrows et al. ■ l2000h . later also in other 
regions. By now all along th e ring one can find X-ray emit- 
ting regions (iRacusin et all 12009^ see Fig. [311). 

The complexity of the interaction of the blast wave 
with the dense inner ring do not fit into the simplified de- 
scription of SNR hydrodynamics as sketched in Sect. [4] 



which was later c onfirmed by a, now ico nic, Hubble Space iMcCravl l2007h . The basic structure of a blast wave heat- 



Telescope image ( IBurrows et al.Lll995h . The fact that the 



of three rings dCrotts et all ll989l : IWampler et all Il990h . dBorkowski and SzvmkowiakL 1 19971 : iMichael et all 12002 



progenitor e xploded as a blue su pergiant, the chemical 
abundances ('Podsiadlowski, 1992 ). and the triple-ring sys- 
tem may be explained by the merging of a 15-2OM0 star 
with a 5 M© star, about 20,000 yr prior to explosion; see 
the hydrodynamical simulations of Morris and Podsiad- 
lowski (2007, 2009). 

SN 1987A is of great importance, since it is the first 
supernova for which one has been able to follow the evo- 
lution from supernova to SNR. Moreover, our knowledge 
about the supernova type and progenitor is much more de- 
tailed than for any other SNR. One can say that the super- 
nova stage has given way to the SNR stage, when the emis- 
sion is no longer dominated by radiation from the cooling 
ejecta, but by radiation from shock-heated CSM/eiectarn 

A few months after the explosion SN 1987A was de- 
tected in hard X-rays ( Sunvaev et al. . 1987; Dotani et^ll 
19870 . The hard X-ray emission was probably the result 
of multiple Compton scatterings of 7-rays generated by 
the decay of radioactive material. Nuclear decay-line emis- 
sion from ^^Co was detected suprisingly early on, indi- 
cating considerable mixing of inner ej ecta toward the sur- 
face of the supernova envelope (e.g. iLeising and Sharj, 
19901 and references therein). Moreove r, the nuclear de- 
cay lines were redshifted (Sandie et al.^988; Teegarden 
et al.. ll989l : lTueller et al.U l990r This is an indication for 



an asymmetric explosion, bringing to mind the evidence 
for an asymmetric explosion of the Cas A supernova in 
the form of jets and high velocity iron knots (Sect. |9Tt . 

The onset of the SNR stage of SN1987A, as far as the 
X-ray emission is concerned, was marked by the gradual 
increase of the soft X-ray emissio n as seen early on by 
ROSAT, stalling in February 1991 dHasinger et allll996l) . 
Since then all major X-ray observatories have observed 
SN 1987 A on a regular basis, following the X-ray evolu- 
tion of the SNR in considerable detail (Fig.l30]l. 

The X-ray observations show that after an initial phase, 
in which the outer shock wave was moving through the, 
relatively tenuous, medium associated with the progeni- 
tor's fast wind, between 1995-1999 the shock wave started 
to interact with the central dense ring seen in the optical. 
This ring has a size of 1.7" by 1.2", consistent with a ring 
of radius 0. 19 pc seen under an angle of ^ 45°. The dense 



'^^ Note that stiictl^_s£e^ingjhisjiieans^hatradi^^ 
et al.^jipO?), should be considered SNRs, as their radio emission is due 
to interaction of the ejecta with the dense CSM. 



ing the CSM, and a reverse shock heating of the ejecta is 
still valid. But a range of shock velocities must be consid- 
ered, as the encounter of the blast wave with the density 
enhancements in the ring and its protrusions lead to a sys- 
tem of transmitted and reflected shocks. The transmitted 
shocks heat the material of the ring, whereas the reflected 
shocks go back into the plasma behind the blast wave. 
These reflected shocks are low Mach number shocks, as 
they go through already shock-heated plasma, and result 
in additional heating and compression of the akeady hot 
plasma. The transmitted shocks will have a range of ve- 
locities, depending on the density of the material and the 
obliquity with which the blast wave hits the protrusions. 
The slowest transmitted shocks (V^ < 200 kms) become 
radiative (Sect. 14.1b . resulting in the bright optical radia- 
tion that accounts for the bright knots that lie as beads all 
along the optical ring as o bserved by the Hubble Space 
Telescope dKirshner , l2007'.Fig.[3ni. 

Because the blast wave reached the ring in about 5-10 yr, 
the average blast wave velocity was initially ^ 15,000 — 
30,000 kmi^'. In contrast, recent expansion measure- 
ments indica te an exp a nsion velocity of ^ 4000 kmi~' 



in the r adio dNg et al 



X-rays dRacusin et al 



2008h and 1800 ±600 kms-' in 
^009). Until 2004 (6000 days af- 
ter the supernova) the expansion velocity in X-rays was 
similar to that in the radio, indicating that the blast wave 
reached the denser region of the inner rings around that 
time. The X-ray emission a lso became brighter around 



that time dPark et all 2005|). Currently , the light curve 
has flattened again dRacusin et aL , 2009L see Fig.l30]l. The 



strong deceleration of the blast wave from 15,000-30,000 
km.?^^ to 1800-5000 km means that a strong reverse 
shock must have developed with shock velocities in the 
frame of the ejecta of ^ 10,000 km 5^'. Note that the 
difference between the radio and X-ray expansion may 
be partially due to the fact that the X-ray emission de- 
pends on density as °^ n\ (Sect. 16. Il l, whereas the radio 
synchrotron emission scales probably as °^ rie (with some 
uncertainty, because synchrotron radiation also depends 
on the electron acceleration properties of the shock and 
on the magnetic-fie ld strength, see also the discussion in 
Zhekov etalll2010l) . The X-ray expansion may, therefore. 



be more skewed toward higher density regions, than the 
radio expansion measurements. 

From 2000 to 2010 SN 1987 A was observed several 
times with both the CCD and grating spectrometers on 
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Fi gure 30: Left: Th e evolution of the X-ray flux of SN 1987 A (iRacusin et al. [ l2Q09h . Right: Two frames taken from Fig. 1 
of Park et al. I (l2006h . showing Chandra images of the SN1987A, which indicate that it did not only substantially brighten, 
but also changed morphology as the shock heated more parts of the circumstellar ring. The first regions to be heated were 
in the northeast, but by 2009 most of the ring is lit up in X-rays, althou gh the brightest regi on is still the northeast. In 
addition the radius of X-ray bright torus has increased from 0.6" to 0.78"( Racusin et al. , 2009h . 



o 




Figure 31: Left: Composite optical/X-ray image (Credit: X-ray: NASA/CXC/PSU/S.Park & D.Burrows.; Optical: 
NASA/STScI/CfA/P.Challis). Right: Illustration of the interaction of the blast wave of SN1987A with the circumstel- 
lar ring. The blast wave shock heats the "fingers" on the inside of the ring, and as a result both in X-rays and in the optical 
the emission comes from discrete bright spots rather than smooth structures. (Source: NASA/CXO/M. Weiss) 
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board Chandra and XMM-Newton. The CCD spectra are 
well fitted with a two component model, a soft compo- 
nent with kTf^ = 0.2 — 0.3 keV, which is close to CIE, and 
a harder compo nent with kTg = 2 — 3.5 keV an d _ne£~ 
3 X 10^' cm-^s dPark et al.1 120061: IZhekov et all l2010h . 



natio ns from various studies are in good agreement ( Heng 



The soft CIE component is associated with the densest re 
gions, heated by the transmitted shock, whereas the hot 
ter, NEI component is gen erally attributed to the plasma 
heated by reflected shocks (iPark et al. [ bood l IZhekov et al 



12010.) . However, it is not clear whether also reverse shock 
heated plasma is responsible for some of the harder X-ray 
emission. 

Note that even the reflected/reverse shock heated plasma 
must be relatively dense as the ionization parameter indi- 
cates He ~ 950 Jjq' cm^^, with f lo the time since the plasma 
was shocked in units of 10 yr. As the SNR ages, the low 
temperature component is becoming hotter (kT^ w 0.2 keV 
in 1999 versus kT^^ ss 0.3 keV in 2005), whereas the 
high temperature component is becoming coo ler {kTg 
3.3 keV and kT^: » 2 . 3 keV , respe ctively), see IPark et al 
(l2006l) : IZhekov et al.l (I2009l l2010'). 



SN 1987 A is an ideal target for the grating spectrom- 
eters of Ch andra, due to it s small angular extent and its 
brightness ("Burrows et al.', '2000; 'Michael et al.l 120021: 
Dewey et al., 2008: Zhe kov et al ., 2009, Fig.[32]i, which 



allows for resolving the He-like triplets (Sect. l6T3T l. The 
earliest high resolution spec trum, taken with the Chandra 
HETGS, dates from 1999 dBurrows et all (2000). It re- 
veals prominent emission lines from H-like and He-like 
transitions of O, Ne, Mg, and Si, and line widths consis- 
tent with an expansion of ^ 4000 km ^ . Recent HETGS 
and LETGS spectra indicate a line broadening ranging 
from 1000-5000 km^^^ (FWHM), but with a bulk veloc- 
ity that is surprisingly low, ^ 100 — 2 00 km ' , consistent 



with the radia l veloc ity of the LMC dDewev et al.l 120081 



Zhekov et al.L l2009h . This is a strong indication that the 



hottest component is not due to the main blast wave, for 
which the plasma velocity should r nove with |Vi; , but rn ust 
result from either reflected shock fehek ov et al. , 2009h or 
from the reverse shock heated plasma. This would recon- 
cile the large line widths and hot temperatures, indicating 
large shock velocities, with a relative low velocity in the 
observers frame (see Sect.|4]and Sect. 15. lb . 

X-ray spectra taken with different instruments do not 
agree in all details. For example the Chandra CCD 
spectra indicate a cool component consistent with CIE 
dPark et al. , 2006 ). whereas the high spectral resolution 
HETGS/LETGS spectra indicate a plasma with a higher 
electron temperature, 0.55 keV, and out of ionization equi - 
librium, with n^t^ 4 x 10^' cm^^s dZhekov et al.1 l2009h . 
The best fit RGS model gives a solution that lies some- 
where in between nef~ 8 x 10^' cm^^s(Heng et al.,2008). 

In reality a range of temperatures and n^t values are ex- 
pected, given the complexity of the ring structure. Even 
the lowest n^t values reported, imply densities in ex- 
cess of He ~ 2000 cm^^. Nevertheless, these differences 
in fit parameters may also result in differences in abun- 
dance determinations. In general the abundance determi- 



et al.. I2OO8I: IZhekov et al.L l2009h . indicating LMC abun- 
dances except for oxygen, which seems underabundant, 
and nitrogen, which is mildy overabundant. The RGS and 
LETGS/HETGS mainly disagree on the iron abundance, 
which is underabu ndant according to fits to the RGS data 
dHeng etallEoOsl) . 

SN 1987A will remain an important object for X-ray 
astronomy in the coming decennia. It will brighten in 
the coming years as the transmitted shock will penetrate 
deeper and deeper into the ring. The main blast wave, 
which will pass around the ring, will at some point en- 
counter the shell that separates the tenuous wind of the 
progenitor from the dense red supergiant wind. This shell 
is likely to have a complicated shape, as indicated by the 
triple-ring system seen in the optical, and the hydrody nam- 
ical simulations of Morris and PodsiadlowskH ( 2007 ). The 
SNR will brighten as a result of encountering the shell, 
first just outside the equatorial region, later also in the po- 
lar regions. As the reverse shock penetrates deeper into 
the ejecta one will be able to probe the composition of the 
inner ejecta of SN 1987A. 

The evolution of SN 1987 A has already shown us that 
the early evolution of SNRs can be much more compli- 
cated than described by the standard evolutionary models 
described in Sect.|4] SN 1987 A may be a special case, with 
its triple-ring system and the possibility that its progeni- 
tor was the result of a recent common envelope merging 
of a binary system. However, such mergers are far from 
unique; even for the text-book shell-type SNR Cas A it 
has been argued that its progenitor was a merger product 
(Sect.l9j]l. 



9.3 



Core collapse supernova remnants and 
neutron stars 



Core collapse SNRs are apart from their abundance pat- 
terns, also distinguished by the presence of a neutron star. 
In the past it was assumed that this would automatically 
imply the presence of a pu lsar wind nebula (PWN, see 



Helfand and Becked 119841 for an early review on neu- 



tron stars in SNRs), but over the last 20 years it has be- 
come clear that not all neutron stars manifest themselves 
through the creation of a PWN dGaensler and Slane , I2OO6L 
for a review). 

There is good evidence that 5-10% of all neutron stars 
are born as magnetars, neutron stars with very high mag- 
netic fields (lO''* — 10'^ G), w hich emit pr edominantly in 
the soft and hard X-ray bands dMereghettil l2008t. for a re- 
view), and whose emission is not powered by pulsar spin- 
down, but by magnetic-field decay. These magnetars ro- 
tate slowly (P = 2 — 12 s), presumably due to a very rapid 
initial spin-down. As a result, their spin-down power is 
currently too low to drive the creation of a PWNr^l 



There is evidence that the most rapidly spinning magnetar 
(1E1547.0-5408 P = 2), located inside a faint SNR (G327.24-0.13 
[Oelfand and Gaensler, 2007) is_sumunde^^^_a^rrtPWN^Vink and 
Bamba, 2002), but this is disputed bv lOlausen et all i201 ll) . 
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Figur e 32: Combined Chandra HETGS count spectrum of SN 1987 A, based on a total exposure of 360 ks (iDewev et al. 
20081). The positive grating order is colored red, the negative order blue. (Credit: D. Dewey) 



Some neutron stars merely seem to manifest themselves 
as unresolved X-ray sources. By lack of a better name, 
they are usually called Compact Central Objects (CCOs, 
de Lucal l2008l for a review). By definition they are as- 
sociated with SNRs. It is not quite clear whether CCOs 
form a real class on their own. The lack of associated 
PWNe around them argues either for low surface mag- 
netic fields, hence small period derivates P, or long ro- 
tational periods, as the pulsar wind luminosity scales as 
E = Ak^IP/P^ (see Sect. |3]l. It seems that at least some 
of the CCOs are characteri zed by an unusual low mag - 
netic field ("anti-magnetars" lGotthelf and HalpernL 2008). 
Examples are PSR J 12 10-5226 in SNR PKS 1209-51/52 
{P = 0.42 s, Bp < 3.3 X 10^' G, iGotthelf and Halpern , 
2007h . PSR J1852-H0040 in SNR Kes 79 (P = 0.10 s. 



Bp < 1. 5 X 10^' G. Gotthelf et al., 2005; Gotthelf and 
Halpern, I2OO8I), and PSR RX J0822-4300 in the oxygen- 
rich SNR Puppis A rP = 0.11 s.Bp <9.8x lO'^ G. Got- 
thelf and Halpern, |200^ Apart from their low magnetic 
fields these stand out because of their relatively slow rota- 
tion period, which implies that they are born slowly rotat- 
ing. 

No pulsation period has yet been found for the point 
source in Cas A, but since it seems not surrounded by a 
PWN, it must either be a magnet ar or an anti-magneta r, 
and given its stable, but declining tHeinke and Ho , 2010t) . 
X-ray emission, it is likely to be an anti-magnetar. 

Another interesting CCO is 1E161348-5055 in RCW 
103. This source stand s out because it is highly variable 



dGotthelf et al.L Il999h and has an unusually large pulsa- 



tional period of 6.7 hr jGarmire et al. [ l2000l: iDe Luca et al 



2006h . This long period, either suggests that its is not the 
neutron star spin period, but instead the orbital period of 
the neutron star and an, as yet undetected, low luminosity 
companion, or the neutron star is magnetar-like and has 
been slowed down due to magnetic braking in the pres- 
ence of supernova fall-back material. 
So what can SNRs tells us about the birth properties 



and progenitor stars of the various types of neutron stars? 
In particular, a) what is the reason that neutron stars have 
a wide variety of properties, and b) how SNRs can help 
to infer the initial birth properties of neutron stars? These 
are questions for which observational data is scarce, as we 
know little about neutron star progenitors, and the birth 
of neutron stars is hidden to us by the supernova event it- 
self, with the exception of the detection of neutrinos from 
SN 1987 A, which marked the birth of the neutron star 
(Sect. |972] i. Strangely enough, no neutron star has yet been 
detected in SNi987A. 

From a theoretical point it is not clear what kinds of pro- 
genito rs produce di f ferent types of neutron stars or black 
holes dHeger et al.L l2003h , and wha t the in i tial spin pe- 
riods of n eutron stars ar e ( Ott et al. , 2006h . Heger et al, 
(l2003h and lNomotoeTaP (l2010h liave argued that less mas- 
sive progenitors (Mms ^ 25 M©) produce neutron stars, 
whereas more massive progenitors result in the creation 
of a black hole, either directly or due to l ate time super- 
nova material fall back. lHeger et alj (l2003h stress that this 
situation is modified by strong stellar winds, which reduce 
the massive hydrogen envelope. As a result very massive 
stars may create again neutron stars. This also depends on 
the initial metallicity of the star, as wind loss increases 
with metallicity. Indeed, there is evidence that magnetars 
are associated with very massive stars, based on the young 
ages of the stellar cluster in which thev are found ( Figer 



et al.. l2005l:lGaensler et alll2005l : lMuno et alll2006l). But 



the ev idence for this is not unambiguous (IDavies et al 
2009h and one should also keep in mind that even if some 
magnetars are associated with the most massive stars, this 
does not necessarily imply that all magnetars are. 

The timing properties of pulsars can be used to deter- 
mine their characteristic age. 



IP 

IP' 



(54) 



which is, however, only a reliable age estimator if the ini- 
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tial spin period Pq <^ P, and if the spectral braking index 
is n = 3 3 which is the value for braking due to a ro- 
tating magnetic dipole. The presence of neutron stars in 
SNRs helps to verify these assumptions as Tc can be com- 
pared with age estimates of the SNR. This helped establish 
that for the neut ron star in Puppis APn^i P ( Gotthelf and 
Halpern,i2009). 

The presence of a SNR surrounding a PWN can also be 
used to constrain the initial period in another way, namely 
by considering the pressure c onfinement of the PWN by 
the shock heated SNR shell (Ivan der Swaluw and WuL 
200ll:IChevalie3.l2005 ^. The total pressure inside the PWN 
is determined by the integrated energy loss of the pulsar, 
which in itself is mostly determined by the initial spin pe- 
riod. A surprising result of these studies is that there is a 
large range in initial spi n periods, which is theoretically 
not well understood fe.g. lOtt et al. [ l2006h . Moreover, and 
this brings us to the above formulated question a), there 
is no clear relation between the initial period of a neu- 
tron star and the type o f supernova as d etermined from 
the SNR characteristics /Chev alieiiboOSh FI 

The same conclusion can be drawn when also consid- 
ering the other classes of neutron stars, CCOs and mag- 
netars. First of all, the discovery of these neutron stars 
owe much to the fact that unresolved X-ray sources were 
discovered in SNRs. Most CCOs and magnetars that are 
more than a million year old are very difficult to discover 
The ceo will not be surrounded by a SNR anymore (and 
hence no longer deserve their name), and will have cooled 
substantially. Magnetic field decay makes that old magne- 
tars become radio-quiet dim X-ray sources. Neutron stars 
in SNRs are therefore important to estimate the the neu- 
tron star birthrate, especially of radio-quiet neutron stars, 
and to estimate the initial spin-period and magnetic-field 
distributions. 

As for the characteristics of the SNRs with CCOs and 
magnetars. At least two CCOs are found in oxygen-rich 
SNRs (Sect. l9Tt . the remnants of the most massive stars. 
For one of them (in Puppis A) we know that the neutron 
star has a low magnetic field. For the other one (Cas A), no 
pulsation has yet been found, and its nature is still unclear. 
For magnetars, there is so far only circumstantial evidence 
linking them to the core collapse of the most massive stars. 
There is also evidence that a perfectly normal neutron star 
can be created from a very massive star: The oxygen-rich 
SNR G292,P-Hl.8 , has a pulsar withf g^O 135 s (Hughes 
et al. T2003bh . and B w lO'"' G fchevaheri l2005h . So the 



conclusion must be that the most massive stars make neu- 
tron stars with any type of magnetic field (low, normal, or 
magnetar-like). 

This SNR/neutron star combination of G292.0H-1.8 
(Fig. |26] Sect. 19. It stands out for another reason as well. 
Its X-ray spectrum has been used to determine an abun- 
dance pattern that suggests that G292.0H-1.8 is the result 
of the explos ion of a star with an in itial mass of M = 
30 - 40 Mq jGonzalez and Safi-Harbl 120031) . This is the 



mass range for which theorists have argued that core col- 
lapse should lead to the creation of a black hole, instead 
of a neutron star. Clearly this is important for our under- 
standing of neutron star and black hole formation, but the 
evidence needs to be further studied. 

In contrast, no neutron star has yet been found in the 
Cygnus Loop (Fig. |4] see Sect. 110. lb . If this means that 
the event did not produce a neutron star but a black hole, 
we have an example of a less massive core collapse event 
(Mms ~ 15 Mq, Sect. 110.2b that did not produce a neutron 
star. But the Cygnus Loop is very big and perhaps a neu- 
tron star is inside or in the vicinity of the SNR and has not 
yet been identified. 

As for magnetars, the origin of their magnetic fields is 
not clear An interesting hypothesis is that the magnetic 
fields result from a d ynamo mechanism in rapidly ro tating 
proto-neutron stars (IDuncan and Thompson! 1 1992h . This 
requires rather short initial spin period of fo ~ 1 — 3 ms. 
Since the rotational energy of such a rapidly rotating neu- 
tron star is £^01 = 5^^^ w 3 x 10^^ (y^) ^ erg, one ex- 
pects that a lar ge part of the r otational energy will powe r 
the supernova dAronsl l2003l : lAllen and H orvathl |2004 . 
But an X-rav study of SNRs w ith magnetars bv Vink and 
Kuiper shows that the SNRs are consistent with 

typical explosion energies of 10^' erg. This is more consis- 
tent with the alternative hypothesis concerning the origin 
of magnetar fields, namely that magnetars are born from 
the core collapse of stars with very large magnetic fields 



(the f ossile field hvpothesis lFerrario and Wickramasinghe , 
2006). There is another, albeit non-conclusive, argument 
against the idea that magnetars are born rapidly rotating. 
As explained above, there is a large spread in initial spin 
periods of normal neutron stars. In fact there is very little 
evidence for pulsars born with periods f n < 10 ms ( van 
der Swaluw and Wu. l200ll:IChevalieill2005h . How likelv 
is it then that 5-10% of the neutron stars are born with 
extremely short periods, fo ~ 1 ^ 5 ms? 

The debate on the origin of magnetar magnetic fields 
has not yet settled, but it is cl ear that SNRs play an im- 
portant role in this debate (e.g.jPurant and van KerkwiikL 



20061: iHorvath and AUenL l20m 



2*" n is defi ned as a -kCl " 
^■'See also lChevaliel <2010h . 



10 Supernova remnants in or ap- 
proaching the radiative phase 

10.1 From non-radiative to radiative shocks 

A couple of thousand years after the explosion SNRs 
have swept up considerably more mass than the super- 
nova ejected. SNRs in this phase of their evolution are 
often referred to as "mature SNRs" (Sect. |4|i. Once their 
shocks have slowed down to below Vs < 200 km ' , the 
shocks are becoming radiative, and the SNR evolution can 
no longer be described by the Sedov-model (Eq.|2|. Slow- 
ing down even more results in a post-shock region that 
is too cool to emit substantially in X-rays. Nevertheless, 
some SNRs with slow shocks have bright interior X-ray 
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emission, as will be discussed in Sect. 110.3] 

The prototypical "mature" SNRs are the Cygnus Loop 
(Fig. |4^) and the Vela SNR, which are on the border be- 
tween phase II and phase III (Sect. |4|i. They are both con- 
sidered to be core collapse SNRs, as I will describe below. 
They have earned their prototypical status predominantly 
by being nearby examples of mature SNRs. Their small 
distance also means that their X-ray emission can be stud- 
ied on very small physical length scales. And because the 
absorption column is low, one can also study their soft X- 
ray and UV emission. 



The Cygnus Loop is at a distance of 540 ± 1 00 pc (IB lair , 
2005h . and has an angular radius of 1.4°, corresponding 
to 13 pc, with a break-out region in the south. Its age 
has been estimated to be 8,000 to 14,000 vr (Levenson 
et al., 1998, 2002; Katsuda et al., 2008c), and its shocks 
have velocities in the range of 150-400 kmjf~^ ("Lev enson 
et al.. ll998ll2002l : lBlaiii l2005': Salv esen et alll200 9^. This 
range in shock velocities indicates that the Cygnus Loop 
has both radiative and non-radiative shocks. In fact, the 
Cygnus Loop derives part of its fame from the beautiful 
optical nebula associated with the radiative shocks, the so- 
called Veil Nebula. The large contrasts in shock proper- 
ties are best explained by the idea th at the Cygnus Loop 
is evolving in a wind-blown cavity ( McCrav and Snow l 
1979l: ICharles et al.lll985l:lLevenson et al.ll 19981: Mivata 
and Tsunemi,'l999'). Earlier in its evolution the SNR must 
have expanded rapidly within the tenuous medium inside 
the wind-blown bubble. But currently most of the shock 
is interacting with the dense shell swept up by the progen- 
itor's wind (see also Sect. |4]l. In such a case parts of the 
shock may reach the radiative phase at a relatively young 
age. Another example is the proba bly even y ounger SNR 
RCW 86 (probably SN 185. Stephenson and Gr een. 2002), 
which also shows a mixture of radiative and non-radiative 



shocks (e.g. ISmithl 119971 : iGhavamian et al.L 1200 Ih . and 
whose plasma ionization ages and emission meas ures tes- 
tify of large densitv contrasts dVink et all 1 19971: Bamba 



et al..l2000':'Bo cchino et a l.'.'2000':'Rho et al. ll2002l: Vink 
et al., 2006). As a "cavity SNR" the Cygnus Loop may 
therefore not deserve its status as a prototypical mature 
SNR. On the other hand, all SNRs when investigated in 
detail appear to have some unique features. 

X-ray spectroscopy of the Cygnus Loop has provided 
some peculiar results. One example is that imaging spec- 
troscopy with ASCA, Suzaku, Chandra and XMM-Newton 
has revealed that the abundances of the bright X-ray shell 
are sub-solar, with typical depletion factors of 5 (Miy- 
ata et al, 1 994, 2007: Katsuda et al., 2008b: Nemes 
et al.7 E008t) l^ The low abundance of the inert element 
Ne shows th at this is not due to dust depletion. In contrast, 
Katsuda et al. (2008b) found that near the rim of the bright 
shell there are abundance enhancements of about a factor 
2. The same team later attributed these enhanced abun- 
dances at the very edge of the SNR to charge exchange 
reactions between unshocked neutral hydrogen and shock 



ionized O VIII (IKatsuda et al.L 1201 ih . However, as they 
point out, further investigation is clarify the possible role 
of resonant-line scattering on the line equivalent widths. 
As discussed in Sect. 16.1.31 this process results in scat- 
tering of photons with energies corresponding to resonant 
lines out of the line of sight. As the photons will eventually 
escape the SNR, it will result in enhanced line emission at 
the edges of emission regions, and an attenuation of reso- 
nant line emission from the regions projected toward the 
interior Indeed, there is some evidence for resonant line 



scattering (IMivata et al.Ll2008[) . but currently it is not clear 
whether charge exchange and/or resonant line scattering 
can explain the X-ray line emission from the rim of the 
Cygnus Loop, or whether the metallicity of the plasma is 
really depleted. Future high resolution X-ray spectroscopy 
may be needed to clarify this issue, as it is needed to re- 
solve the resonant lines from forbidden lines. Resonant 
line scattering will boost the G-ratio toward the center of 
the SNR, and lower it at the edge of the rim. Charge ex- 
change will only enhance the G-ratio at the edge of the rim, 
but without affecting the G-ratio in the rest of the shell. In 
Sect. 110.21 1 will discuss the X-ray spectroscopic results of 
inner regions of the Cygnus Loop. 

One striking feature of mature SNRs like the Cygnus 
Loop is the complexity of the emission, which is caused 
by the interactions of the shocks with density inhomo- 
geneities (clouds, or cloudlets), and by a variety of shock 
velocities giving rise to both radiative and non-radiative 
shocks. Moreover, these shock-cloud interactions lead to 
the creation of reflected and transmitted shocks (Levenson 
et al., 1996; Patnaude et al., 2002), not unlike the situation 
described in connection wifli SN 1987A (Sect.l9^. 

A SNR that has been used quite extensively to investi- 
gate the interaction of the blast with the inhomogeneous 
CSM is the Vela SNR. The Vela SNR contains a pulsar 
with a rotational perio d of 89 ms, which is surrounded by 



a pulsar wind nebula (IHelfand et al.Ll2001h . Its character- 
istic age is ~ 1 1 , 000 yr. The ac tual age of the pulsar and 
SNR may be older, - 20, 000 yr dAschenbach et allll995h . 
VLBI parallax measure ments of the pu l sar sh ow that its 
distance is 287 ± 19 pc (jPodson et al. , 2003 ). at which 
distance the angular size of the SNR, 4°, corresponds to a 
physical radius of 20 pc. The SNR is located in a complex, 
starforming region, which also contains the 7^-Velorum 
massive stellar binary. Both of them are embedded within 
the hot bubble/HII region that is known as the Gum-nebula 
dAschenbach et al.lll995l:rSushch et al.ll201 ll) . 

ROSAT and XMM-Newton X-ray imaging spectroscopy 
of isolated regions, in combination with optical narrow fil 
ter imaging show a complex struc t ure dBocchino et aT 



2OOOI: ILu and Aschenbachl I2OOOI: iMiceU et all l2005h 



Although this is not confirmed bv lLevenson et al.l )2002l) . who only 
analyzed a small portion of the SNR with Chandra. 



that is best described in terms of three components: 1) 
a tenuous inter-cloud medium, which is difficult to char- 
acterize in X-rays due to its diffuse nature and low sur- 
face brightness, 2) soft X-ray emission from cloud cores, 
heated by slow transmitted shocks, and with temperatures 
around fer w 0.1 keV, and 3) hotter regions surrounding 
the clouds, consisting perhaps of plasma evaporated from 
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Figure 33: False color map of the Vela SNR as observed in X-rays with ROSAT ( [\schenbach et al. , 1995 ). The figure 
points out the location of the "shrapnels" and their designations. (Source: Max-Planck-Institut fr extraterrestrische Physik 



(http : //www .mpe . mpg . de I.) 



the clouds (see also lKahn et alill985h : this plasma could 
be heated by thermal conduction from the hot intercloud 
medium, but may otherwise cool by radiation and ther- 
mal conduction with the cool cloud cores. The cloud cores 
have a moderately high densities of n w 5 cm^^. 

10.2 Enhanced metal abundances mature 
SNRs 

One striking result of imaging X-ray spectroscopy of ma- 
ture SNRs has been the identification of metal-rich plas- 
mas even within mature SNRs. In Sect. [8] and Sect. l9.1l a 
few mature SNRs were mentioned that could be identified 
as Type la (e.g. Dem L71) or core collapse/oxygen-rich 
SNRs (e.g. N132D and Puppis A) based on their abun- 
dances. However, even older SNRs such as Vela and the 
Cygnus Loop have regions with high metal abundances. 

Despite its large angular size, a significant fraction of 
the Cygnus Loop has now bee n mapped with Suza ku and 



XMM-Newton CCD-detectors jUchida et al.L 12009 ). This 



has revealed that metal-rich plasma is associated with the 
low surface brightness interior of the SNR, which has 
a temperature that is considerably hotter {kT k, 0.6 keV, 



Tsunemi et al.L 120071) than the bright shell {kT k, 0.2 keV). 



Even within this hotter region there is some variation in 
abundance patterns, with O, Ne, and Mg, more abundant 
in the outer regions, and Si, S, and Fe mo re abundant in 
the c entral region (Fig. 13 41 Tsune mi et all 2007: Uchida 
et al.. l2008h . Recently also the detection of Ar lines was 
reported, indi c ating high Ar abundances (8-9 times solar. 



Uchida et all 1201 Ih . This suggests a layered explosion. 



with the more massive elements situated in the center This 
is in contrast to Cas A, for which the Fe-rich ejecta has 



(partially) overtaken the Si-rich ejecta (Sect. lgTI ). 

The overall abundance pattern of the inside of the 
Cygnus Loop is consistent with that expected from a 
core collapse supernova with an initial progenitor mass of 
M ss 15 M©, but with 5-10 tim es more Fe th an predicted 
by supernova explosion models (ITsunemi et al., 2007). As 
noted in Sect.|2]the predicted Fe yield has considerable un- 
certainty. The idea that the Cygnus Loop progenitor was 
not a very massive star (M > 20 Mq) is also consistent 
with the radius of the stellar wind bubble of 1 3 pc, which 
is more in line with an early type B-star than with an early 
type 0-stai- (Chevaher, 1990). 

It is tempting to speculate that the more layered metallic- 
ity structure of the Cygnus Loop, as opposed to the more 
radially mixed metallicity of Cas A, has to do with the 
different progenitor types. The Cygnus Loop progenitor, 
with its relatively low initial mass, probably exploded as a 
red supergiant, but may still have had a significant fraction 
of the hydrogen envelope at the time of explosion. The Cas 
A progenitor was almost completely stripped of its hydro- 
gen envelope at the time of explosion (Sect. 19.11 ). which 
made it perhaps easier for convective motions inside the 
supernova explosion to affect the layer ing of the ejecta. In- 
deed the simulations of lKifonidis et al.l (12003) show a qual- 
itative difference in core ejecta velocities between Type II 
and Type lb (i.e. stripped progenitor) supernovae. 

The metal-rich plasma in the Vela SNR is not associated 
with the interior, but rather with a number of protruding 
plasma clouds, usua lly referred to as the Vela "bullets" or 



"shrap nels" (Fig. l33llAschenbach et al.Lll995l : IStrom et al 



1995h . The name shrapnels is justified by their high abun- 



dances (iTsunemi et al. . Il999h . which indeed suggest that 
the clouds consist of supernova ejecta. Although it is not 
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Figure 34: The emission pattern of oxygen (left) and iron line emission (right) from the Cygnus Loop as a function 
of position wi t hin th e SNR, following a track from roughly nort heast to southwest. T he green sym bols are taken from 
Tsunemi et alJ OOOTh and black and red points were determined bv lUchida et alJ (l2008h . (Figure from lUchida et all 120081) 



clear whether this also means that the clouds were ejected 
as such by the supernova, or whether they were formed 
due to some hydrodynamical instabilities at the contact 
discontinuity in the early SNR phase. The morphology of 
the SNRs are suggestive of clouds surrounded by bow- 
shocks, whose opening angles indicate Mach numbers 
M « 3 — 4. For the h ot, tenuous enviro nment of the Vela 
SNR (n 0.01 cm'^ IStrom et alllT995i) . this coiTesponds 
to velocities of V « 300 - 700 km I The average shrap- 
nel densities are ~ 0.1 — 1 cm^^ with masses for shrap- 
nel A and D of, r espectively. 0.005 MrT^and 0.1 Mf7i( Kat- 
suda and Tsunemi l2005Ll2006h . 



Given their projected distance from the main shell of 
~ 1°, corresponding to / w 5 pc, the shrapnels must have 
survived passage through the hot bubble for quite some 
time. Even for a mean average velocity difference between 
shock and shrapnels of Av w 1000 km 5^', the travel time 
is ftrav ^ //Av ^ 5000 yr This sho uld be co mpared to the 
typical cloud-crushing time fee (.Klein et al 



occurring in shrapnel D dKatsuda and Tsunemil l2005h . 
which wiU lead to the destruction of this cloud. 

As already in dicated, imaging spectroscopy with A SCA 
dTsunemi et alld999). Chandra jMiyata et al.[.l200lh . and 



XMM-Newton dKatsuda and Tsunemill2005L l2006h have 
shown that the clouds represent supernova ejecta. These 
studies also indicate that the shrapnels originate from dif- 
ferent layers inside the supernova, with Shrapnel A being 
more rich in the oxygen-burning product Si, and shrapnel 
D more rich in Ne and the carbon-burning products O, Ne, 
Mg (with enhanced abundances by a factor 5-10 with re- 
spect to solar). The shrapnels are apparent due to their pro- 
jected distance from the Vela SNR. However, more shrap- 
nels may be project ed unto the main SNR shell. Indeed, 
Mice li et al. ( 2008 ) found various regions observed by 
XMM-Newton that have overabundances of O, Ne, Mg and 
Fe, which may be either indicative of other shrapnels, or 
which may be regions containing the remains of destroyed 
shrapnels. 



X 



Vb 



100/ VlOOOkms-' 



-1 



R 

1 pc 



= 980|t^| I—-— rl 1^ I yr, 

(55) 

with X the density contrast between the shrapnel and its 
ambient medium. Comparison of these time scales, assum- 
ing Vb ~ Av, shows that the apparent survival of the shrap- 
nels indicates that the density contrast between the shrap- 
nels and the ambient medium must be high % ~ 1000. A 
higher average velocity does not help, as it lowers both 
the travel time and the destruction time. A higher density 
contrast could either mean that the initial density of the 
knots was higher than it is now, or that the ambient den- 
sity is lower, ^ 0.001 cm^^. Either way, the Vela shrapnels 
may owe their survival to the low density of the ambient 
medium and a fortuitous timing concerning the phase in 
which we happen to observe the Vela SNR. There is some 
morphological evidence for hydrodynamical instabilities 



10.3 Mixed-morphology supernova rem- 
nants 

Among the mature SNRs, one class of objects sticks 
out: the mixed-morphology (or thermal-composite) SNRs 
(Sect. O. They are characterized by centrally dominated, 
thermal X-ray emission, whereas their radio morphology 
is shell-like. But apart from the curious radio/X-ray mor- 
phology they have a number of other characteristics in 
common (e.g. ICox et all Il999l) . They tend to be older 
SNRs (> 20,000 yr) and are associated with the denser 
parts of the ISM. Many of them are associated with OH 
masers, which is an indication of interacti ons between the 
SNR shocks and molecular cl ouds (e.g. Iciaussen et al 



1997l:lFrail and Mitchellll 19981 Table gll. Several of them 
are also sources of GeV and even TeV 7-ray emission (Ta- 
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ble|4]i, which likely arises from cosmic ray nMcZeO inter- 
acting with dense gas in the SNR shell, or from escaping 
cosmic rays interacting with unshocked parts of molecu- 
lar clouds. These interactions produce pions, of which the 
neutral pions decay into 7-ray photons. Finally, imaging 
X-ray spectroscopy of mijxed-morphology SNRs shows 
that their X-ray emitting plasmas have more or less homo- 
geneous temperatures. 

The morphology of these SNRs is difficult to explain 
with standard SNR evolution models (Sect.|4]l. There have 
been severals models to explain both the centrally en- 
hanced X-ray emission, which contrasts with their shell- 
type radio appearances (Fig.l4H). The model of White and 
Long (1199 11) assumes that the central thermal emission is 
caused by dense cloudlets that have survived the passage 
of the forward shock, but are slowly evaporating inside 
the hot medium due to saturated thermal conduction. They 
found a self-similar analytical model for the evolution of 
thermal-composite SNRs, which in addition to the Sedov 
solutions (Sect. HI has two additional parameters, namely 
C, the ratio between the mass in the clouds and the mass in 
the intercloud medium, and T, the ratio between the evap- 
oration time of the clouds and the S NR age. 



1999')( 



An alternative model proposed bv lCox et al. 
centrates on the contrast between the radio and X-ray mor- 
phology. According to their model the SNRs evolve inside 
medium density CSM. The forward shock has decelerated 
to velocities below ^ 200 km.s^', resulting in a strongly 
cooling shell that is too cool to emit X-ray emission. Be- 
cause the shocks are radiative, the shock compression fac- 
tors are large, giving rise to strongly compressed magnetic 
fields, and enhanced cosmic-ray electron densities. As a 
result the radio synchrotron emission from the shells is 
strongly enhanced. The electrons are either accelerated 
by the SNR, or are simply ambient relativistic electrons 
swept up by the radiative shock. The X-ray emission from 
the interior has a homogeneous tem perature due to ther- 
mal conduction ( Cui and CoxL 1 1 992h and turbulent mixing 
( Shelton et al. . 2004a'). Numerical hydrodynamical mod- 
els based on the model of Cox; et al. ( 1999) indeed show 



centrally enhanced X-ray m orphologies fshelton et al 
19991: 1 Velazquez et al.l |2004'). 



Although the nature of mixed-morphology SNRs is 
still debated, there has been considerable p rogress in our 
knowledge of them over the last ten years. ^ Shelton et al. 
(|2004b) showed that SNR G65.2H-5.7, which is clearly in 
the radiative evolutionary stage 0, shows an X-ray filled 



^'As opposed to electrons. See also the discussion in Sect. 1 11.21 
^"Xhis explanation for synchrotron emission from SNRs is called the 
Van der Laan mechanism (van der Laan, 1962). Uchiyama and Aharo- 
nian (2008) recently proposed that the radio emission may be addition- 
ally enhanced by the presence of secondary electrons/positrons, i.e. the 
electrons/positrons products left over from the decay of charged pions, 
created due to cosmic ray nuclei colliding with the background plasma. 
The presence of secondary electrons/positrons may also exp l ain the flat 
spectral radio indices Uj^ in Table |4] Finally, iBvkov et alj )2000l) dis- 
cusses the possibility of electron acceleration by magnetohydrodynamic 
turbulence behind high density shocks. 

^*See for example the well defined shell in optical [OIII] emission 
iBoumis et aill2004i) 
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Figure 35: Suzaku SIS spectrum of G359.1-05 ( Ohnishi 
et al., 201 1), showing in green the Si and S Lya line and 
in blue the radiative-recombination continua. 



morphology. Due to the low absorption column toward 
this SNR, ROSAT was also able to reveal a shell-like outer 
part in X-ray emission below 0.3 keV. These two charac- 
teristics are in qualitative agreemen t with the explan ation 



1999h . 



for mixed-morphology remnants bv lCox et al 

Due to progress in X-ray imaging spectroscopy over 
the last 15 year, several other clues about the nature of 
mixed-morphology SNRs have appeared. Many of the 
mixed-morphology SN Rs show evidence for metal -rich 
plasmas in the interior: Lazendic and Slane ( 2006h list 
10 out of a sample of 23 thermal-composite SNRs with 
enhanc ed abundances. Examples are W44 Is helton et all 
2004al). HB 21. CTB 1 and HB 3 jLazendic an d Slaii^ 
20061: IPannuti et al.ll2010bh. IC 433 dTroia et al. . 2008.) 



and Kes 27 (iChen et al.L 120081) . Kawasaki et"aiT (120051) 



noted another feature that may be generic for this class of 
SNRs: their ionization state is close to ionization equilib- 
rium (Sect. 16.1.21 ). unUke young SNRs and old SNRs like 
the Cyg nus Loop. Moreover, the thermal-composite SNR s 
IC 443 dKawasaki et all I 2OO2I: lYama guchi eTall l2009h . 
W49 B JOzawa et all 120 09^1 and G3 59.1-0.5 ( Ohnishi 
et al.,'2011') show signs of overionization, in the form of 
radiative-recombination continua (RRCs, Sect. I6.1.2I I as- 
sociated with metals Uke Si, S, and Fe (Fig.l35ll. 

It is not clear whether metal-richness is a generic feature 
of mixed-morphology SNRs, but it seems at least to be a 
common feature. There are several ways in which metal- 
richness helps to explain the characteristics of thermal- 
composite SNRs. Fkstly, the uniform temperatures in their 
interiors are ofte n attributed to thermal conduction (e.g. 
Cox et al. . Il999h . A problem for this explanation is that 
magnetic fields Umit thermal conduction across field lines 
(Sect.|53|i. However, if the interiors consist predominantly 



W49B is a young ejecta dominated SNR, and it is not clear whether 
it belongs to the class of thermal-composite SNRs. However, apart from 
its age it has thi'ee characteristics of thermal-comp osite SNRs: a high am - 
bient density, centr ally enhanced X-ray emission jRho and Petrdfl998l) . 
and 7-ray emission iAbdo et al.Lr2010j) . 



50 



1 0.3 Mixed-morphology supernova remnants 



J.Vink Supernova remants 



Table 4: Mixed morphology SNRs. 
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1.9 
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X(n) 


G3 1.9+0.0 
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24 


0.49 


8.5 






X(k) 




X(n) 
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This list is based on the list by lLazendic and Slanel |2Q06|) . Sizes, distances, and radio properties (radio flux Fr, and spectral index a^) were obtained from Green D. 
A., 2009, "A Catalogue of Galactic Supernova Remnants (2009 March version)". Astrophysics Group, Cavendish Laboratory, Cambridge, United Kingdom (availa ble at 
[http://www.nirao.cam.ac.uk/surveys/snrs/ The column names indicate: Enh. metals, the presence of plasma with overbaundances iLazendic and Slane. 20061 
unless indicated otherwise); RRC, the presence of radiative recombination line continua; GeV, the detection in the GeV band (Fermi); TeV the detection of associated sources 

in the TeV band with Cherenkov telescopes; OH Maser, the detection of associated OH masers. The letters in parentheses refer to the following references^ 

a)[Abdo et al. i2009); b) Abdo et al. 1 2010a); c) Abdo et al. 1 2010b); d) Abdo et al. i2010d); e) Abdo et al. i2010e); f) Abdo et al. 1 20101); g) Aharonia neTan |2008tft ; h) 
Albert et al. i 2007a); i) Bocchino et al. i2009); j) Brun et al. ( 2011); k) Castro and Slane ( 2010); I) Chen et al. i2008) m) Claussen et al. 1 1997); n) Frail et alj 119961) : o) 
Green et al. 1 1997) p) Koraleskv et al. 1 1998a); q) Lazendic et al. i 2002); r) Ohnishi et al. i 201 1); s) Ozawa et al. 1 2009); t) Park et al. i 2009); u) Yusef-Zadeh et alj |1995l) ; v) 
Yusef-Zadeh et aU i 1 996i) ; w) iYusef-Zadeh et aljU999i) ; x) iYamaguchi et alit2009i) ; v) >Feinstein et alJ<2Q09ft ; lKrause et alj<2Qlli) . 
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of supernova ejecta, then the magnetic fields may be very 
low, simply because the stellar magnetic field has been di- 
luted by the expansion. Magnetic flux conservation gives 
B = B^,{R^,/R)^, with/?* the stellar radius, and the radius 
of the SNR interior Taking typical values of lO'^ cm, 
R ss 10^" cm, Bh. = 1 G, shows that the interior magnetic 
field may be as low as lO^'"* G, approaching the value 
where conduction across field lines becomes comparable 
to conduction along field lines. Hence, the conditions are 
optimal for thermal conduction, but only up to the bound- 
ary between ejecta dominated and swept-up matter, where 
Rayleigh-Taylor instabilities at the contact discontinuity 
may have enhanced the magnetic-field strengths and may 
have made the magnetic fields more tangled, preventing 
strong thermal conduction across the contact discontinu- 
ity. 

Secondly, metal-rich plasma emits more X-ray emis- 
sion, further increasing the contrast in X-ray emission be- 
tween metal-rich interior and the cool shell of swept-up 
matter. An important additional aspect to consider here 
is whether the interiors consist of enhanced metal abun- 
dances, i.e. a mix of ejecta and swept-up matter, or of pure 
ejecta. For SNR W49B, a rather young member of the 
mixed-morphology class, a pure metal abundance seems 
likely, but for older SNRs the plasma is probably best de- 
scribed by hydrogen/helium dominated plasma with en- 
hanced abundances. One can illustrate the key features of 
the cool-shell/hot-interior scenario for mixed-morphology 
SNRs with some simple numerical calculations, ignoring 
for the sake of argument the more rigorous analytical treat- 
ment of the SNR evolution by Cox et al. ( 1999). 

Firstly, most mixed-morphology SNRs appear to be 
in the snowplough phase of their evolution, with Vs < 
200 kmi^', which results in cool. X-ray dim, but opti- 
cal/UV bright regions immediately behind the shock front. 
The interior, presumably consisting of supernova ejecta, 
is more or less homogeneous in temperature and, because 
of pressure equilibrium also reasonably homogeneous in 
density. Given the age of the SNR and the relatively high 
sound speed in the interior, one can assume approximate 
equilibrium betwee n the ram pressure at the forward shock 
and the interior (c.f. lCui and Coxll 19921) . This gives: 



'intenor'^^ mtenor 



(56) 



off, given the crudeness of the approximation. At least it 
suggests that a substantial fraction of the interior plasma 
consists of ejecta, and that the interior magnetic field may 
be low due to large expansion of the ejecta material, or 
early shocked swept-up material. 

Note that if we assume that the interior consists of met- 
als only, this does not lead to a consistent result. For exam- 
ple, if we assume O VII to be the dominant ion, then we 
have «ovii — n/% k 0.29 cm^^ and Minterfor ~ 100 Mq. 
This is clearly much more oxygen than a massive star 
can produce. So pure metal abundances seem unlikely, al- 
though perhaps there are some regions in the interior that 
have unmixed eject a. For a discu s sion on the metal abun- 
dances in W44, see lShelton et all ( l2004ah . 

A SNR will reach the snowplough phase when its age 
is larger than f^ad ~ 44,6000(£5i/nH)^''^ yr. For W44 
(jiH = 6 cm^^), this corresponds to 25,000yr. Hence, an 
ionization age of 4 x lO'^ cm^^s is expected (Sect. l6TT2T l. 
This is relatively high, but below the ionization age for 
ionization equilibrium (lO'^ cm^^s). However, the situa- 
tion is more complex, as the interior density and tempera- 
ture may have been higher in the past, and, therefore, the 
degree of ionization may be higher than indicated by the 
simple multiplication of the present day age and electron 
densit v. This could explain the observation bv Kawasaki 
et al. ( |2005|) that most mixed-morpholopv SNRs are near 
ionization equilibrium. Note that the high density in the 
interior of W44 and other mixed-morphology SNRs is a 
direct consequence of the high ISM density (Eq. l56T l and 
the thermal conduction, which results in a more uniform 
interior density with medium hot temperatures, rather than 
very high temperatures with very low interior densities. 

In the context of mixed-morphology SNRs, it is of in- 
terest to show why s ome old SNRs are ex pected to have 
RRCs. As noted by Ohnishi et al.l (2011), one may ex- 
pect RRCs whenever the plasma cools faster than the re- 
combination time, for e xample due to thermal conduction 
dKawasaki et al.1 2002) or adiabatic expansion (Itoh and 



Masai. Tl989T l^ima^uch^^^!n2009l) nYamaguchi et al.l 
(2009^ calculated that the time scale for cooling by con- 
duction is too long, so adiabatic expansion is probably the 
appropriate explanation. The time scale for adiabatic ex- 
pansion can be calculated from the thermodynamic rela- 
tion TV^^^ = constant, which means that the adiabatic 



with 



the total interior plas ma density . If we coohng time scale is given by (|Broersen et al.|,|20U: 



use the values for W44 listed by iCox et al 
Vs = 150 kms~\ ^Tinterior = 0.6 keV and po 
10^^^ gcm^^, we find a density of njnteiior 
which translates into an electron density of n 



(Il999h 
1.0 X 

2.3 cm-3 
0.8cm-3 



in rou gh agreement with X-ray measurements dCox et al 
[1999'). W44 has a radius of R ^ 10 pc, for an estimated 
distance of 2.5 kpc, but the X-ray emitting interior has a 
radius of /^interior ~ 6 pc, corresponding to an interior vol- 
ume of Vinterior ~ 2.7 X 10^^ cm^^. This means that the 
total internal energy is U — nkTVintehov ~ 0.9 x 10^" erg 
(about 10% of the explosion energy). The X-ray emitting 
mass is Minterfor ~ 28 Mp;,; somewhat higher than the ex- 
pected ejecta mass of a Type II supernova, but not too far 



Tad = - 



T\ _ 1 V 



1 R 



3(1 -r)R 



2R 2^ 
(57) 



with j3 the expansion parameter (Sect. l4!2T i. 7=5/3, and 
f the age of the SNR. For SNRs in the snowplough phase 
we can use j3 = 0.25 (Eq|8]i. This suggests that Tad ~ r/8. 
For overionization to occur it is first required that the ion- 
ization age of the plasma must be large enough (otherwise 
there could underionization). The second condition is that 
the timescale for recombination is longer than the adia- 
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batic cooling time scale: 



11 Shock heating and particle accel- 
eration: observations 



^ad Tree 



1 



(58) 



with ttrec the recombination rate. The mixed-morphology 
SNRs show evidence for overionization of Si XIII, for 
temperatures around kT = 0.2 — 0.7 keV. For these tem- 
peratures the recombination rate for Si XIV to Si XIII 
is approximately Ctr pr = 5.9 x 10~^^ cm^s~' ( Shull and 
van Steenberg, 198^!" Using now that «e — O.Sninterior 
and Eq. |56] the condition for overionization in mixed- 
morphology SNR is: 



/ < 



8 8/trin 



OiecWe O.SOrecPoV/ 



(59) 



■-1 irv6 — 1/ ''s \ / ^ ^interior \ 

2 X 10 n„ r 

" V150kms-iy VO.SkeV/ 



6x 10-14 cm^s-i 



yr, 



with «H the pre-shock density. This equation shows that 
overionization should be common, as most SNRs are 
much younger than a few million year. More important is 
whether the average plasma ever reached ionization equi- 
librium, i.e. whether / «eiir > lO'^ cm^^s. The discovery 
of RRCs in mixed-morphology SNRs has been relatively 
recent. Given Eq. |59] one expects that many more SNRs, 
with relatively dense interior plasmas, should have RRCs 
in their X-ray spectra. 

Finally, some word of caution considering the interpre- 
tation of mixed-morphology SNRs. The explanation of 
their properties sketched above is relatively simple, but 
nature is often not that simple. For example, one may ask 
how it comes that SNRs with maser emission, and which 
are interacting with molecular clouds, have densities of 
only MH = 5 cm^^. This is relatively high, but not nearly 
as high as the dense molecular clouds usually associated 
with m aser emission, nw > 1000 cm^^ jFrail and Mitchell , 
1998). But in fact, millimeter- wave and near infrared ob- 



servations of W44 and W28 bv lReach et al.l (l2005h reveal 
the presence of dense molecular gas. However, this dense 
gas has only a small filling factor. iReach et al.l ( l2005h even 
note that some of the dense clumps may survive the shock 
passage, which brings u s back to the scenario proposed by 
White and Long (Il99lh . The interclump densities, which 
applies to 90% o f the volurn e are, h owever, consistent with 



the estimates of Cox et al. (Il999h 



5 cm ^. These 



infrared observations show that old SNRs have a com- 
plex structure, which, for a proper understanding, requires 
multi-wavelength observations. 



Electron-ion temperature equilibra- 
tion 



11.1 



As explained in Sect. 15.31 from a theoretical point of view 
collisionless shocks do not naturally lead to equilibration 
of electron and ion temperatures in the immediate post- 
shock region. The degree of post-shock temperature equi- 
libration by SNR shocks is therefore an important ob- 
servational issue. Unfortunately, X-ray measurements of 
plasma temperatures usually refer only to the electron tem- 
perature, as it is the electron temperature that determines 
X-ray line ratios (Sect. l6T3] ), and the shape of the thermal 
X-ray continuum (Sect. lOI ). 

A first indication that the electron temperature is not 
equal to proton temperature, or at least that Eq. |20] 
(Sect. 15.1b does not hold for young SNRs, is the fact that 
for young SNR no X-ray temperatures above kT^. « 5 keV 
have been reportedly This is in contrast to the expected 
temperatures for some young SNRs with high shock ve- 
locities. For example, Cas A and Tycho have shock veloc 
ities of, respectively, 



5100 kms- 



19981; 



. '(|Vinketal._ 

iKoralesky et al.', 'l998bl; 'Pelanev and Rudni ckL l2003l: Pat 
naude and Fesen, 2009) and K w 4000 km s ^^ ( Revnoso 
et al.. ll997l : IKatsuda et aUl2010at) . They should, accord- 
ing to Eq.|20l have plasma components with kT^. > 19 keV. 
Even though part of the X-ray continuum of those two 
SNRs is due to X-ray synchrotron emission (Sect. II 1.21) . 
plasmas as hot as 20 keV should easily be detectable as 
they give rise to relatively flat spectra with power-law in- 
dices around 1.5. In contrast, the observed spectral indices 
in the 4- 6 keV continuum d ominated band are typically 
w 3 (e.g. iHelder and \^ l2008h . Also the hard X-ray 
spectra do not show signs of thermal X-ray emission frorn 
plasmas with kT > 5 keV ( Allen et al. , 1997 :lFavata et al 



ll997yAllen et al.Lll999l : IVink et allbOOlU Renaud et^ 



2006b). 

Apart from non-equilibration of electron-ion tempera- 
tures there may be another reason for the low electron 
temperatures in young SNRs: non-linear cosmic -ray accel- 
eration may be so efficient that the overall plasma temper- 
ature sarever^Jow_(SectJ^^ Drury et aU, 20091: Helder 



et al.. l2009l : IVinketalll2010[) .' 



In order to estimate the proton or ion temperature, as 
opposed to the electron temperature, one needs to mea- 
sure the thermal Doppler broadening. Thermal Doppler 
broadening gives rise to a Gaussian line profile with line 
widtll^ 



(60) 



One of the highest temperatures reported is kT,^ = 4.3 keV for the 
SNR IC 131 in M33 tTullmann et al. . 20091. 

The lack of a factor 2 in the equation may be confusing. Note that 
iRvbicki and Lightmaj il979h show that the profile is proportional to 
exp[— miC^(£ — £'o)^/(2£o^^i)]- Comparing this to a gaussian exp[— (a: — 
xo)^/{2cj^)] shows the correctness of Eq. l60l 
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The line width is often given in terms of Doppler veloci- 
ties, in which case (7v = ^kTilm\. In order to measure the 
thermal Doppler broadening one best uses a spectrum of 
a region that is expected to be devoid of bulk line of sight 
motions. This is expected at the edges of SNRs, where 
the expansion of the shell only results in velocities in the 
plane of the sky. 

Most measurements of thermal Doppler broadening are 
based on optical and UV observations of Balmer and Ly- 
man lines of hydrogen. A comprehensive l i st of th ese mea- 
surements is given by iGhavamian et al.l (l2007h . Balmer 
and Lyman emission lines from SNR shocks require that 
the shocks move through (partially) neutral gas. As the 
neutral hydrogen penetrates the shock heated region the 
neutral atom may get excited, giving rise to narrow line 
emission in the optical/UV. In addition to excitations also 
charge transfer reactions may occur, i.e. a shock heated 
proton may pick up an electron from a neutral particle. 
This process gives rise to t hermally broadened Balmer or 
Lyma n Une emission ("e.g. IChevalier et al. , ll980l: iHend. 



region very close to the shock front, and the current grat- 
ing spectrometers do not easily al low to isola t e part icular 



I2OIO') from which the hydrogen tem perature immediately 
behind the shock can be measured (Ivan Adelsberg et al 



2008). The ratio between the narrow and broad line emis- 



sion is a function of shock velocity and the degree of 
electron-proton temp erature equilibration. 

Ghavamian et al. (l2007h presented evidence that for 
low shock velocities the proton and electron tempera- 
tures are equilibrated, whereas for shock speeds above 
~ 500 km 5^ ' the temperatures become increasingly non- 
equilibrated, with a rough dep endency of T^ jT^ °^ 1 /V^. 



The w ork bv lvan Adelsberg et a l. (2008) and lHelder et al 



(1201 Ih does not corroborate the dependency, but does con- 
firm that slow shocks produce thermally equilibrated plas- 
mas, w hereas fast shocks d o, in general, not. Note that the 
work of Helder et"al] (|20 lib, and earlier work of Rakowski 
et al. ('2003'), is based on electron temperatures measured 
using X-ray spectroscopy, whereas proton t emperatures 



were measured from broad Ha line emission. Helder et al 



(1201 Ih used the SNR RCW 86, which has large density 



and s hock velocity contrasts ( Vink, 2006; Hel der et al 



2009h. and may be e fficient in accelerating cosmic rays 



dHelder et af, '2009'; 'Aharoni an erall l2009b. Rakowski 
et al. (2003) measured temperature equilibration in the 
LMC SNR DEM L71 with infeiTed shock velocities of 
500-1000 kms-K For DEM L71 the T^/Tp oc l/y2 rela- 
tion seems to be valid. 

In addition to the hydrogen lines, also lines from other 
ions can be used to measure ion temperatures. A note- 
worthy example is the UV line emission from the north- 
western region of SN 1006 as measured by the Hopkinson 
Ultraviolet Telescope (HUT), which shows that the lines 
from H I, He II, C IV, N V, and O VI all have the same 
width of - 2300 km5"'(FWHM iRavmond et"an. Il995h . 



This indicates that the ions are not equilibrated, since there 
is no dependency of the line width on the mass of the ion 
(see Eq.l60b. 

Thermal Doppler broadening measurements in X-rays 
are rare, as one needs high resolution spectroscopy of a 



regions of an extended object. But Vink et al. ( 2003) mea 



sured the line broadening of the O VII Hea line emission 
from a bright knot that stands out at the northeastern rim 
of SN 1006. The size of the knot is approximately 1'. This 
measurement was possible for SN1006, because this rem- 
nant is rather large for a young SNR (30'). As a result 
only a small amount of the overall X-ray emission from 
the rest of the SNR is contaminating the XMM-Newton 
RGS spectrum. The expected line profile for a narrow line 
consists of a sharp rise, as for a normal point source, but 
it has a long tail, caused by emission from the interior 
of SN 1006 (Fig. 136b. The shape of the short wavelength 
side of the line profil e is, therefo re, sensitive to the ther- 
mal line broadening. I Vink et al.l (2003) reported a line 



profile with a width of <7 = 3.4 ±0.5 eV, corresponding 
to (7,, = 1777 km.?^'. Using Eq.|60] this translates into a 
temperature for O VII of kToya = 528 ± 150 keV. This 
seems extremely hot, but using Eq.|20l with nii = Ibm^ in- 
stead of /imp, shows that it corresponds to a shock veloc- 
ity of ^ 4000 km ^ , provided that oxygen has not been 
equilibrated with protons. The high O VII temperature 
should be compared to the electron temperature as mea- 
sured with the EPIC instruments on board XMM-Newton, 
kTf^ w 1 .5 keV. This is a clear indication that in this elec- 
trons and oxygen ions are not equilibrated. 

The nature of the X-ray emitting knot in SN 1006 is 
not quite clear It seems to have enhanced O and Si abun- 
dances, and may therefore be an ejecta knot, in which case 
it was heated by the reverse shock. Another important is- 
sue is whether non-thermal components, such as acceler- 
ated particles may have influenced the thermal balance of 
the shock, as discussed in Sect.. 15. ll The northeastern part 
of SN 1006 seems, in that sense, a relatively good place 
to measure equilibration processes in the absence of cos- 
mic rays, as this region shows no evidence for X-ray syn- 
chrotron emission (see also Sect. 111. 2] |. 

Although the RGS observations of SN 1006 are perhaps 
the most direct measurements of thermal Doppler broad- 
ening in X-rays, thermal Doppler broadening is important 
to take into account when using the overall line broaden- 
ing of SNRs. An example is th e Doppler line broade ning 
of LMC S NRs like 0509-67 5 dKosenko et all[2008h and 
0519-69.0 dKosenko et all I2OI0I) as measured with the 
XMM-Newton RGS. Thermal line broadening may even 
affect the line w i dth as measured with CCD instruments. 
Furuzawa et al. 1(12009') reported that Fe-K line broaden- 



ing of X-ray line emission from Tycho's SNR, and its 
spatial distribution, cannot be easily modeled with bulk 
velocities alone, but that thermal Doppler broadening is 
needed as well. From this they infer an iron temperature of 
(1 - 3) 10^° K (860-2600 keV)! So there are several pieces 
of evidence that electrons and ions may indeed have differ- 
ent temperatures. In X-rays, however, the evidence seems 
at the moment to pertain to plasma heated by the reverse 
shock. The evidence for non-equilibration of temperatures 
of plasma heated by the forward shock is mostly based on 
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Wavelength (A) 

Figure 36: Left: Map of OVII emission made from several Chandra- ACIS observations (see also Fig. 137b. The lines 
indicate the region observed by the XMM-Newton RGS instrument. The target was the bright knot in the northeast inside 
the RGS region. Right: Detail of the RGSl spectrum of the northeastern knot, showing O Vll Hea line emission. The 
dotted line is the best fit model without lin e broadening, whereas the solid line shows the model including thermal line 
broadening (adapted from lVink et all 120031) . 



Ha measurements. For these measurements the amount of 
equilibration is based on the narrow to broad line flux ra- 
tio, but the interpretation is m odel dependent and does no t 
always give consistent results (Ivan Adelsber^ et alil2008h . 
Moreover, current models do not take into account the ef- 
fects of cosmic -ray precursor physics, which is likely im- 
portant for the correct interpretation of the line ratios (van 
Adelsberg et al., 20 08; Rakowski et al., 2008: Held er et al., 
201 itlRavmond et al.ll201 ll) . 

It could be that the dichotomy between equlibrated 
shocks below 500 km 5^', and unequilibrated shocks 
above 500 kms^^ is somehow related to cosmic -ray ac- 
celeration, possibly due to an increased role of plasma 
waves (e.g. Rakowski et al., 2008). However, this different 
equilibration behavior between slow and f ast shocks may 
also have another origin: as pointed out bv' Bvkov (2004), 
electrons have thermal velocities larger than the shock ve- 
locity for Mach numbers M < yjm^jmt = 43 . The veloc- 
ity with which the electrons are therefore scattered behind 
the shock is largely determined by the pre-shock tempera- 
ture rather than by the shock velocity. This may lead to ad- 
ditional electron heating mechanisms. A Mach number of 
^ 43 corresponds, for a typical sound speed of 10 kmi^^ 
to a shock velocity of 400 km ' . Interestingly, low Mach 
numbers may also occur for high shocks speeds in the 
presence of cosmic rays. This happens when the cosmic- 
ray precursor significantly preheats the incoming plasma. 
The different equilibration properties of the TeV source 
RCW 86 (Helder et al., 2011) and DEM L71 (Rakowski 
et al.. l2003i) may be explained by the more prominent role 
of cosmic -ray acceleration in RCW 86. This is suggestive 
of an intimate relation between cosmic-ray acceleration 
and temperature equilibration at the shock. 

In the future, with the availability of microcalorimeter 
spectrometers on Astro-H and IXO, it will be possible to 
routinely obtain line broadening measurements from re- 



gions close to SNR shock fronts. The advantage that X-ray 
measurements have over Ha measurements is that no neu- 
tral hydrogen is needed upstream of the shock. Right now 
the sample of young SNRs with thermal Doppler broad- 
ening measurements is skewed toward Type la SNRs, as 
they seem to be more likely to have Ha emission from the 
shock region (Sect. [HJ. The problem with this is that the 
presence of a large fraction of neutral atoms may dampen 
Alfven waves. This may result in different shock-heating 
properti es and may also que nch efficient cosmic -ray accel- 
eration dDrury et al. . Il996h . X-ray measurements of ion 
temperature, therefore, provides a more unbiased probe of 
shock-heating properties. The disadvantage is that on av- 
erage the X-ray emission comes from further behind the 
shock front than Ha emission. 

High resolution X-ray spectroscopy may also be used 
to gather information on the post-shock equilibration 
process itself, because different plasma ionization ages 
should also correspond to different degrees of equilibra- 
tion (i.e. full/partial non-equilibration and equilibration, 
see Sect. 15.2b . 
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Figure 37: Chandra image of SN 1006 based on a mo- 
saic of several observations. The three colors correspond 
to the 0.5-0.6 keV band (O Vll He-a, red), 1- 1 .7 ke V (con- 
tinuum with some Mg XI emission, green), and 2-5 keV 
(continuum dominated, blue). SN 1006 was the first SNR 
for which imaging spectroscopy with the ASCA showed a 
clear spatial separation of continuum emission, identified 
as synchrotron emission mostly from near the rims of the 
SNR, and thermal emission (dominated by O VII), which 



dominates the emission from the interior ( Kovama et al 



1995) 



11.2 X-ray synchrotron emitting filaments 

Although SNRs are thought to be the main sources of 
Galactic cosmic rays (Sect. 15.5b . until 1995 the obser- 
vational evidence for this consisted mostly of the detec- 
tion of radio synchrotron radiation from SNRs. The elec- 
trons responsible for radio emission have energies in the 
GeV range. So radio synchrotron radiation informs us that 
SNRs can accelerate electrons, but does not proof that 
they can accelerate particles to energies comparable to the 
"knee" in the cosmic ray spectrum (^ 3 x lO'^ eV). 

In 1995 it was for the first time convincingly demon- 
strated tha t a SNR, SN 1006, e mitted X-ray synchrotron 
radiation / Kovama" et al.L [l995h . That the X -ray spec- 
trum of this source was unusual was already clear from 
observations with the solid-sta te spectrometer on board 
the Einstein X-ray observatory (iBecker et al. . Il980h . The 
spectrum appeared featureless, and was completely dif- 
ferent from the line-rich X-ray spectra of a young SNRs 
like Tycho or Cas A. This was interpreted by Revnolds 
and Chevalier tl98lh as the result of synchrotron emis- 
sion fro m shock accelerated electrons. However. Hamil- 
ton et al. (Il986h proposed an ingenious interpretation in 
which most of the emission was coming from carbon- 
rich, shocked ejecta, which would provide a large ther- 
mal bremsstrahlung emissivity (Eq. l43l l without much 
line radiation in the 0.5-7 keV range. The case for X-ray 



synchrotron radiation by iKovama et al.l (Il995h was set- 
tled thanks to the spectral-imaging capabilities of ASCA, 
which established that thermal radiation was being emit- 
ted from all over the face of SN 1006, but that the ther- 
mal emission was overwhelmed by non-thermal contin- 
uum emission coming from two regions close to the for- 
ward shock (Fig. [37T i. This had important consequences 
for the idea that SNRs can accelerate particles, as X-ray 
synchrotron emission implies that electrons are acceler- 
ated up to energies of 10 — 100 TeV This requires rather 
small diffusion coefficients (see Eq.|52]for the loss-limited 
case). If indeed maximum electron energies are limited by 
radiative losses, then ions can in principle be accelerated 
to even higher energies, possibly even beyond 10^^ eV. 
But the X-ray synchrotron radiation does not proof that 
this is the case. 

Around the same time that X-ray synchrotron radiation 
from SN 1006 was established non-thermal hard X-ray 
emission was det ected from the young SNR Cas A by 
the CGi;0-OSSE dThe et al.Lll996h. KXTE-YiEXTE (Allen 



et al.. 'l997^ and Beppo-SAX-VT)¥ ("F avata et all Il997h 
instruments. This hard X-ray emission is Ukely caused 
by synchrotron radiation as wel l, although o t her ex plana- 
tions have been proposed (e.g. Allen et al., '1997'; Vink 
et al ., 1999; Laming, 2001 a; Vink, 2008b^Helder and 
Vink, IIOOST 

Moreover, ASCA observations also estab- 
lished that two newly discove red SNRs RX J1713 . 7-3946 
(G347.3-0.5, Pfeffermann and AschenbachLll996h and RX 
J0852.0-4622 (G266.2- 1 .2. fAschenbachLI 1998h were simi- 
lar to SN 1006 in that their X-ray emi ssion was completely 
dominated by synchrotro n emission (IKovama et al.Lll997l ; 
ISlane et all[l999ll200lh . In fact, no thermal X-ray emis- 
sion from these SNRs has vet been detected (e.g. Hiraga 
et al. , 2005; Acero et al. 2009; Slane et al. 2001 ; Pannuti 
et al. j2010al) ! 

The spatial resolution of Chandra helped to establish 
that not only rather large SNRs are source of X-ray syn- 
chrotron radiation (SN 1006, RX J17 13.7-3946, and RX 
J0852.0-4622 are all larger than 30'). Also Cas A (Got- 
thelfet al, 2001 ) Tych o (Hwang et al, 2002) and Kepler 
dRevnolds et al. . l2007h have regions close to the shock 
front that are continuum dominated. These regions are, 
however, very thin (1-2"). It was soon shown that the nar- 
rowness of these filaments must be the result of relatively 
large magnetic fields (^ 100 — 600 AlG,'Vink and Laming', 
2003; Berezhko et al.ll2003; Bamba et al, 2004, 2005; 



fVolk et alLl2005UBalletnB06l; |Parizot et~aril2006l). 

Vink and Laming ("2003") determined a high magnetic 
field for Cas A by assuming the width of the filament is de- 
termined by the radiative loss time of X-ray synchrotron 
emitting electrons. As the plasma is advected away from 
the shock front, a radiative loss time also corresponds to a 
physical length scale. The plasma velocity with respect 
to the shock front is given by V2 — Vs/Z (cf. Eq. fT2l) . 
whereas the synchrotron loss time is given by Eq.|49] The 
synchrotron loss time, therefore, corresponds to an advec- 
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tion length scale of 



^adv — 



X 



(61) 



Since Tsyn depends on both the magnetic field and electron 



energy (Tsyn °= B^^E^^), one needs additional information 
to determine B and E separately. For that one can use the 
dependence of the typical photon energy on the magnetic- 
field strength and electron energy (oc E^B, Eq. |47j. This 
is graphically depicted in Fig. |38] showing that for Cas 
A B ^ 100 jtiG near the shock front. This method is not 
the only way to determine the magnetic field. One can 
also assume that the width of the synchrotron emitting fila- 
ments are equal to the diffusion length scale (Zdiff °^ B^^E, 
Eq. l42b appropriat e for the energ i es of the synchrotron 



emitting electrons (IBerezhko et al. 



Both methods give roughly the same magnetic-field es- 
timates ( Ballet . 2006). In principle, the X-ray synchrotron 
rims should be always wider than, or comparable to, the 
diffusion length scale, as this is the width of the region 
from which electrons are actively being accelerated (indi- 
cated by the shaded region in Fig. 138b). But the diffusion 
length scale is always comparable to the advection length 
scale f or electrons with energies close to th e maximum 
energy (IP arizot et al. , 20061 : Vink et al. , 2006). This is eas- 
ily demonstrated using the basic condition for determining 
the maximum electron energy, i.e. the acceleration time 
equals the synchrotron cooling time, Tacc ~ fsyn (Eq.lSlb. 
The acceleration time scale is within an order of magni- 
tude given by Tacc ~ D/v2 (c.f. Eq.l4TI). Using Eq.|6T]one 
can rewrite the condition for reaching the maximum elec- 
tron energy, Tacc ~ '^syn? ill "^acc ~ 4dv/v2- Eq. l42l shows 
that Zdiff ~ D/v2 « 1^accV2' Hence, the condition Tacc ~ Tsyn 
also implies Zjiff ~ kdv 

By combining Eq. |49] and Eq. |5T] one obtains a rela- 
tion between the typical width of the X-ray synchrotron 
emitting region, assumed to be /adv, and the magnetic - 
field strength, under the assumption that we observe syn- 
chr qtron emission frorn electrons near th e cut-off energy 
(c.f. lParizot etan.l2006l:lvink et al.ll2006h : 



B7 «26 



adv 



Ox lO'^cm 



-2/3 



J?'/'(Z4-^) '^%G, (62) 



with B2 the average magnetic field in the shocked plasma, 
X4 the total shock compression ratio in units of 4, and 
T] as defined in Eq. [39] Note that this magnetic field es- 
timate is independent of the shock velocity and depends 
only weakly on the shock compression ratio and 77. Since 
the actual width is determined by a combination of diffu- 
sion and advection the actual observed width is larger than 
either /adv or Zdiff, and one can take Zobs ~ %/24dv 

It is worth mentioning that not all X-ray synchrotron 
radiation may come from loss-limited spectra. For large 
SNRs like RX J1713.7-3946, with a rather broad X-ray 
synchrotron emitting shell, one may wonder if the mag- 
netic field is not so large that the X-ray synchrotron spec- 
trum is age limit ed. Also for the youngest known Galactic 
SNR, G1.9H-0.3 jRevnolds et al.[l2008l: iBorkowski et"ar 



2010h . it is not clear whether the X-ray synchrotron emis- 
sion is age or loss limited. One possible way to test this 
is too monitor the non-thermal X-ray emission for some 
time. For example, recently it has been shown that the 
X-ray synchrotron brightness of Cas A has declined be- 
tween 20 00 and 2010. At the same time the spectrum has 



softened dPatnaude et al.L 1201 1) . This can be understood 



assuming a loss-limited spectrum, as the exponential cut- 
off energy depends solely on (Eq.l52]l. Since the shock 
is decelerating the cut-off energy will decrease, thereby 
making the spectral index larger (corresponding to a softer 
spectrum) and the brightness smaller. For an age-limited 
synchrotron spectrum the cut-off photon energy depend 
on the age of the SNR. Thus the cut-off energy will in- 
crease with time. Depending on the magnetic-field evolu- 
tion in the shock region this could i n principle make the 



X-ray synchrotron spectrum harder (ICarlton et al.L 1201 Ih . 



The secular evolution of the X-ray synchrotron emission 
does, however, also depend on other variables, such as the 
number of electrons accelerated as a function of age, but 
the evolution of the exponential cut-off is likely to be the 
dom inant contribution to the spectral evolution ( Patnaude 
et al. I2OI iLlKatsuda et al.Ll2010bh . 



11.3 X-ray based evidence for efficient 
cosmic-ray acceleration 

One of the surprises that came out of the magnetic-field 
determination was that the magnetic fields were typically 
100 — 600 /xG, larger than might be expected, if they were 
merely caused by the compression of the interstellar mag- 
netic field (BisM ~ 5 fiG). Additional evidence for these 
high magnetic fields is that the X-ray synchrotron emitting 
regions seem to vary on time scales of several months to 
years, which may be indicative of the synchrotron cooling 
time scales, and, hence high magnetic fields (Patnaude and 
Fesen, 2007; Uc hivama et al.. 2007; Uchiyama and Aharo- 
nian, 2008)0 The high magnetic fields suggest that some 
magnetic-field amplification mechanism is operating near 
the shocks of young SNRs. 

There are several theories that Unk magnetic-field ampli- 
fication to efficient cosmic-ray acceleration. A particular 
promising mechanism is that discovered by Bell (.20041) . 
in which cosmic rays streaming outward result in an elec- 
tron return current. These electrons are then deflected by 
the magnetic field, thereby creating a very turbulent, am- 
plified magnetic field. This mechanism is expected to give 
ub ~ VsUcr °= pV|, with ub and Ucr respectively the mag- 
netic field and cosmic-ray energy density, and po the am- 
bient plasma density. 

Combining the magnetic-field estimates of several 
young SNRs shows that there is indeed a clear trend that 
SNRs evolving in low density environments seem to have 
lower magnetic fields than those evolving in high density 
media (Bamb a et al., 2005; Volk et al., 2005; Ballet, 2006 ; 
Vinki l2008al) . This also explains the trend that the wider 



However, an alternative explanation ma y be varying magnetic fields 
due to passages of large scale plasma waves jBvkov et alll2008bh . 
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Figure 38: Left: The maximum cosmic -ray electron energy versus magneti c-field strength for the reg ion just downstream 
of Cas A's shock front, as determined from the thickness of the filaments (IVink and LamingLl2003r) . The shaded area is 
exclu ded, because the filament width cannot be s maller than the minimum possible diffusion length (c.f. lVink and Lamina . 
2003h . (These figures were published before in (l\any,|2006|). ) Right: The power slope of the X-ray continuum in the 4-6 
keV band based on the 1 MS Chandra observation of Cas A (IHelder and V ink. 2008). The color scale indicates the slope 
difference with respect to the mean slope of F = 3 .2. Brighter color mean harder spectra, and signifiy a larger contribution 
of synchrotron radiation. Note that there seem to be two general regions of synchrotron radiation, the outer edge and the 
inner region near the reverse shock. A high resolution Chandra X-ray image of Cas A is shown in Fig. |26] where the 
X-ray synchrotron can be seen as a faint blue wisp at the edge of the SNR. 
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Figure 39: The dependence of the post-shock magnetic 
field energy density B^/(8;r) on the shock velocity V,. The 
dashed line shows a V} dependency and the dotted line a 
de pendency. The input val ues can be found in ( Vink , 
120061), except for SN1993J (iFransson and Biornssonl 
[1998) and Kepler's SNR, for which the shoc k velocity in 
the southw est was taken f rom (Vink, 2008c') and the den- 
sity from (iCassam-Chenai et al..,2004,) . 



X-ray synchrotron emitting regions, hence lower magnetic 
fields, are found in larger SNRs: these SNRs have a com- 
bination of high shock velocity and large radiu s, because 
they evolve in a low density region (lUeno et all 12006) . 



The relation between magnetic-field energy density and 
shock velocity is less clear The reason is that for X-ray 
synchrotron emitting SNRs the dynamic rang e in veloc- 
ity is relatively small 2000 - 6000 km.y~h. lV51k et al. 
12005) find that °^ PoV^ , but they used an outdated mea- 
surement of t he shock velocity of Cas A 2500 kmi^'). 
VinkI (12008a') instead favors a dependency of / (87r) °= 



P0V5 , but noted that the dynamic range makes the depen- 
dency on Vs uncertain. Including the magnetic-field deter- 
minati on for SN 1993J (Vs » 2 0,000 kmi_^nd_B« 
64 G, IFransson and Biornssonl |i998; Tati schefj . l2009h 
indicates that B^ °^ seems more likely (Fig. [39T l. But 
it should be noted that it is not quite clear whether the 
magnetic-field amplification mechanism in SN1993J has 
a similar origin as in young SNRs. 

There has been some concern that amplified, turbulent 
magnetic fields dampen rapidly, even to the point that the 
width of the narrow X-ray synchrotron filaments in young 
SNRs may not be determined by the advection or diffusion 
length scale, but by the magnetic-field damping time scale 
dPohl et all l2005h . However, this would affect the radio 
and X-ray synchrotron profile similarly, and this does not 
appear to be th e case (Cassam-Chenai et al., 2007; Got- 
thelfet al.r EoOlb . 

The X-ray synchrotron emission from SN 1006 has an- 
other characteristic that may be of importance for a better 
understanding of cosmic-ray acceleration: it appears that 
the synchrotron radiation is coming from just two regions 
of the SNR, i.e. there seems to be very little synchrotron 
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emission from the center of the SNR (IWillingale et al 



10 



1996l:lRothenflug et aLlEoollCassam-Chenai et allboosl) . 
A likely explanation is that this is due to initial the mag- 
netic field configuration. Given the inferred polar geome- 
try, the implication is that X-ray synchrotron radiation is 
favored by a magnetic field that is parallel to the shock nor- 
mal. On the other hand, the presence of X-ray synchrotron 
radiation requires a turbulent magnetic field. This seems 
a contradiction, but the two conditions can be combined 
by assuming that efficient cosmic-ray acceleration favors 
a parallel magnetic field to start, but on ce it has started i t 
creates its own turbulent magnetic field (' Volketal.lEooi) . 
SN 1006 may be special in that its background magnetic 
field was much more regular than that of other X-ray syn- 
chrotron emitting SNRs. 

One question that may be raised is whether the magnetic 
fields in young SNRs are caused by amplification, or that 
even creation of magnetic fields is possible. There is at 
least some indica t ion for magnetic-field creation, because 
Helder and VinkI (l2008h showed, on the basis of a spec- 
tral analysis of Chandra X-ray data (Fig. [38t)). that most 
of the X-ray synchrotron emitting filaments are not asso- 
ciated w ith the forward shock, but with th e reverse shock 
(see also lUchivama and Aharonian , boos '). Since these fil- 
aments are also narrow, the local magnetic fields must be 
of the order of 100 - 500 jtiG. A problem with this is that 
the ejecta, due to their large expansion must have very low 
magnetic fields. So this suggests that Bell's mechanism, or 
any other ampUfication mechanism, even operates for very 
low seed magnetic fields, or perhaps does not even require 
a seed magnetic field. Alternatively, cosmic rays diffus- 
ing from the forward shock to the reverse shock may have 
provided a seed population of cosmic rays that resulted 
in magnetic-field arn plification through Bell's mechanism 
( Schure et al. , 2010l) . Note that there are also suggestions 
that t he reverse shock region^ of RCW 86 (Rho et al., 
20021) and RX J1713.7-3946 ( Zirakashvili and Aharonian. 
2010l) are sources of X-ray synchrotron radiation. 

Irrespective of the exact magnetic-field dependence on 
shock velocity, the evidence for magnetic-field amplifica- 
tion in young SNRs has important consequences. The pres- 
ence of amplified magnetic fields imply large cosmic -ray 
energy densities. This makes it likely that cosmic rays are 
indeed being efficiently accelerated, i.e. a large fraction of 
the shock energy may be transferred to cosmic rays. This 
in turn suggests that the post-shock temperatures are ex- 
pected to be low and the shock compression ratios should 
be larger than four (Sect. lSTt . 

As argued in Sect. 111.11 the lack of very hot electron 
temperatures in young SNRs may be caused either by a 
lack of electron-proton temperature equilibration, or by ef- 
ficient acceleration. Proton temperatures are less affected 
by equilibration effects, and can be determined fr o m Ha 
line widths. Using this technique iHelder et al. (l2009l) 



showed that in the northeastern part of RCW 86 the pro- 
ton temperature is w 2 keV, much lower than the ex- 
pected kT > \Q keV. They derived from this a post-shock 
cosmic -ray pressure contribution of > 50 %. Note that 
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Figure 40: Radio to X-ray s ynchrotron spect rum of the 
northeastern part of RCW 86 ( IVinketal.Ll2006h . The mod- 
els show that a curved synchrotron spectrum is needed. A 
simple synchrotron spectrum with a spectral radio index 
a = 0.6 (corresponding to an electron particle index of 
q = 2.2, Sect. 111. 2T i and an exponential cut-off in electron 
energy cannot simultaneuously fit both the radio and X- 
ray data points. 



RCW 86 is a TeV source dAharonian et all l2009l) and 
that t he northeaste r n part emits X-ray synchrotron radia- 
tion dBa mba et all I2OO0I: Ifiorkowski et all l2001bl: Vink 
et al., 2006|). A similar study of the rapidly expanding 
LMC SNR B0509-67.5 (Sect. lU indicates that there in 
the southwest the post -shock pressure is f or > 20% pro- 
vided by cosmic rays dHelder et alll2010h . This method 
can also be employed using X-ray spectroscopy once high 
resolution, non-dispersive, spectrometers will be available 

Another effect of efficient cosmic -ray acceleration is 
high shock compression ratios (Sect. 15.1b . There is in- 
deed evidence for high shock compression ratios in Ty- 



cho (Warren et al., 2005: Cassam-Chen ai' et all 2007 ) and 
SN 1006 ( Cassam-Chenai et al.„.2008.) . Although the com- 
pression ratio cannot be measured directly. X-ray imaging 
spectroscopy of both SNRs shows that the ejecta that are 
located very close to the forward shock, which is an ob- 
servational con sequence of a high compression ratio ( De- 



2000I) . 



courchelle et al., 

Efficient acceleration should also lead to a curvature 
of the relativistic electron spectrum. The reason is that 
particles of different energies experience different com- 
pression ratios during the acceleration process, depending 
on whether they experience the overall shock compres- 
sion ratio, the gas-shock (sub-shock) compression ratio, 
or something in between (see Sect. 15.51 ). Radio spectra of 
a few SNRs do indeed show evid ence for spectral hard- 
ening with increasing frequencies ( Revnolds and EUisonL 
19921). For X-ray spectroscopy the consequence of spec- 
tral curvature is that the X-ray synchrotron brightness can- 
not be simply estimated from an extrapolation of the radio 
synchrotron spectrum, but it should be brighter than ex- 
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pected. This has indeed been demo nstrated to be the c ase 
for the northeastern rim of RCW 8 6 (fvink et al.Ll2()()6l see 
Fig.gOli and SN 1006 dAUen et al.ll2008h . 



11.4 The relation between X-ray syn- 
chrotron and 7-ray emission 

The field of GeV and TeV 7-ray astronomy has made 
tremendous progress over the last decade, thanks to tech- 
nical progre ss in building sensitive TeV C herenkov tele- 
scopes (see iHinton and Hofman for a review) 
and, very recently the launch of the NASA Fermi satel- 
lite, which is sensitive to ^ 0.1 — 100 GeV 7-rays. To- 
gether these telescopes have greatly expanded the number 
of 7-ray sources, and, over the last decade have resulted in 
the first firm detections of continuum 7-ray emission from 
shell-type SNRs. The first such detection in the T eV band 
was Cas A with HEGRA dAharonian et alll200il) . but the 
list of TeV detected telescopes has greatly expanded since 
then@ 

The 0. 1 — 100 GeV band was sampled by 7-ray satellites 
since the 1970ies, but the first evidence for associations 
between GeV 7-ray sources and SNRs were found with 
CGRO-EGRET, and co ncerned mature SNRs interacting 
with molecular clouds (Espos ito et al. . Il996h . Fermi has 
confirmed these early associations, and greatly enhanced 
the number of 7-ray sources among mature SNRs interact- 
ing with molecular clouds (Sect. 110.3b . In addition, also 
a number of younger SNRs have now been detected. The 
first catalog of F ermi listed 4 1 source s as possibly associ- 



2010al) and this Ust is still 



ated with SNRs dAbdo et al. 
expanding 

E 

The great importance of GeV and TeV 7-ray astron- 
omy is that it can give us a direct view of the acceler- 
ated cosmic-ray nucleQ which make up 99% of the cos- 
mic rays observed on earth. These nuclei produce 7-ray 
emission if they collide with atomic nuclei, thereby cre- 
ating, among others, neutral pions, which decay in two 
7-ray photons. However, two other important 7-ray radia- 
tion processes originate from relativistic electrons: interac- 
tions with background photons result in inverse Compton 
up-scattering, whereas interactions with ions in the SNR 
result in bremsstrahlung. X-ray observations are crucial 
for disentanghng all these contributions. The thermal X- 
ray emission is important for determining the local plasma 
density, which determines the 7-ray non-thermal brems- 
strahlung and pion-decay contributions. The X-ray syn- 
chrotron emission can be used to determine the local mag- 
netic field. The combination of magnetic field and radio 
and X-ray synchrotron flux can then be used to infer the 
inverse Compton contribution to the 7-ray emission. 



At the moment of writing 1 3 shell-type SNRs are now claimed to be 
TeV sources, see http://tevcat.uchicago.edu 

'^Note that the 7-ray may not necessarily come from the SNR shell, 
but in some cases from an embedded pulsar wind nebula. 

^^'^Protons and other ions. These are often named "hadronic cosmic 
rays", as opposed to leptonic cosmic rays, which are mostly electrons 
and positrons. 



Although X-ray observations seem to indicate that 
cosmic-ray acceleration is efficient in young SNRs, the 
recent observations of SNRs in GeV and TeV 7-rays pro- 
vide a more ambiguous story. Although many SNRs have 
been det ected in TeV 7-rays w ith TeV telescopes like 
HEGRA (Aharonian et al.,'200l'), HESS (e.2. Aharonian 
et al. . 2004). MAGIC (Albert et al., 2007b) and Veritas 



(e.g. lAcciari et al.L 1201 Ih and the Fermi satellite, the 7- 



ray observations do not provide clear evidence that cosmic 
rays are being accelerated efficiently, nor that they SNRs 
accelerate ions to beyond "the knee". 

There has been a strong debate about the nature of TeV 
emission X-ray synchrotron radiation dominated SNRs. A 
central role in this debate concerns the brightest TeV emit- 
ting SNR, RX J1713.7-3 946. For example Berezhko and 
Volk (2010) advocates a dominant pion-decay contribu- 
tion to the TeV y-rav emission, whereas Katz and Wax- 
man ( 2008) and'Eflison et alJ tOl ^ argue that the 7-ray 
emission is dominated by inverse Compton scattering. The 
lack of thermal X-ray emission is an important ingredient 
in this debate, as to some it indicates low densities, and 
hence low pion-decay l uminosities, whereas others (no- 
tably Drurv et al. , 20091) have argued that it could be the 
result of extremely low temperatures caused by efficient 
cosmic-ray acceleration (see Sect. 15. 11 1. An argument in 
favor of an electron origin for the 7-ray emission of RX 
Jl 7 13.7-3946 is the close resemblance between the TeV 
7-ray and X-ray brightness distribution, which suggests 
that both originat e from the same po pulation of particles: 
electrons (F ig.l4n lAceroetalll200 9'). 

However, there are two arguments in favor of pion de- 
cay. Namely the relatively hard spectrum in TeV, which ex- 
tends up to 100 TeV, an d the inferr ed high magnetic field 
in the SNR, > mG (Uchi vama et al .. 2007)- Higher mag- 
netic fields imply lower relativistic electron densities for 
a given X-ray synchrotron flux. Hence, the electron con- 
tribution to the TeV flux is expected to be low, and pion 
decay becomes the most likely dominant source of TeV 7- 
rays. The high magnetic-field estimates are based on the 
the observation of rapid X-ray flux changes and compar- 
ing them to the synchrotron loss time scales (see Eq.|49]l. 
But the fluctuations may also originate from fluctuations 
of due to long wavelength Alfven waves (Bvko v e t all 



2008bh . Another argument in favor of relatively high mag- 



netic fields are the small widths of the X-ray synchrotron 
rims (Sect. [TL2l i. But for this measurement the problem is 
what region one should measure the width of: the whole 
shell or just a bright filament? For example, one can take 
the overall shell size of RX J1713.7-3946 to be the size of 
the filament (since for this SNR the whole shell seems to 
emit synchrotron emission), or use some filamentary sub- 
structures. This problem is illustrated in (Fig. |42] |. If one 
takes the whole shell region, at best one comes up with a 
length scale of 100" corresponding to 25 — 30 jiG. 

The most serious challenge to the pion-decay model 
for R X J17 13.7-3946 comes from recent observa tions by 
Fermi dAbdo and Fermi LAT Collaborationll201 1 "), which 
show that the shape of the broad band 7-ray spectrum 
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Figure41: Left: 7- ray excess count map of RX J 17 13.7-3946 as observed by HESS. Overlayed are 1-3 keV X-ray (A5CA) 
intensity contours (lAharonian et all 120070 . Right: A co mparison of the X -ray flux (1-10 ksW,XMM -Newton) and HESS 
y-ray integrated flux (1-10 TeV band, Figure taken from lAcero et al.L 120091) . 
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Figure 42: X-ray emission profile of the southeastern part 
of RX J1713.7-3846 (extracted by the author fro m the 
XMM-Newton image reduced bv lAcero et alll2009l) . It il- 
lustrates the difficulty in determining for this SNR what 
the width of the X-ray synchrotron emitting region is. The 
red line indicates a simple model of a projected sphere 
with a uniform emissivity with a shell width of 1 16". For 
a distance of 1 kpc this corresponds, according to Eq. |62] 
to 25 liG. 



is more consistent with inverse Compton scattering than 
pion decay. The implication of this result is surprising, be- 
cause RX J17 13.7-3946 was often used as a key object in 
the argument that SNRs can accelerate ions to very high 
energies. 



An important result of Fermi dAbdo et al.Ll2010d) is the 



detection of GeV 7-ray emission from Cas A, which indi- 
cates that the total cosmic-ray energy content is less than 



4% of the total explosion energy. This is disappointing if 
one thinks that young SNRs should be able to accelerate 
cosmic rays efficiently, and convert > 10% of the explo- 
sion energy to cosmic rays (Sect. 15.51 ). The case for pion 
decay as origin for the 7-ray emission is stronger for Ty- 
cho's SN R, given the recent de tection with the Veritas TeV 
telescope dAcciari et al.M201 ih . 

The apparent lack of evidence for a high cosmic ray ac- 
celeration efficiency for Cas A is also somewhat surpris- 
ing, given the high magnetic field of Cas A, as argued in 
Sect. II 1.21 which is suggestive of efficient cosmic -ray ac- 
celeration. It could be that a large part of the energy in 
cosmic rays has already escaped the SNR. Since Cas A 
evolves inside the wind of its progenitor, the background 
plasma that the escaped cosmic rays are in likely has a 
lower density than the SNR shell. 

Indeed, several theories about non-linear cosmic ray 
accelera tion point toward t he importance of cosmic ray 
escape dBlasi e t al., 2005; Vladimiro v et all I2O O8': Re- 
ville et al., l2009u ,Vink et al.. .2010') . Future very sensi- 
tive T eV telescopes like the C herenkov Telescope Array 
(CTA, CTA ConsortiumL 2010l) could in principle be used 
to search for "cosmic-ray haloes" around young SNRs. 
Hopefully, the improvement in sensitivity in TeV astron- 
omy can be met with a similar improvement in sensitivity 
and spectral resolution in X-ray spectroscopy, which will 
help to detect thermal X-ray emission and provide temper- 
atures for X-ray synchrotron emitting plasmas. 

Interestingly enough, the case for a pion-decay origin of 
the 7-ray emission is most unambiguous for some of the 
mature SNRs that are interacting with or are near mo lec- 



ular clouds. An example is W44 (.Abdo et alil2010eL see 



also Sect. 110.31 1. for which it is argued that inverse Comp- 
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ton scattering, or bremsstrahlung do not provide an ade- 
quate description of the spectrum. But the y-ray spectrum 
of this source steepens around 10 GeV, which is consis- 
tent with the idea that cosmic rays with energies above 
lO'^ eV are present. This break in the spectrum could 
be evidence that cosmic rays with larger energies have 
already escaped the SNR in the past. This is also con- 
sistent with several SNRs that have 7-ray sources nearby, 
which could be caused by cosmic rays that escaped from 
the SNR and are now interacting with the high density 
medium of molecular clouds. The GeV/TeV sources as- 
sociated wi th SNR W2 8 is an exam ple for such an associ- 
ation tAharonian eTallliboSk.Abdo et aLlHoiOb.) !^ 



12 Concluding remarks 

Our understanding of SNRs has greatly increased over the 
last decade due to the X-ray observatories Chandra, XMM- 
Newton, and Suzaku. The second half of this review es- 
sentially summarized the advances in SNR research made 
with these observatories. 1 think that a number of these re- 
sults deserve more recognition outside the field of SNR 
research. I will summarize some of them here. 

X-ray spectroscopy has proven to be a promising tool to 
separate young Type la (Sect.[8]i from core collapse SNRs 
(Sect. 19. 11 1. X-ray emission for Type la SNRs has given 
us important clues toward the nature of the progenitor sys- 
tems: the progenitors did not produce fast winds, and de- 
layed detonation models seem to explain the abundance 
patterns well. Type la SNRs are also well stratified, unlike 
core collapse SNRs. 

An important question concerning core collapse super- 
novae is what drives the explosion (neutrino absorption, 
or other mechanisms?) and what types of massive stars 
produce neutron stars and what types make black holes. 
X-ray observations show that the core collapse process is 
chaotic, with pure metal plasma obtaining high velocities 
in seemingly random directions. But for one SNR, Cas A, 
there is evidence for an ejecta jet/counter jet. These are in 
Si-group material, not in Fe-rich material. With future ob- 
servations we may learn more about rare elements, which 
may indicate the neutron content of the inner most ejecta. 
This may help us to identify the role of neutrino absorp- 
tion in the explosion. 

As for the question what determines the creation of a 
black hole. There is at least one SNR, G292.0H-1.8, for 
which there is an indication that it produced a neutron 
star, but had nevertheless a progenitor mass above 30 Mq 
(Sect. |93T l. It is important to confirm this results with fur- 
ther X-ray studies, as a popular hypothesis is that stars 
above 25 Mq leave behind a black hole. 

Another uncertainty in the theory of core collapse su- 
pernovae is how much iron they produce. The iron yield 
is greatly affected by the details of the explosion, such 
as explosion asymmetries and fall back on the neutron 



star There is evidence that the Cygnus Loop (Sect. 110.21) 
contains more iron than predicted by explosion models. 
Whereas for Cas A, iron seems to have obtained a higher 
velocity than predicted by explosion models. 

One surprising result of X-ray imaging spectroscopy 
is that even mature SNRs still have regions with en- 
hanced metal abundances. This is true for the Cygnus 
Loop (Sect. 110.2b and many mixed-morphology SNRs 
(Sect. 110.3b . For the Vela SNR X-ray imaging spec- 
troscopy confirmed that the bullets outside the main shock 
are indeed ejecta "shrapnels" that have enhanced abun- 
dances. 

Finally, X-ray imaging spectroscopy has played an im- 
portant role in increasing our understanding of cosmic- 
ray acceleration. A role it shared with the rapidly evolv- 
ing fields of TeV astronomy with Cherenkov telescopes, 
and GeV astronomy with the Fermi satellite. In particular. 
X-ray detections of narrow width synchrotron filaments 
have provided evidence for magnetic-field amplification 
near SNR shocks (Sect. ll 1.2b . This implies that SNRs are 
efficiently accelerating cosmic rays. However, the overall 
interpretation of both X-ray and 7-ray observations is still 
ambiguous, in particular concerning the total energy in 
cosmic-rays contained in young SNRs (Sect. [TTT4b . 

As is clear from the above, our understanding of SNRs 
has greatly benefitted from the coming of age of X-ray 
imaging spectroscopy. The impact of high resolution X- 
ray spectroscopy has been more modest, as Chandra and 
XMM-Newton have grating spectrometers, which are not 
ideal for observations of extended sources as SNRs. Nev- 
ertheless, for bright sources of small angular extent impor- 
tant results have been obtained. 

In the near future high resolution X-ray spectroscopy 
of SNRs will likely become much more important as mi- 
crocalorimeters will become operational. These will have 
a spectral resolution of Afi w 5 eV at 6 keV {R w 1200) 
and also have some limited spatial resolution. The first op- 
portunity to use them will start after the launch of Astro- 
After Astro-H the next step forward will hopefully be 
ATHENA, a mission concept for an European X-ray obser- 
vatory that will offer a large effective area (> 1 m^ at 1 
keV) and a spatial resolution of < 10", which should be 
compared to 1' for Astro-H. 

The tantalizing glimpse of high resolution X-ray spec- 
troscopy of SNRs offered by Chandra and XMM-Newton, 
shows that SNR research beyond 2014 will enter a new 
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